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The 


Quasi-Steady Air Forces for Use in 


Low-Frequency Stability Calculations 


MARTIN GOLAND* 
Midwest Research Institute 


SUMMARY 


Simple, approximate expressions are derived for the air forces 
acting on an airfoil executing oscillatory motions with arbitrary 
damping. These expressions are in forms convenient for use in 
nalyses of aircraft stability and response. The derivations are 
of a simple character, consisting merely of power series approxi- 
mations for the known results from unsteady aerodynamic theory. 
It is shown that the present approach enables the evaluation of 
ircraft stability derivatives with better theoretical accuracy than 
do the methods now conventionally employed, while not appre- 
tiably increasing the labor associated with their use. In the pres- 
ent paper, detailed attenticn is concentrated on the determina- 
ion of the stability derivatives for low-frequency motions; pre- 
liminary consideration is also given to methods for dealing with 

e higher frequency range. 


INTRODUCTION 


N A RECENT STABILITY and response investigation of a 

special nature conducted by the author,f it was 
ound that the value of the stability derivative c,,,5, the 
otal mement on the craft due to the elevator rotational 
elocity, entered into the analysis in a most critical 
ashion. In particular, if this derivative were assigned 
A negative value (presuming the usual N.A.C.A. sign 
onventions, to be given in detail later), the aircraft 
isplayed a static divergence that is contrary to physi- 
al experience. Furthermore, the Cmns calculated for 


€ system by means of the techniques in current use, 

Received January 19, 1950. Revised and received April 3, 
950. 
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such as those described by Jones and Cohen! and White 
and Klampe,’ does show such a negative value. This, 
of course, suggested that the conventional methods 
for calculating the c,,,; derivative might require some 
modification. 

It readily became apparent that this is actually the 
case. This can be seen as follows: The total aircraft 
moment arising from motion of the elevator is, to a very 
good approximation, equal to the negative of the prod- 
uct of the tail length of the craft and the tail force due 
to elevator motion. The tail lift coefficient arising 
from elevator motions is given by White and Klampe? 


in their Eq. (21) as 


Two 
m 
T 
where m is the slope of the tail lift coefficient vs. angle 
of attack, and the various 7 terms are defined by 
Theodorsen.* Theodorsen points out in his Table | 
that both 7\) and 7), are positive, while the value of 
T; is negative. Hence, for elevator oscillations of the 
type e”, it follows from Eq. (1) that both the real and 
imaginary parts of the tail lift due to elevator motion 
are positive in value. This amounts to the statement 
that both c,,; and C»,»s5 are always negative. 


b dé 
1 dt 


Th 


2r 


(1) 


9 


re 


3 + = Ty; +m 


The lift coefficient given by Eq. (1) is presumably an 
approximation of the coefficient calculated exactly by 
Theodorsen.* However, if Theodorsen’s Eq. (XVIII) 
is used and if the lift coefficient due to elevator motion 
is computed for values wb/v in the range from 0 to 0.1, 
which is of interest in stability problems, it is found that 
the real and imaginary parts of the coefficient have 
opposite signs, the real part being positive but the imag- 
inary part having a negative value. This is confirmed 
by a glance at the table of oscillating airfoil coefficients 
published by Smilg and Wasserman;‘ at v/bw values of 
10 and 16.67, the real and imaginary parts of the co- 
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efficient L, are of opposite sign for all elevator con- 
figurations. 

This contradiction between the signs of the imagin- 
ary parts of the lift coefficients calculated by the exact 
procedure and from Eq. (1) leads to the conclusion that 
Eq. (1) requires modification. Hence, the methods 
used to calculate ¢,,,, in such sources as references | 
and 2 are also in need of alteration. 

The source of the discrepancy in Eq. (1) will now be 
explored, and a more suitable expression will be derived. 
It will further be shown that the procedures in current 
use for calculation of other stability derivatives should 
also be modified, and appropriate forms for these cases 
will also be given. 


THE DERIVATIVES Cy, AND Czp5 FOR LOW-FREQUENCY 
OSCILLATIONS 


It is instructive to continue consideration of the air- 
foil lift due to control surface motion on a more rigor- 
ous basis. Consideration is restricted here to con- 
figurations where the elevator has no aerodynamic 
balance. For an airfoil of semichord 6 and an in- 
finite aspect ratio, Theodorsen,* in his Eq. (XVIII), 
finds the lift per unit span due to oscillation of the ele- 
vator to be 


dé d*6 
Mg= —pb?(vT +6b7T%—)+ 
L, ” ( ‘at tiie =) 


ee) 

(2) 
2m dt 
where 6 represents the elevator deflection, positive for 
a depressed trailing edge; all other notations are taken 
directly from Theodorsen. Note that the elevator 
configuration determines the values of the various 7 
parameters. By conversion from physical units of 
time to the dimensionless time 7, where 


F Ti 
2a pvbC (« =% +b 
us 


r = (v/b)t 


Eq. (2) becomes 


To 
6+ 


T 2a 


L; ai ps| —Tid — 7,6 + 2rC ( 


where a dot denotes differentiation with respect to r. 
The term 27C represents the slope of the lift coeffi- 
cient vs. angle of attack for the airfoil under the dynamic 
conditions of oscillation; Theodorsen shows that C is, 
in general, a complex number and is a function of the 
dimensionless frequency k = wb/v, where w is the fre- 
quency of oscillation in physical units. 

It should be noted that expressions equivalent to 
Eq. (3), but for control surfaces with aerodynamic 
balance, are given by Theodorsen and Garrick.® This 
permits direct extension of the present analysis to 
cover such cases. 


AERONAUTICAL 


AAC pe 


SCIENCES—OCTOBER, 1959 
It follows from Eq. (3) that the lift coefficient fo, th , 
. ° . . . . - . .s 
airfoil section is the portion of the right-hand Side oj 
the brackets. Fo, ‘ 
elevator oscillation according to e’“’, the lift coefficien: 


the expression contained within 


per unit elevator deflection is 


k* jer — jkT, + 2nC (= + jk | 


T 2r 


The problem now is to arrive at a simplified forn , 
Eq. (4), which is valid for low frequencies of oscillatio, 
In current practice, this ; 


I 


i.e., small values of k. 
done by replacing 27C by its value when k = (0, it the 
being equal to the steady-state lift coefficient vs. angk 
of-attack slope m and by retaining terms in Eq. ({ 
which are of the order k. 
leads to the Eq. (1) statement. 
however, that this approach to the problem does no) 
yield results that are entirely correct. 

The most direct method for approximating Eq. ({ 


It is seen that this procedur 
It will now be showy 


for low values of k would appear to be to expand th: 
expression as a Taylor series about k = 0 and retajy 


terms up tothe orderk. Thus, 


l; ~ (15)o + k (dl; dk) 


where the subscript zero denotes the value whe 
k=0. Toattempt this, form the derivative 
dls an ern ° Ty ac 
= 2kT, —_ 114 + ] 2rC - + 2r x 
dk : ; 2r dk 


(= +H =) , 
T 2r 


The nature of the derivative dC/dk for the small val. 
ues of &, and in particular when k = 0, has been studied 
For an airfoil of in 





by Heilenday and Rauscher.‘ 
finite aspect ratio they find that, for k small and pos 


tive, 
dC dF dG ‘ ‘ 
= j ees sie. 7 
dk dk ~*~ dk 2 
' 
and this shows that as k approaches zero, the imaginary 
part of dC/dk approaches minus infinity. | 
This at once leads to the conclusion that a Taylor serie} 
expansion of |; about k = O is not permissible. Such 


an expansion about a point k > 0 is possible, but for the 
present purpose this procedure is subject to the greal 
disadvantage that the air forces are not returned pre 
cisely to their steady-state values when k = 0.* 

While the difficulty described above can never b 
overcome theoretically, it is possible to approach the 
problem of approximating /; in the following fashion 
If the function G, the imaginary part of C, is approx 
mated by a best-fit straight line in some small, arbi 
trary region of k around zero, then the slope of this lint 
can be used as the value of dG/dk at k = 0. Tht 

* The question of expansion about a point & different from ze! 
is discussed further later in this paper. 
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straight-line fit of G should, as noted above, preferably 
return to the value zero when k = 0. 

For example, Fig. 1 shows a plot of G vs. k for small 
values of the abscissa. Also shown are two straight- 
line fits of this curve, calculated by the method of least 
squares. The line of higher slope on the plot is based 
on a five-point calculation using values of k between 0 
and 0.05 in steps of 0.01. The equation for this line is 


G'00.01)0.05 = —2.922 k (S) 


The second line is fitted by ten points in the region from 
0 to 0.10 in steps of 0.01 and has the equation 


Go 0.01)0.1 = — ?.086 k (9) 


Based on this approach to the problem, the Taylor 
series expansion for l;, accurate to order k, is readily 
shown to be 


Tw iF ; 
vo. eee ( ) a] at x 
T dk / 0 


dG’ To 


2m (10) 
dk wr; ' 


- Ty 
—T, + m—+ 
2x 
A comparison of Eqs. (10) and (1) shows that they dif- 
fer by the last terms in each of the parentheses con- 
tained in the former equation. 
It follows from Eq. (10) that the stability derivatives 
¢,,and ¢.,,;for harmonic oscillation are given by 


‘ Tio | 4 > (Tr) | 1 
| ; = 7 e ( ) 
; Tv 7 ~ dk 0 


tie 1G’ Tw 
Cia, = (-7. +m = + 2 ") (12) 


1 
2r dk 


It is shown in Appendix I that it is valid to use these 
derivatives for the case of either decaying or growing 
oscillatory motions, as well as harmonic motion, by 
permitting the value of the ‘frequency’ k to be com- 
plex. Thus, ¢.; and ¢.,; are stability derivatives that 
entirely conform to the usual connotations. For com- 
plex values of k, Eq. (7) continues to hold, provided the 
real part of k is positive.* 

As an example of the results obtained from Eqs. (11) 
and (12), Fig. 2 shows a plot of these derivatives as a 
function of & in the range from 0 to 0.10. An elevator 
with hinge line at the mid-chord is presumed, and the 
following coefficients, valid for an airfoil of infinite as- 
pect ratio, are used: 


(2rC)o = m = 20, (dF /dk)y = —7/2 


Also shown in Fig. 2 are values of the derivatives as 
calculated exactly from Eq. (4) and as computed from 

“The need for stipulating that the real part of k must be 
Positive if Eq. (7) is to hold arises from the fact that the slope of 
the function F(k) is discontinuous at the origin. If the deriva- 
tive dF/dk is formed by approaching the origin first from the 
left and then from the right, values of +72/2 and —7/2, respec- 
tively, are obtained. 
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Fic. 1. The function G plotted vs. the nondimensional fre- 
quency & for the low-frequency range. Also shown are the 
straight-line curve-fits of Eqs. (8) and (9). 
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Fic. 2. The lift coefficient derivatives cx5 and ¢rp5 as a func- 
tion of the nondimensional frequency k, calculated by three dif- 
ferent methods. The coefficients refer to a control surface with- 
out aerodynamic balance and with hinge line at mid-chord. 


Eq. (1). It is seen that the approximation suggested 
during the present analysis fits the exact results fairly 
well in the range where Eq. (8) is applicable. 

The derivatives currently in conventional use are, 
however, not entirely correct. The cz; from Eq. (1) 
fails to reflect a 10 per cent reduction as k goes from 0 
to 0.05. Also, it is seen that the c,,; from Eq. (1) 
shows severe discrepancies, both with regard to sign 
and magnitude. This same inaccuracy will be reflected 
in the value of the stability derivative ¢,, ... 

Before concluding this section of the paper, it should 
be remarked that it is possible to approach the stability 
problem on a completely rigorous basis, making use of 
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the exact expressions for terms such as /;, given in this 
case by Eq. (4). The stability roots must then be ex- 
tracted by iteration of the stability determinant, as 
was done by Goland and Luke’ and Goland.* * It is 
probable that the added accuracy achieved by this 
procedure, which is of a tedious nature, will be war- 
ranted only when studying aircraft of marginal stability 
where great precision of the results is required. 


THE REMAINING AIRFOIL SECTION DERIVATIVES 


Proceeding in the same fashion as in the case of the 
Cys and c,,5 derivatives, it is a straightforward matter 
to derive the approximate forms for the various re- 
maining airfoil section derivatives, accurate to the 
order k. In the development that follows, in addition 
to the use of the dimensionless time 7 for reference, 
all physical lengths are referred to the semichord b. This 
causes physical velocities to be referred to the forward 
velocity v. 

In addition, physical forces are referred to the product 
pu*b?; hence, physical moments, whether a total airfoil 
moment or a hinge-line moment, are referred to the prod- 
uct pv*b*, This reference procedure for forces and 
moments is somewhat at variance with usual stability 
practice but is particularly convenient for the present 
presentation. 

The coordinates used to describe the motion of the 
airfoil are the same as those of Theodorsen but with 
slightly changed notation. Let y be the absolute, 
downward velocity of the airfoil at a chordwise point 
located a distance a behind the mid-chord. The ab- 
solute airfoil rotation about the point a, positive in the 
nose-up sense, is written as 6. As before, the control 
surface deflection, positive for a depressed trailing 
edge, is 6. The instantaneous angle of attack of the 
airfoil is, then, 


a=y7+0 (13) 
and it is to be noted that a is identical with the instan- 
taneous downward velocity of the craft relative to a 
frame of axes fixed on the airplane. 

As in the last section of the paper, a control surface 
without aerodynamic balance is presumed, so that the 
treatment of reference 3 applies. However, the analy- 
sis can readily be extended to cover cases with aero- 
dynamic balance by utilizing the results in reference 5. 

The remaining lift coefficient derivatives for the air- 
foil section are then found to be 


Cr, = rg = ™ + 2r(dF, dk)o k (14) 


* The author takes this opportunity to point out an error that 


appeared in the manuscript of reference 8. Eq. (11) should read 


Wn = + 16 = Nz + 5) 
and the succeeding development must be corrected accordingly. 
However, the calculations reported in reference 8 were carried 
out by means of the correct analysis and are, to the best of the 
author’s knowledge, correct. 
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Chpy = © + 2m (dG’/dk) se 
Cry = E+ m[(1/2) — a] + 22(dG’ dk) (1p 


The section derivatives for the airfoil moment ¢ 


efficient, referred to the point a and positive for 4 nos. 
are obtained as 


7 l >, (dF 
ak ead Pa ae 


= aw + [a + (1/2)]2r(dG’ /dk) (1g 


l l 
-(5 — a) w+ ( — a’) m + 


up moment, 
( m 7 = 4 ma 


: "py 


mye se 
( ') dG’ 
a+ 2 (19 
2 dk 
, . oo 1\ 71 
ms = — (1; + T 10) + (« + *) 1 x 
é v 
E ‘s (“r) 
ay (9 
aj,\ * 
Om pb ” -[t- == ar+ 5 ta] + 


l Ti T16 dG’ 
a+ 5)|™ a + 27 ~ dé 


The parameters c and the various 7 are defined in refer 
ence 3. 

The section derivatives for the hinge-moment « 
efficient, positive for a moment that tends to depres 
the trailing edge, are found to be 


T2 dF 7 
Cry = Oy, = — E + 2r e ) | (2? 





T, — (Ty/2mr) 2x(dG' /dk) (23 


Choy 


. a l 
Ch = 27. + 7, - (a = ) r| = 
pe 2 


seep 


Tt fi 1G’ 
#1 (; - a) m + ant | (24 
_— TT w ; (<) | 
Cre = _ See — 74Tw) on? |» + 27 dk h 
(25 
is Vu _ Tu (Fu m = in on) (26) f 
Qa 2a \ 2x 7 dk 


Figs. 3, 4, and 5 show plots of the exact values oi 
various of the section derivatives, as well as the value 
obtained from Eqs. (14) through (26). These plots art) 
for pitching of the section about its mid-chord axis and 
for a control surface hinged at the mid-chord. The 
value of dG’/dk used is that appropriate for 0 < k <{ 
0.05. The agreement between the exact and approx! 
mate values over this range is seen to be acceptable. 

Also shown in the plots are the derivative values 
according to the treatment of reference 2. These att 


seen to show serious departures from the exactly calcu) 
| 
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Fic. 3. Total airfoil moment coefficient derivatives Cmy and 


pp a8 a function of the nondimensional frequency k, calculated 
Cmpy 28 é 


by three different methods. The coefficients refer to the moment 
tbout the mid-chord, induced by rotation about the mid-chord. 


lated derivatives. The discrepancies in the hinge- 


moment derivatives are particularly severe. 


THE DERIVATIVES FOR HIGHER VALUES OF k 


Considering complex k = ¢e”, Fig. 6 shows a plot of 
C(k) vs. ¢ (the absolute magnitude of k) for various 
values of the argument v. The data for this plot are 
taken from a paper by Y. L. Luke and M. A. Dengler, 
both of the Midwest Research Institute, in which the 
behavior of C(k) for complex k is comprehensively 
investigated; publication of this paper in the near 
future is contemplated. The important point to be 
gathered from an inspection of Fig. 6 is that, while a 
power series for C(k) written about k = 0 may suffice 
for 0 < ¢ < 0.1 (approximately), it cannot possibly 
afford acceptable accuracy when ¢ falls outside of this 
region. 

While the region 0 < ¢ < 0.1 is generally sufficient to 
embrace all cases of stick-fixed stability and many 
cases of stick-free stability, modes of the latter type are 
sometimes encountered which involve values of ¢ con- 
siderably above 0.1. The forms developed in this paper 
are evidently of no use for an analysis of such cases. 
The derivation of acceptable procedures for use in such 
instances entails great complexity and is presently 
under study. 

One obvious approach is to expand C(k) as a Taylor 
series about a point k different from zero; it is found 
that an expansion about k = 0.4 affords fairly good 
accuracy over a wide region of ¢ above 0.2 and for prac- 
tical values of vy. However, then the expansion is not 
about the point corresponding to the pseudostatic case, 
and this entails considerable practical disadvantage. 
Little can be said at this time about the problem of 
dealing with stability modes involving large ¢, except 
to point out the limitations of the present analysis. 

THE DeRIVATIVES FOR AIRFOILS OF FINITE SPAN 

Up to the present, the development has been based 
on the Theodorsen analysis,*? which presumes an air- 
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Fic. 4. Hinge-moment derivatives c,, and Cap, as a function 
of the nondimensional frequency k, calculated by three different 
methods. The coefficients refer to a control surface without aero- 
dynamic balance and hinged at the mid-chord; the airfoil rota 
tion is about the mid-chord. 
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Fic. 5. Hinge-moment derivatives crs and cays as a function of 
the nondimensional frequency k, calculated by three different 
methods. The coefficients refer to a control surface without 
aerodynamic balance and hinged at the mid-chord. 





Plot of C(k) vs. k for complex values of k = fe». 


Fic. 6. 
foil of uniform chord and infinite aspect ratio. The 
question now arises as to what modification of these 
forms is required in order to make them applicable to 
airfoils of finite aspect ratio and with chord dimensions 
that vary along the span. 

Consider first the spanwise distribution of the param- 
eter m, the section lift coefficient vs. angle-of-attack 
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slope when k = 0. The problem of determining the 
spanwise variation of m is the same as that of comput- 
ing the steady-state loading on airfoils of finite span. 
This problem is of a well-known character, and a variety 
of techniques, thoroughly described in the literature, 
are available for its solution. 

With regard to the calculation of the spanwise dis- 
tribution of the F + jG function, however, a different 
situation exists. The problem is discussed by Reissner,’ 
Reissner and Stevens," and Wasserman,!'! and a further 
bibliography is contained in those papers. A study of 
these references shows that a precise treatment of the 
problem entails great complexity. Not only must the 
spanwise distribution of the F + jG function be de- 
termined anew for each airfoil plan form and. each 
frequency of oscillation, but the nature of this function 
also apparently depends on the, character of the airfoil 
motion—1.e., whether it is of the y, 6, or 6 variety.* 

Suppose, however, that a new function f + jg is de- 
fined so that 

f+nz= 
where, it is recalled, both F + jG and m refer to sec- 
tion coefficients of the airfoil. When k = 0, since 
(2rF)) = mand G = 0, it is evident that the function 
f + jg has the value unity across the entire span, re- 
gardless of the airfoil plan form; however, when k dif- 
fers from zero, this function has values that, in general, 


(24F/m) + j(24rG/m) (27) 


vary across thé span. 

It will how be assumed, for the sake of conciseness 
in the present treatment, that for each value of k and 
for each airfoil plan form an average value of the func- 
tion f+ jg can be deduced which can be taken to act 
across the entire span. Formal methods for arriving 
at this average f + jg can be prescribed on the basis of 
treatments such as references 9 through 11, but these 
are of a relatively tedious nature. 

Since for stability work the function f + jg need be 
defined only for low values of k, it is highly probable 
that a convenient and concise procedure to permit its 
determination for a particular airfoil can be deduced. 
However, the answer to whether or not this is so must 
await the results of a further theoretical study. Actu- 
ally, the parameters of interest to the present treat- 
ment are not the f + jg directly but are the slopes 
(df/dk)o and dg’/dk, where g’ has a definition analogous 
to the G’ of Eqs. (8) and (9). It will be possible, in 
any event, to tabulate values of these slopes for airfoils 
with various standard plan forms, and interpolation can 
be used in the cases of airfoils with plan forms differing 
slightly from these. 

In the absence of accurate values for the slopes, it is 
suggested that those appropriate to an airfoil of in- 
finite aspect ratio be used—i.e., for the average func- 
tion, 

* The present discussion does not touch on the effects of com- 
pressibility, which add further severe complications to an under- 
standing of the unsteady aerodynamic mechanisms. 
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dk)yp = feo = —2/2 


(dfar 
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tf 
UL a0. hese a ae G29 nm 
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0< k< 0.05} 


= —2.086; 0<k< 0.10) 


The chief error arising from the use of these values ;,| 
airfoils of finite span is in the g,’ parameter, wh 
should. have’ smaller absolute values than those 
Eq. (29). 

Finite span corrections must also be applied to th 
apparent: mass terms in Eqs. (11) through (26): they! 
are the terms that. have neither m nor F + jG slope| 
attached. Adequate methods for determining the 
te is . , : . 
form suitable for design-office calculations, and the ys 
of the two-dimensional values appearing in the equ 
tions is therefore recommended. 

The next question to be settled is that of determi; 
ing the lift and moment coefficient derivatives for th 
complete airfoil. In order to account for airfoils wit 
variable chord, it is necessary that a reference lengt 
for the entire system be established; this is general 
taken as the mean aerodynamic semichord of the wing 
However, any other reference length that afford! 
greater convenience for a particular study can be chose 
with equal facility. In the present treatment, th 
reference length is notated 3,. 

All physical dimensions of the airfoil are then x 
all physical times, to v/b,; all physic 
forces, to pv*b,*;— and all physical moments, to po’) 
The nondimensional frequency parameter that describ: 


ferred to 5b,; 


the motion of the system is, then, 


k, = wh,/v 





In the }, system, let y be a spanwise coordinate wit 
its origin at the wing root; let 6’ be the semichord a! 
station y, so that b’ = b/b,; and let 2/ be the win 
span. The local frequency parameter k at station 


is then equal to 


k = b’k, (3 


Consider now the wing lift due to a displacement 
the coordinate y—i.e., a vertical translation of the en} 
The wing lift per unit span at station rf 


tire wing. 
Recalling thei 


first obtained from Eqs. (14) and (15). 
the section derivatives in Eqs. (11) through (26) ar 
nondimensionalized in terms of the /ocal semichord | 
the wing lift per unit span in physical coordinates 1s 


puv"b [¢1,¥ a CLy, (dy dr)| 


and, hence, in the 6, system the wing lift per unit span! 


DD 
o- 


b’ [er, ¥ + &’crp, (dy/dr,)) 


The total wing lift is, then, 
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X where a dot now denotes a derivative with respect to 


time in the b, system 
The lift coefficient derivatives for the entire wing are 
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td | 
moni Ku, “ / b’c, ,dy (34) 
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Kipy = J {0)*tmrdy (35) 
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Consider next the airfoil lift derivatives with respect 


> these : ‘ : ‘ m 
7 2 | to rigid body pitch of the entire surface (6) and with 
; ps respect to symmetrical operation of the control sur- 
é 6 faces (6). Both of these derivatives are of interest in 
St in sgugnhe a ; oy 
th longitudinal stability calculations. It foliows that these 

e us 6 7 . ‘ " . 

‘ derivatives are given by expressions identical with 
e equa ~- - . . 

} Eqs. (34) and (35), save that the subscript y is replaced 

by the subscripts 6 and 6, respectively. 

termir ; . ae si 
for t If use is now made of Eqs. (11) through (16), it can be 
OF the 
4. -.f shown that 
Is wit shown tl 
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t, thie and analogous results follow for the 6 and 6 airfoil deriv- 

atives. 

lk In the case of the airfoil moment coefficient deriva- 
1YSice 


tives, proceeding in the same fashion as for the lift 
pn", derivatives, it is readily shown that 
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tion }} in these equations, expressions for the remaining longi- 

| tudinal airfoil moment derivatives result. The sec- 
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"| by Eqs. (17) through (21). It then follows, for example, 
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Finally, recalling that hinge-moments are nondimen- 
sionalized by pv’b,*, it is seen that 
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Ky, = foray (42) 
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with analogous results for the remaining hinge-moment 
derivatives. The pertinent section derivatives are 
given by Eqs. (22) through (26). 

While attention has thus far been concentrated on 
the airfoil derivatives that enter into longitudinal sta- 
bility calculations, it is evident that the lateral stability 
derivatives can be compounded with equal facility. 
The procedures used for their derivation are directly 
analogous to those employed in the longitudinal cases. 
Thus, for example, the derivative for the airfoil rolling 
moment due to rate of roll ¢, with all units in the 6, 
system, is found to be 


bs 
Kp; = -2 f yer, b’ dy (44) 
0 
Ke; = -2 f ¥eu,,(b')* dy (45) 
0 


The spanwise distribution of m and the slopes fio and 


and also 


g,’ used for the calculation of the section derivatives 


Sk 


must, of course, be appropriate for the rolling motion. 
In the stability equations for an aircraft, the forces 
and moments on both the wing and tail surfaces ap- 
pear; these must be nondimensionalized according to 
a scheme consistent for the entire aircraft system. 
If the semichord b’ and the spanwise coordinate y for 
both wing and tail are referenced to the same 0, (usu- 
ally one-half the wing mean aerodynamic chord), all 
such considerations are automatically taken care of. 


CONCLUSIONS 


It is a natural conclusion of the present study that 
all stability derivatives should be based on the treat- 
ment outlined here, rather than on the procedures re- 
ported in the past literature. At the same time, it is 
important to note that the use of the more correct 
derivatives proposed in this paper does not increase the 
degree of complexity required to carry out a stability 
solution over current levels. 

Wherever possible, of course, stability derivatives 
should be evaluated by accurate prototype or model 
testing; in the usual instance, only the pseudostatic 
(k = 0) values of the derivatives can be conveniently 
determined in this manner. These pseudostatic values 
must then be extrapolated to the case of unsteady 
motion by the methods described in this paper. For 
example, a static wind-tunnel test would provide the 
value of A,, at k = 0; according to Eqs. (42) and (22 
this value is 


; "Tw. 
(Kn,)yo = — [ =m (b’)? dy (46) 
l 


9 


. 7 


and the value of the derivative when k, is not zero must 
then be taken as 
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1a 


kon 
Kn, - (Kn,)o — ky fro / == m (b’)* dy (47) 
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As a final remark, it is of interest to note that certain 
of the findings in this paper were anticipated in prior 
studies conducted by Gruschwitz,'* Sdhngen,'* and 
Smilg.'* These authors point out that the stability 
derivatives in conventional use do not adequately de- 
scribe the aerodynamic damping present in the system. 
They also point out that the negative damping pre- 
dicted by the Theodorsen derivatives for certain os- 
cillation conditions will probably account for the 
lateral ‘‘snaking’’ phenomenon encountered in some 
aircraft designs. The work of Gruschwitz is of fur- 
ther interest in that he brings out the unusual behavior 
of the damping derivatives called for by the Theodorsen 
theory as the pseudostatic condition is approached; it 
is shown here that this is associated with the infinite 
value of the slope of the G vs. k curve as the frequency 
k approaches zero. 
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APPENDIX | 


The Steady-State Air Forces Induced by Oscillatory 
Nonharmonic Airfoil Motions | 
The airfoil motions to be considered here are thy, 

that proceed according to the schedule e’’”, where 


A=k — jé (A 


and where, by definition, both & and & are real qua 
tities. The motions are seen to be oscillatory with, 
nondimensional frequency k, but they are nonharmoy) 
when ¢ is different from zero. A negative value for! 
defines a stable motion (amplitude decreasing wit 
time), while a positive & implies unstable oscillation 
The natural oscillatory motions of an aircraft m 
be stable, neutrally stable, or divergent; hence, a st; 
bility analysis of the system must employ a statement 
the air forces which is capable of embracing all thre 
such types of motion. 

In the body of this paper, a study is made of th 
steady-state air forces induced by sinusoidal motions, 
the airfoil. The term steady state is used to denote tk 
fact that only those components of the air forces a 
considered which, like the motion of the airfoil, pr 
ceed according to the schedule e”’. Physically, thi 
has the implication that the air forces are dealt wit 
after the wake is infinitely long and after the effects: 
aerodynamic starting transients are no longer present. 

These steady-state air forces are, of course, shownt 
be a function of the frequency k. It will now be show 
that, by the simple replacement of k by the paramete 
\ of Eq. (A1), the harmonic air-force statements becont 
applicable to cases of nonharmonic oscillatory motion 

Consider motions of the airfoil which are defined bj 
a particular time history of the displacement coordinati 
x(r); this coordinate may be a translational velocity] 
the angle of attack, the control surface displacement, 
the like. Furthermore, let «(7) represent the circulator 
air forces that arise from a unit step-function displace 
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ABSTRACT 


Based on the assumption that the aerodynamic forces acting 
at any point on a lifting surface are linearly dependent upon the 
local angle of attack at that point, the theory of pure rolling mo- 
tion is applied to the experimentally observed motion of a finned 
missile in supersonic flight in order to (1) determine how well the 
actual motion is represented by the theory; (2) determine the 
suitability of the Aerodynamic Range at the Ballistic Research 
Laboratory, Aberdeen, Md., to the proposed free-flight roil 
technique; and (3) determine certain aerodynamic coefficients 
associated with the motion over a range of Mach Numbers from 
1.67 to 2.7. Values for the aerodynamic coefficients are also de- 
rived using the linearized supersonic theory and are compared 
with those obtained from experiment. 


The results indicate that the actual motion is well represented 
by the theory to within the small experimental errors and that 
excellent reproducibility of the aerodynamic coefficients in roll 
is obtained. Furthermore, the results show a fairly good corre- 
lation with the linearized theory when applied to aerodynamic 
surfaces of 16 per cent thickness as employed on the missiles used 
in the reported tests. 


INTRODUCTION 


[" THE PAST, the activities of the exterior ballistician 





= those of the aerodynamicist have been widely 

separated. The ballistician was primarily interested in 
rpwingless and finless spinning bodies of revolution travel- 
ing at supersonic speeds, and the aerodynamicist was 
interested in subsonic aircraft whose primary elements 
However, in recent years, a 
change has taken place—the ballistician is now 
equipping his bodies of revolution with wings and/or 


at, o fins and calling them missiles, and the aerodynamicist 


ator 
lace 


has vigorously extended his field of interest and ac- 
ivity into the supersonic régime, with the consequent 


his (@@PPlication to supersonic aircraft and guided missiles. 


irc} 
x(rf 


he two fields have thus merged and the ballistic and 
erodynamics people can mutually contribute and 
benefit from their respective work and tools. 


Presented at the Aerodynamics Session, Eighteenth Annual 


i Meeting, I.A.S., New York, January 23-26, 1950. 


* Part of the material presented in this paper appeared in the 
experimental section of a Doctoral Dissertation for Yale Uni- 
Wersity (June, 1949) by Ray E. Bolz. The authors wish to 
acknowledge the contributions of Dr. A. C. Charters to the gen- 
tral program and the assistance of Miss Karolyn Skegas in the 


’ Preduction of the data. 


t Associate Professor of Aeronautical Engineering. 
+ ° - ~ Ts . 
t Aerodynamicist, Free-Flight Aerodynamics Branch. 
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A Method of Determining Some 


ents from Supersonic 


Free-Flight Tests of a Rolling Missile’ 


HN D. NICOLAIDES? 


Case Institute of Technology and Ballistic Research Laboratories 


One of the experimental tools developed by the ex- 
terior ballistician appeared particularly applicable 
to the general field of experimental aerodynamics 
namely, the Free-Flight Range at the Ballistic Research 
Laboratory, Aberdeen, Md., under the direction of 
Dr. A. C. Charters.'? In order to explore a particular 
phase of this applicability, the authors, with the willing 
consent and cooperation of Charters, several years ago 
embarked upon an experimental program in the Aero- 
dynamic Free-Flight Range to evaluate the usefulness 
of the equipment for accurately yielding not only static 


but more important dynamic force and moment 
coefficients associated with given winged and/or 
finned aerodynamic configurations in supersonic 
flight. 


The general experimental technique is: The motion 
of the missile is determined by measurements taken 
from successive photographs of a configuration during 
free-flight down the range. The integrated equations 
of motion are then “fitted” to the measurements by a 
least squares procedure, and the constants of the final 
equation yield the experimental values of the aero- 
dynamic coefficients. The accuracy of the procedure 
then depends upon the validity of the equations of 
motion, and a correct theory must necessarily result in 
a final equation that fits the measurements within the 
experimental error. 


The first project in the overall program is presented 
in this report, wherein a simple motion—namely, that 
of pure roll—is imposed on the missile and the theory 
of pure rolling motion is applied to the experimentally 
observed motion in order to obtain the associated 
aerodynamic coefficients. The fit of the equation 
of rolling motion to the range data is obtained, and a 
comparison is made of the experimentally obtained 
coefficients, one static Ci, and one dynamic C,,, to those 


predicted by a linearized theory. 


It is not intended that this paper present the results 
of a finished research program. However, it does in- 
clude the basic theory of motion, the specific range- 
test technique adopted, the method of data reduction, 
some specific results over a series of Mach Numbers, 
the experimental accuracy and reproducibility of the 
tests, and the overall evaluation of the range as applied 
to the technique devised. 
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Superscripts 


additional test data at the higher Mach 


SYMBOLS 


radius of body of missile 
area of fin or wing panel (see Fig. 1) 


arbitrary constants depending on boundary 
conditions 

wing span 

VM? - 1 

wing chord 

damping constant |(A7,/7) — (Kr/m) 


Ki;/1 

drag coefficient 

roll moment derivative due to canted surface, 
(OC)/06), 6in rad. 

roll moment derivative due to rolling velocity, 
dC1/d(pb/2V) 

roll moment coefficient, L./gA p> 

diameter of body 

slope of the lift curve at any point x, y on fin 

L3/V? 

L,/pV 

(p/2) SCp 

rolling moment due to angle of cant 

rolling moment due to rolling velocity, damping 
moment 

axial moment of inertia (slugs-ft.*) 

Mach Number 

mass (slugs) - 

number of canted fins 

total number of fins 

rolling velocity (¢), rad. per sec. 

steady-state rolling velocity (¢).., rad. per sec. 

probable error 

fin-tip helix angle 

steady-state fin-tip helix angle 

(3 /4)d? ¢ 

steady-state rolling velocity (deg. per ft.) 

linear velocity of missile along trajectory (ft. 
per sec.) 

Cartesian coordinates: y along fin span; x along 
fin chord 

linear distance along the free-flight range 

angle of attack 

angle of cant on one fin 

tip Mach cone angle 

polar coordinates shown in Fig. 3 

free-stream air density, slugs per cu.ft. 

angle of roll 

i= 


Dot over symbol is time derivative of variable 


Prime mark is the derivative with respect to 


range 


s, distance along 


GENERAL EXPERIMENTAL TECHNIQUE 


A suitable aerodynamic test vehicle is gun-launched 


down the Aerodynamic Free-Flight Range. 


The spe- 


cial firing technique employed results in small pitching 
and yawing displacements of the test vehicle during the 


subsequent flight. 


This vehicle may be a model of any 
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1950 


aerodynamically stable missile configuration having 
one or more pairs of fins set at a differential angle 
incidence (cant) to induce a rolling velocity, Duriod 
flight, as it passes each of 25 stations accurately logy, 
along a 218-ft. distance, the missile is Photograph, 


(shadowgraphs) simultaneously in the vertical 
horizontal planes by means of an extremely gh, 


iT 


The time of occurrenc 


meme 


duration spark discharge. 
the discharge for nine of the stations is measured ;| 
within 10-6 From the photographs, the Spat 
orientation of the missile is obtained, from which ; 


sec. 


linear displacements (to 0.001-ft. accuracy) and angy 
displacements in roll (to +0.7° accuracy) are dete 
mined at the 25 discrete distances along the trajectoy 
as well as at nine discrete values of time. 

With these data obtained from the aerodynay 
range and the analyses in the following sections, 
aerodynamic coefficients associated with the rolj 
motion may be determined and compared with ¢ 
theoretical values calculated on the basis of linear 


supersonic theory. 


— 


THEORY 


(A) Equation of Motion of a Missile in Pure Roll 


If a symmetrical winged and/or finned missile wif 
canted surfaces is in free flight and has an angu 
motion of pure roll (yawing and pitching displaceme 
negligible), the missile may be considered as being act 
upon by two roll moments. One moment, tending 
increase the roll velocity, results from the lift-o 
couple produced by the canted surfaces in linear forw 
motion; the second moment, tending to decrease | 
roll velocity, is the damping moment produced by 
lift-force couples induced on all the aerodynamic si 





faces by the roll velocity itself. 

(1) Roll Moment Due to Canted Fins.—The 1 
moment induced by n, differentially offset (cant 
wing or tail panels may be expressed as the integral 
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Fic. 1. Schematic diagram of canted tail surfaces. 
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1 having product of the aerodynamic lift force and the moment 
angle : arm of this force about the missile axis. Consequently, 
= from Fig. i, 
Y Toca nzpV° "b/2 (dC, , 
ograph| lL, = — / J . & ) yidxdy (1) 
1cal rs . came P 
LY shor L, = Ki, V? (2) 
rence ‘ ; ; ° 
sured (2) Roll Moment Due to Rolling Velocity, Damping 
€ spat Moment. When a missile rolls about its own axis, 
hich i) the lifting surfaces are subjected to a normal velocity 
langyg vector induced by the rolling velocity. The induced 
re dee} velocity results in an induced angle of attack on the 
ajectonp surface which varies over the surface as a linear func- 
tion of the radial distance from the rotational axis. If 
dynapt the roll damping moment, neglecting the viscous damp- 
ions, ty ing on the fuselage,* is defined as the roll moment in- 
» roll, duced by and opposing or damping the rolling velocity, 
vith i then this aerodynamic damping moment (with refer- 
nearing ence to Figs. | and 2) may be expressed by a double 





| integration as follows: 


* *h 2 fie ) 
Pn L ‘ 
_—T - y? dx dy 3 
Ly P 2 / a) a = zy aot sa 









or 

sile wf L, = K lp P V (4) 
— (3) Dynamical Equation of Roll.—-Applied to a 
ceme® symmetrical missile in pure rolling flight, subject to 
18 iH small pitching and yawing displacements, the dynami- 
ding ' cal equation in roll becomes, simply, 

tf t-Lor 

forw Ig =L,—L, (9) 
“ast "§ Here, the rolling moments exerted by the wing and tail 
| by! surfaces arising from the pitching and yawing angular 
nic S displacements and angular velocities have been neg- 
| lected because (1) for small pitching and yawing dis- 
he 1 placements (mean square yaw less than 1.5°) these 
cant’® moments are extremely small and of second order com- 
€gru® pared with L; and L, (see Appendix I) and (2) they are 


oscillatory with respect to their direction of action and, 
as such, any small effect that they may otherwise exert 
on the resulting rolling motion tends to cancel out. 


(B) Solution of Dynamical Equation of Roll 


For the missiles traveling the length of the range 
between stations (218 ft.), the decrease in velocity is 

/ | less than 7 per cent, and, for values of flight Mach Num- 
ber in excess of 1.6, the variation of dC,/da, the slope 

) of the lift curve, resulting from this small deceleration 

is less than about 3'/. per cent from the average value. 
\ || Therefore, little error is introduced if the values of 


ae ; ? — 
lhe effect of viscous damping of a missile fuselage may be 
estimated from the performance of spinning bodies of revolution 


as determined in the aerodynamic range. The numerical value 


NES of the viscous damping coefficient of our missile is approximately 
0.16 per cent of the aerodynamic damping coefficient obtained 
Consequently, the effect of body damping may then be neglected 

3. for the purposes of this report 
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Fic. 2. Cross-sectional view of a canted tail surface at a given 


span location showing the effect of rolling velocity. 


Ky, and Ki, in Eqs. (2) and (4) are considered as con- 
stant for any given flight. 
Insertion of Eqs. (2) and (4) into Eq. (5) gives 


I¢ + Ki ¢oV = Ki, V? (6) 
where 
. n c *b/2 (dC 
Ki, = ss ( ‘) y6 dx dy (7) 
= owe sue GOT as ; 
. n (¢ (/2 (dC, 
Ky = ws ( :) y? dx dy (8) 
2 0 y=a da ry : . 
and @ = d@/dt = p, the rolling velocity. 


The approximate linear equation of motion of the 
missile along the trajectory may be written 


mV = —K,V? (9a) 


where 


Ke => (p 2) SCp (9b) 


is considered constant for the small variation in ve- 
locity V. 

Eqs. (6) and (9) are now transformed to the inde- 
pendent variable, z, the horizontal distance along the 
range, by 


d dzd ai 
dt dtdz ~~ dz | 
d? a _dv d ” 
pm” ee 


This transformation as written is valid, since in the 
range the distance along the trajectory may be as- 
sumed equal to the distance along the range and the 
missile velocity becomes V = dz/di, with negligible 
error. If a prime is used to denote a derivative with 
respect to z, then Eqs. (6) and (9a) become 


"+ 9’ (— + =") ee 1 
“ r* 2 lI ie 
and 
V’/V = —Kp/m (12) 
Substitution of Eq. (12) into Eq. (11) gives 

Ky Kr Ki, 
¢” + o'( ? | ——'=0 (13) 

I m I 


a linear differential equation with constant coefficients. 
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Fic. 3. Schematic view of tail surface illustrating integration 
procedure 


Defining 
Ci = (K,,/1) — (Kr/m) (13a) 


Cy = K I (13b) 


'5 
Eq. (13) becomes 
¢” + Cd’ — G = 0 (14) 


The general solution of Eq. (14) may be written im- 
mediately as 


@ = B+ sz + Ae ‘" (15) 
where 


5 = C C; (16) 


steady-state rolling velocity. 

A and B are arbitrary constants that depend upon 
the boundary conditions. A curve given by Eq. (15) 
relating the missile roll angle and the distance along 
the trajectory is fitted to the data points by means of 
the method of differential corrections. The process of 
‘fitting’ then results in a determination of the roll con- 
stants C; and C; (and therefore the constants K,, and 
Ki,). 

Thus the damping and rolling moments can be 
experimentally determined as a function of missile de- 
sign and Mach Number. 


(C) Theoretical Determination of the Roll-Moment 
Coefficients 


The values for A,, and Ki, defined in Eqs. (7) and 
(8) are determined experimentally in the aerodynamic 
range by the technique discussed in this paper. How- 
ever, approximate values for these two coefficients may 
also be determined theoretically from the application of 


the existing linearized supersonic theory. 

The theory as applied here assumes the two-dimen- 
sional value of dC,/da in the areas of the fins outside 
the tip region for the calculation of K,, and K,. This 


assumption, of course, cannot give the exact theoretical 


lift distribution on the wing within the Mach qy, 
originating from the intersection of the wing leadin,§ 
edge with the body. However, since this regio: 
(where the lift distribution is only approximately cl a fun 
culated because of the neglect of the body-wing jun K, fo 
ture interference effect) is close to the body, the », 
sultant error in the theoretical rolling moments shoy 
be small. Fre 





The effects of body upwash on wing lift when th, 
body is at an angle of yaw have also been neglecte| 
However, in the roll tests in the aerodynamic range, i 
should be noted that yaw during flight down the ranp 
is always intentionally restricted by the launching teq 


where 


nique to an extremely small angle (mean square yay 
less then 1.5°) in order for the roll Eq. (6) to be near 
exact. Any upwash on the wing resulting from th 
yawed body would result in (1) small symmetric; 
forces (lift forces in the same direction on both halve 
of the fin) when the fins are oriented at 0° or 45° 
the plane of yaw, and (2) possibly very small unsymme 


An 


Frc 
give 
rical forces due only to flow dissimilarities when th 
missile is oriented in roll to any angle, other than thos 
stated above, referred again to the plane of yaw. Sino . 
a roll moment can obviously arise only from the wu} 4) 
symmetrical forces, it may be concluded from th gratic 
qualitative observations above that any roll momen} oytgi 
resulting from body upwash phenomena would not on! 
be small moments but would be rapidly oscillating i 
sign as well as in magnitude. Such moments woul iii 
contribute negligible effect to the rolling motion of th inthe 


missile resulting from the relatively large and contin Re 


ously directed moments, L, and L,. 


One further point may be noted—namely, that th 





forces on the cone-cylinder fuselage cannot contribu! Fr 
to a rolling moment other than the negligible viscou 
damping moment already mentioned. 


(1) Determination of K,, from the Linearized Supe 
sonic Theory.—From the linearized two-dimension and 
theory of Ackeret,’ the value of the lift slope at @ 
points on the surface of an infinite-aspect-ratio liftin 
surface is given by the expression 


() = 4 = 4 a Fr 
da/., WM?—-1 B (K,) 
which is a function of Mach Number only. Howevel to Ec 
for a finite wing or tail surface, the pressure distriby} 


tion over the lifting surface included within the tip Mac} 
cone is influenced by the wing-tip boundary condition 


as well as by the variable angle of attack (py/ V hel (3) 
angle) induced over the surface by the rolling velo Coeff 
ity. expel 

This pressure distribution is evaluated on the basis in te 
the three-dimensional linearized supersonic theory coefh 
reference 5. The value, then, for (dC,/da),., , for at By 


incremental area on a rectangular wing within the t 
Mach cone boundary is 
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(S) 8 —_ B /b ) | GC )( b + x (1 
= sin — —y)}— ~ ,— 
da/yzy B ; x G i y 2 7)? 2 B 


a function of position (x, y) on the wing, as well as of the flight Mach Number. From Fig. 3, then, the value of 
K,, for the missile is 
Ki, = (Ku) area I — (K wy) area nT (A wens WI (19) 


From Eqs. (8) and (18), using d.S in place of (dx dy) for the differential area, 


: or 8h. B (b l b b x 3 » 
(Kip): = > Fa <B sin “i 33 = 2 gInY— > + ) y* ds (20) 


where, from Fig. 3, 


dS = rdudr = (x/cos p) (dx/cos p) dul 


y = (6/2) — x tan p { (21) 
An integration of Eq. (20) using Eqs. (21) results in 
(Kip), = (4en/B?*) [(c?b?/32) — (7c*b/96B) + (17c*/384B"*) | (22) 
From Eqs. (8) and (17), integration for (K,,);, and (K,,);;; is accomplished, and these values, added to Eq. (22), 
give 
. 2pn c* cb cb? ch® ca 
oa bs + 4sB? 16B) 24. 3 ) (23) 


(2) Determination of K,, from Linearized Supersonic Theory.—The expression for A,, is determined by an inte- 
5 - '5 a 


gration of Eq. (7). Here again, Eq. (17) gives the value of the lift slope on the wing or tail surface at all points 
outside of the area bounded by the Mach cone. The lift slope for the area within the Mach cone is given by 
(dC,/da),, y = (8/7B) sin“! V IBl(b 2) — y]}/x (24) 


derived in references 5 and 6. This value differs from Eq. (18) because the lift forces that give rise to the roll 
moment L; are a result of a uniform angle of attack at all points on the surface—namely, the angle of cant. 
Referring again to Fig. 3, 


K, et (Kj, area! 5 (Ki, )area il 5 (Ki) area il 


é 


From Eqs. (7), (21), and (24), 
; pns6 ; “tan '(1/B) \ 8 B[(b 2) _— y] / yx 
(Ki), = 7 = > ——— i ae 
2 J x=0S p=0 lB x \ COS“ 


(Ki), = (pms6/B*) [(c*b/4) — (ce? 4.8B)| (25) 


and 


From Eqs. (17) and (7), the values of (Kadi and = and 
(Kida are determined, and the addition of these values pb 
ie a Ci ( 


2V 


Joa pon = K,,Vp 


to Eq. (25) gives 


: n36 (b% e cb sii 
(K,,) = pns ( 4 ~ @% ) (2g) Then, 


B 4 SB? 4B : 
; , . P af c a 2A 5 yt f om 
(3) Determination of the Conventional A erodynamiu ls = - 
ra ae 5: ae ‘ ; pbA,n56 pbA ,ns6 
Coefficients in Roll.—It is desirable to express any (27) 
experimentally or theoretically determined coefficients C —4K,, —4/] : Kr 
in terms of the conventionally defined aerodynamic _—. pb?A,n ial pb?A ,n m 


coefficients of roll, C,, and C 
. ’ Ip © 15° 7 . " ‘ . 
By definiti Comparisons of the experimental and _ theoretical 
’ definition ares J _ 
: : values of C,, and C,, offer a means of evaluating the 


Ls; = Ci, gApns6b = Ki,V" accuracy of the linearized supersonic theory as applied 
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Fic. 4 


Typical shadowgraph of missile in free flight showing 
locating pins in the tail fins and supersonic flow phenomena. 


It will 
be noted from Eqs. (27) that accurate measurements of 
the actual angle of incidence of each tail and wing 
panel is required for the theoretical determination of 


here when used to calculate these quantities. 


C,, for any missile tested. 


METHOD OF DATA REDUCTION 


The experimental values of the aerodynamic coeffi- 
cients in roll for the missiles tested are determined by 
the following general procedure for reduction of the 
data obtained from the aerodynamic range. 

(A) Accurate measurements of the angular orien- 
tation in roll (@) and the linear position (zs) of the 
missile in space are obtained from each of the photo- 
graphic plates taken during the flight of the missile 
down the range. Accurate measurements of the time 
at which nine of the photographs were taken are deter- 
mined from the chronograph record and megacycle 
counters. 

(B) The form of the solution @ = ¢(s) [Eq. (15)] of 
the theoretical roll equation is ‘‘fitted’’ to the meas- 
ured values of @ and z, and the optimum values of the 
constants C,; and C, in Eq. (15) are determined by the 
method of differential corrections. 

(C) The value of the retardation constant A ,/m for 
the missile is determined by existing techniques, using 
the measured values of z and the time ¢. 

(D) 


are determined and their statistical accuracy (P.E.) 


Values of the aerodynamic coefficients in roll 


is computed. 
The details of these four general steps will now be 
discussed. 


AERONAUTICAL 


SCIENCES—OCTOBER, 1950 


(A) Measurements from the Photographic Plates 


(1) Angular Orientation in Roll. 
mine the roll angle from the shadowgraphs of the | 
missile (see Fig. 4), it is first necessary to identify ¢ 
quadrant in which a particular reference fin OF Win 


In order to dete ' 


occurs. 
nizable in both the vertical and horizontal phot 


Therefore, one particular fin must be rege. | 
graphs taken at each station. Identification of the fi, 
was obtained, after considerable experimentation, } 
locating a small pin in the trailing edge of each fin ne, 
the tip. The shape of the four pins differed as ge, 
in Fig. 4, and it was always possible to distinguish » 
least one of the shapes in both photographs and there} 
definitely fix the quadrant of roll orientation. 1), 
roll angle of the missile at any station may be calculate 
from the coordinates of each of the identification pigs] 
[Eq. (28)] as obtained from the horizontal and vertic; 
spark photographic plates using the usual methods 
reduction of range data for obtaining the spatial pos 
tion of a point from its images projected on the hor 





zontal and vertical plates. Then, from Fig. 5, 


tan! lh = 

AD 

ar ~ EF 
EH 


AC — CD _ 
AD 
EG — GH _ 

EH 


Psta = tan ' 


tan an 
The angular accuracy of the roll angle at a station i 
estimated to be +0.7° from a comparison of the fou 
individual determinations made according to Eq. (28 

2) Linear Position of Missile-—The location of « 
reference point on each station along the range is know 
to 0.001 ft. 
graphs, and the location of the center of gravity of the 


This reference point appears on the phot 


missile in space at each station is accurately an 
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Sketch illustrating the spark station geometry. 





Cc 





Fic. 5. 


eal 








simp! 
refer 


(B) 
TI 
const 
theo! 
g, 2 
ess ( 
appr 
nate‘ 
the 1 
optit 


(1 
roll ¢ 
of le: 
of ec 
be e 
of th 


o= 


in W 
may 


Om 


whet! 


initi: 


is a 
(AB 
stan 


ess. 


may 
may 
resic 
not 
repr 
initi 
dete 

T 
of di 
dete 

“7 
angle 
the 1 
missi 

+] 
succe 
the p 








bs 
to deter : 
S Of th] 
ntify thy 
OF Wing 


recog. | 
| phot 
the fin: 
tion, }y 
fin neg 
aS sey 
Tush a 
thereby 
1 Th 
culate 

on pins 
verticg 
hods , 

a] posi 
1e hor; 





(28 


tion is 
le four 
. (28 
n of a 
SNOW! 
yhoto 
of the 

anc 











SUPERSONIC FREE-FLIGHT 
simply determined by the technique presented in 
reference 3. 

(B) Determination of the Constants in the Roll Equation 


The constants C, and C2, and at the same time the 
constants A, B, and s, are determined by “‘fitting”’ the 
theoretical form of the roll Eq. (15), @ = ¢(2), to the 
, z data obtained from a particular firing. The proc- 
ess of “fitting” consists first of finding the initial or 
approximate values of the constants in Eq. (15) (desig- 
nated as Ao, Bo, So, and C,) and second of applying 
the method of differential corrections to determine the 
optimum values of the constants. 

(1) Method of Differential Corrections.—Since the 
roll equation contains an exponential term, no method 
of least squares has been developed to handle this type 
of equation directly. However, the roll Eq. (15) may 
be expanded by a Taylor series about the initial values 
of the constants 
74+ AB + 2As + 

ze “AC, AA + 
“#1(AC,)?/2]) +.. 


@ = Bot soe + Avge 
e ‘AA —sAve “FAC, 


2*Ave 


in which all terms in AC; and AA of order 2 or greater 


may be neglected if AC; and AA are small. Then, 
dm —~% = AB+2As+e AA — Ayse~°* AC, 
where 


—Cig2 
Po = (Bo oa S02 + Aoe -? 
initial value of roll angle, and 
om = measured value of roll angle 
is a linear relation, and the differential corrections 
(AB, As, AA, and AC,) to the initial value of the con- 
stants may be obtained by the usual least squares proc- 
The corrected constants 


ess. 


51; = So As 
C = Ci, + AC, 


A, = Ao + AA > 
B, = By + AB; 


(29) 


may in turn be used as the initial values, and the process 
may be repeated until the sum of the squares of the 
residuals* approaches a minimum. If the process is 
not convergent, then either the roll equation does not 
represent the motion with sufficient exactness or the 
initial values for the constants A and C, are not well 


determined. 


Therefore, the successful application of the method 
of differential corrections depends upon a fairly accurate 
determination of the initial values of the constants. 


* Residual here is defined as the difference between the roll 
angle of the missile as measured and the angle calculated from 
the roll equation at any distance, z, along the range where the 
missile was photographed. 

t If the sum of the squares of the residuals decreases with 
Successive applications of the method of differential corrections, 


the process is convergent. 
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The technique employed to obtain them is discussed in 
the following section. 

Determination of the Initial Values for the Roll 
Differentiation of the roll Eq. (15) yields 


(2) 


Constants.t 


¢’ = s — CjAe “ (30) 


and substituting the value of Ae~“* from Eq. (15) into 
Eq. (30) 


o’ = (s + CB) + (CGs)z + (Cio (31) 


The value of ¢’ may be computed at each station 
by the method of divided differences,’ using the known 
values of ¢’s and 2’s in Eq. (15), @ = f(z). 


f'(2)o = (’)o = f(Z0, 21) + (20 — 21)f (20, 21, 22) 


where 0, 1, and 2 are consecutive stations along the 
range and f(Z, 2:) and f(2o, 21, 22) are defined in reference 
7 ** 


The ¢’ at the last two stations, (x — 1 and n), how- 
ever, cannot be computed by the above equation be- 
cause the second difference cannot be obtained. The 
¢’ at (n — 1) station is given by 


2)f (Sn 2) on ly Zn) 


, 
(@ )n [= I(Sn 2) Sn ay + (Zn po Be 
and at the last station by 
, %e e - le ° « 
(p lj, = S (Sn 2) Zn-1) -+- (2s. “— «n—2 ~~ #n i)f (Sn 2) ~n—ls Zn) 


With the ¢, ¢’, and z values known for each station, 
Eq. (31) is solved for the constant C, by the method 
of least squares. Using this value of C,, Eq. (15) is in 
turn solved by the method of least squares and values 
of s, B, and A are determined. These four values 
are then the initial values C),, so, Ao, and Bo as required 
by the method of differential corrections. 

(3) “Fit” of the Roll Equation.—The final equation 
representing the rolling motion of the missile is Eq. 
(15), written with the optimum values of the constants 
derived from the data by successive applications of the 
method of differential corrections until the sum of 
the squages of the residuals approaches a minimum. 
The fit of the final roll equation to the range data may 
then be evaluated by a comparison of the statistical 
magnitude of the final residuals with the magnitude of 
the estimated error in the measurement of the roll angle 
from the photographic plates. The statistical magni- 
tude of the final residuals is given by the usual statis- 
tical quantity Probable Error,’ called here ‘‘the Prob- 
able Error of the residuals’’ 


fxn. = 0.6745 Vs Adq?/(n — 4) (32) 
where 
A@ = residual 
n = number of observations 


t Procedure suggested by Dr. Zdenek Kopal, Professor, Massa- 
chusetts Institute of Technology. 
** Page 21. 
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Aerodynamic Free-Flight Range showing the spark 
stations located along the flight trajectory. 


Fic. 6. 


If the P.E. of the residuals is of the same order of 
magnitude as the estimated error in measurement of 
the roll angles from the shadowgraphs (+0.7° approxi- 
mately), then the final roll equation represents the 
best equation that can be written to fit the data and, 
as such, represents the true rolling motion of the missile 
within experimental accuracy. 


(C) Determination of the Value of the Retardation 
Constant Kr/m 


The retardation constant defined by Eq. (9b) is pro- 
portional to the coefficient of drag of the missile. The 
constant is then determined in the conventional manner 
described in reference 3, where it is assumed that the 
time is given by a power series in the distance meas- 
ured from the center of the range (zo) as 


t = do + ay(2 — 20) + ae(z — %)?> +... 


The coefficients do, a1, d2, etc., are determined from 
the experimental measurements by the method of least 
squares, and the drag coefficient is, in turn, determined 
from these coefficients. 





“Railed”’ gun for launching missiles. 


Fic. 7. 


AERONAUTICAL 


SCIENCES—OCTOBER, 1950 
(D) Aerodynamic Coefficients and the Statistical 
Precision 


The values of the constants C; and C; of the roll] equa 
tion, as determined by successive applications of thy | 
method of differential corrections, may be considere; 
as the best values that the experimental data can yielg 
consistent with the assumptions made in writing the 
equations of motion. The roll coefficients C,, and C, 
are then determined simply by Eq. (27). . 


The effectiveness parameter, ~,b/2V, the steady 
state tip helix angle, may be obtained from the linear 


psd = b (5) i 
2V 2(57.8) \K,/ (33 


ized theory by 


and, from the experimental data, by substitution oj 
Eqs. (13a) and (13b) into Eq. (33); whence, 


pb _ _b =| Co | 
2V 2(57.3) LC: + (Kpr/m) (34a) 
and using Eq. (16) gives 
pb sb ( = . 
2V.—_-.2(57.3) Cm (34b 


The accuracy of the final aerodynamic coefficients f 
is again expressed by the statistical quantity, PE, 
dependent upon the P.E. of the constants in the roll 
equation. The general equations for the P.E. of any 
function of independent quantities whose P.E.’s are | 
known is given by’ 


Fue Yon Ow 
R = 24 () 24...() +2 @ 
Fe net dg. ro? + ng 7,° (35 


where 





Q = F(g, 2> qs; ce ey Qn) 


and R, 1, r2,..., 7, are the P.E.’s of Q, qi, gz, -- +5 Qu | 


respectively. 


EXPERIMENTAL TESTS 


The purpose of the experimental program is not to 
present a finished research program but rather to in- | 
vestigate the suitability of the aerodynamic range to 
the free-flight roll technique presented in this paper. | 
Specifically, the purpose may be restated in three 
parts: 


(1) To determine the accuracy of agreement of the 
theoretical roll equation with the experimental data 
measured in the aerodynamic range, as well as to de- 
termine the accuracy of the reduction methods used | 
to obtain the roll constants (aerodynamic coefficients) 
from the data. 

| 

(2) To determine the ability of the combined range | 
and technique to reproduce accurately aerodynamic 
coefficients. 
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SUPERSONIC FREE-FLIGHT TESTS OF A ROLLING MISSILE 617 
om 
_ 7.850 %D. 
7 " 
AN 
‘ . 
20° 3 | 
- - — me ae 
Fic. 8. Design of test missile 
(3) To compare the experimental and _ theoretical en 
values of the aerodynamic coefficients. 
Average Moment of 
; Mis- Angle of Inertia in 
(A) Experimental Apparatus and Procedure sile Cant on Roll, J Mass of 
; - : Num Two Tail (Slugs Missile, 
(1) Aerodynamic Range.—The aerodynamic range ber Fins, Deg. Material per Sq.Ft.) m (Slugs) 
and the test procedure are described briefly in the first I 2.029 — 1.240 X 10° 1.224 XK 10°? 
. i 2 . ’ uselage 
section of this paper, and further details are presented plus 
in reference 3 and Fig. 6. The one unique modifica- bronze 
: , : a ss fins 
tion to the equipment is the special ‘“‘railed”’ gun (shown 2 4288 Aluminum 1.235 1 29] 
in Fig. 7) used for launching winged and/or finned mis- — 
. . - . : ° US 
siles. In the operation of the gun, the projectile is el 
guided by its body in the central cylindrical portion of fins 
. : “ith: . er a _ ; ‘ 3 2.108 All bronze 2.186 3.068 
the barrel, the fins or wings riding freely in the slots. | yt? te |) en et 3 067 
(2) Design of Missiles for Firing Program.—To ful- 5 1.933 Aluminum — 1.266 1.226 
? _ ? ie : fuselage 
fill the purposes of the program, a single missile design plus 
was chosen consisting of the simplest configuration from bronze 
‘ ; . ‘ fins 
a construction standpoint, as well as from a standpoint 6 2 100 Aluminum 1.263 1.223 
of amenability to theoretical supersonic aerodynamic a 
’ _ ae ‘ - : ee ; lus 
calculations. The missile design is shown in Fig. 8 and i 
consists of a cone-cylinder fuselage with four rectangu- _ ? fins . 
bir ee 3 ao eae TI a ee f 4.454 Aluminum 1.253 1.221 
ar fins in crucilorm arrangement. le fin surlaces fuselage 
are single-wedge-shaped with 16 per cent maximum owe 
. - . ~ yronze 
thickness. On each model, two of the tail fins were rs 
canted. 
For five of the models, the fuselage was of aluminum 
with wings of Tobin bronze. Bronze wings were used (B) Results and Discussion 
to obtai a highe : f inertia 1 i . > 4 cae kT 
tain (1) a higher moment of inertia in roll in order Four of the missile models (Table 1) were gun- 


to decrease the rate of acceleration in roll and corre- 
spondingly lengthen the curved portion of the measured 
roll displacement-distance curve, and (2) a higher mass 
in order to decrease the linear deceleration of the missile 
and result in more nearly constant-velocity flight. 
The remaining two models were machined entirely from 
Tobin bronze to increase further the moment of inertia 


launched at an intended Mach Number of about 1.7 in 
order to investigate reproducibility. The other three 
models were fired in order to obtain results at higher 
Mach Numbers. The actual Mach Numbers that re- 
sulted in the seven flights both at the beginning and 
end of the trajectory, as well as the atmospheric air 
density and atmospheric velocity of sound, are recorded 


in Table 2. 


TABLE 2 


in roll. A summary of the missile design is given in 
Table |. 
Mach Number 
at Begin Mach Number 
Missile ning of at End of 
Number Trajectory Trajectory 
l 1.790 1.675 
2 1.752 1.635 
3 1.621 1.575 
! 1.665 1.616 
5 2.722 2. 587 
6 2.301 2.166 


7 2.450 2.327 


‘ ~ 


Mach Number Velocity of 
at Center Sound 
of in Range 
Timing Atmosphere 


Air Density 
of Range 
Atmosphere 
Slugs per Ft.) 


Data (Ft. per Sec.) ( 

1.746 1,135 2.297 xX 107-3 
1.694 1,131 2.296 

1.592 1,121 2.357 

1.637 1,120 2.357 

2.682 1,125 2.342 

2.247 1,132 2.309 

2.416 1,132 2.311 
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7 | | | TABLE 3 
& | Model l 2 3 1 5 6 7 j 
sual | | , | P.E., deg. 0.3 0.7 0.6 0.6 0.7 0.3 
eit rt id 
ERR ae | | equation represents the actual motion of the models 
4000}-+ ) ! tested to within the experimental error. It therefoy, 
4 | rir appears that the rolling motion is unaffected by th 
1 | ee yawing motion for the mean square yaws encountered 
w 3200 io in the tests (see Table 4). 
3 ng ce The optimum values of the constants in the roll equa 
uw : ae tion were obtained by the methods outlined 1n the Datg 
os saa ? eo rit Reduction Section and are given with their P.E.’s jy 
us ‘Mie e% ‘ eae Table 4 as determined from the reduction. : 
z mon as Tt a: (2) Aerodynamic Coefficients.—The experimental 
4 —ss Lt ae values of the aerodynamic coefficients for the models | 
& ee a / | | = ; tested may be obtained from the data in Table 4 and 
soo} | a wes Eqs. (27). The theoretical values of the coefficients 
a ’ t+ for the models are obtained from Eqs. (23), (26), and | 
¥ a; ; ae | | | (27). The values for the coefficients are given in Table | 
0 Litt | Lt 5, and their good quality is again indicated by the sta- | 
™ = 7 wie ed _ tistical factor, P.E., the significance of which is dis 
POE LS TAMEE, €-FCEy cussed in the section on ‘‘Method of Reduction.” 
Fic. 9. Measured values of the missile roll angle at the range 


station locations plotted along with the curve of the derived 
equation. (The graph does not show the quantitative agree- 
ment of the experimental data to within +1° of the derived 
curve because of the insensitivity of the seale 


(1) Rolling Motion.-The angular orientation in 
roll of the first four models, as measured from the 
shadowgraphs taken at each spark station along the 
trajectory, is shown by the points in Fig. 9 to illustrate 
the data obtained. In this same figure are plotted the 
equations of the curves that best fit the data as obtained 
by the methods discussed in the previous section. The 
agreement of the data to within a degree of the fitted 
curve cannot be shown in the figure because of the in- 
sensitivity of the scale. 

The difference between the curves for missile Nos. 2 
and |, as well as Nos. 4 and 3, is a result of the differ- 
ences in the angles of cant employed. The differ- 
ence between the curves for missile Nos. 2 and 4, as 
well as Nos. | and 3, is due primarily to the differences 
in the roll moment of inertia (see Table 1). 

The use of a greater moment of inertia, it was felt, 
would result in a better determination of the roll co- 
efficients because of the slower roll acceleration rate 
and linear deceleration rate. However, no such vari- 
ation in experimental accuracy with moment of inertia 
can be noted in the results, and the accuracy appears 
dependent only upon the number and distribution along 
the trajectory of photographs successfully obtained dur- 
ing a given flight. 

The excellent ‘‘fit’’ of the roll equation to the experi- 
mental points that cannot be observed in Fig. 9 because 
of size limitations is given in Table 3 as the statistical 
quantity, P.E. P.E. is small and of the 
same order of magnitude as the estimated error in roll 


Since the 


angle measurements, it may be concluded that the roll 


The experimental values of C,, and Ci, (Table 5 
are plotted, respectively, in Figs. 10 and 11 against 
Number (Table 2). The 
5, and 6 deviate 


the average flight Mach 


figures show that values for models 3, 


somewhat from the experimental curve. In particular 
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Fic. 10. Experimental values of Cy as a function of Mach 


Number plotted along with the theoretical curve. 
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SUPERSONIC FREE-FLIGHT TESTS OF A 
TABLE 4 
Model l 2 3 1 
A 660 1,386 S82 2,285 
3 —624 — 1,298 — 838 —2,079 
C 0.0166 0.0169 0.0096 0.0098 
P.E., % 0.2 0.3 2.3 0.5 
P 14.90 32.02 15.57 32.64 
P.E., % 0.1 0.1 0.8 0.2 
Q=Cs 0.247 0.542 0.150 0.321 
P.E., % 0.2 0.3 2.4 0.6 
Kp/m 3.7 x 2 2.8 1.2 1.2 
P.E., % 0.7 0.1 LJ 0.4 
VAN 1.48 0.34 0.35 1.15 
No. of observations 14 2] 14 18 
TABLE 5 
Model l 2 3 } 
Mach Number 1.746 1.693 1.592 1.637 
—(), (theory) 0.449 0.460 0.483 0.472 
— Cj, (exp.) 0.548 0.559 0.546 0.564 
P.E., % ‘ 0.2 0.3 2.3 0.5 
Cis (theory) 0.716 0.740 0.794 0.769 
Cis (exp.) 0.778 0.806 0.784 0.802 
P.E., % 0.2 0.3 2.4 0.6 
p b . 
oy (theory) 0.0283 0.0602 0.0303 0.0633 
psb = ss : _ 
oy (exP-) 0.0251 0.0540 0.0264 0.0553 
P.E., % 0.04 0.05 0.6 0.2 
the discrepancy in the value of model 3 is noted. obtained 


From Table 5, however, the P.E. associated with the 
experimental results for models 3, 5, and 6 are three 
or four times greater than the largest P.E. of any of the 
This is positive evidence that the 
reductions for these three The 
weakness is due not only to the relatively small number 
of photographs obtained during flight down the range 
but also to the poor distribution along the range of the 
photographs. In 


other four models. 


models were weak. 


stations yielding these successful 
particular, in Fig. 9 for model 3, note the lack of sta- 
tions at the beginning of the range where data are re- 
quired for a good determination of the constant C). 
The experimental curves in Figs. 10 and 11 are drawn 
through’ the weighted average the coefficients 
(weighted on the basis of the P.E.’s) for models | 
through 4 and very near to the coefficient for model 7. 


of 


This procedure, then, enables the curves to be drawn on 


the basis of the statistical accuracy of the data 
points. 
The theoretical values of C,, calculated from the linear- 
)) 


ized theory as applied in this paper are 16 to 22 per cent 
less than the experimental values, and the theoretical 
values of C), are 7 to 17 per cent less than the experi- 
mental values. 

Concerning this deviation of the experimental values 
from the theoretical, it will be remembered that the 
airfoil section used was of triangular shape (single 
wedge) with about 16 per cent maximum thickness at 
the trailing edge. A comparison of the theoretical lift 
of an infinite-aspect-ratio wing of such cross section 


from the 
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ROLLING MISSILE 619 
5 6 7 
918 925 1,930 

—919 — 900 —1,776 
0.0124 0.0127 0.0128 
1.5 2.0 0.4 
14.17 15.78 33.58 
0.4 0.9 0.1 
0.176 0.200 0.429 
1.6 2.2 0.4 
1.8 2.2 1.9 
1.2 0.1 0.7 
1.69 1.59 0.09 
14 14 19 

5 if) 7 

2.682 2.247 2.416 

0.310 0.362 0.340 

0.401 0.426 0.424 

1.5 2.0 0.4 
0.461 0.549 0.511 
0.570 0.618 0.619 
1.6 2.2 0.4 

0.0251 0.0277 0.0583 

0.0240 0.0261 0.0559 

0.4 0.9 0.1 


by the linearized theory with that obtained 
exact theory using the oblique shock equa- 
































ROLL MOMENT DERIVATIVE OF A SURFACE DUE TO CANT Cy 
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MACH NUMBER, M 


Experimental values of Cig as a function of Mach 
Number plotted along with the theoretical curve 
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STEADY STATE FIN-TIP HELIX ANGLE EFFECTIVENESS PARAMETER os , RAD. 


° 10 20 30 40 
CANT ANGLE 6, DEG. 
Fic. 12. The effectiveness parameter plotted as a function of 


the angle of cant on each of two fins. 


tions results in about a 15 per cent lower value for the 
linearized lift at a Mach Number of 1.7 and about 22 
per cent at Mach Number of 2.6. 
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CANT ANGLE 86, DEG. 


Fic. 13. The roll moment produced by a canted fin surface 
plotted as a function of the cant angle. 


SCIENCES—OCTOBER, 1950 
TABLE 6 
Mach ' 
Model Number G., Cis (p.b/2V)/5 & 
l 1.665 —(0.566 O.820 0.01296 
2 1.665 —0.565 0.821 0.01956 
3 1.665 —0.580 0.746 6.01261 
4 1.665 —Q0.557 0.787 0.01245 
Weighted average —0.564 0.816 0.01238 
P.E. of weighted 
average, % 0.3 0.6 0.5 
Theory —0.466 0.758 0.014 
(3) Effectiveness Parameter.—The — experimental 


values for the effectiveness parameter obtained from 
Eq. (33), Table 
and the theoretical values for this parameter obtained 
from the Eqs. (23), (26), and (34a) are given in Table 5, 


4, and the physical measurements 


together with their P.E.’s. Fig. 12 gives the theoretical 
curve for the effectiveness parameter per degree angle 
of cant (p,b/2V)/6 as a function of Mach Number, to- 
gether with the experimental points. It is noted that 
the experimental curve is approximately 12 per cent 
lower than the theory and does not exhibit the upward 
curvature in the low Mach Number region. 

(4) 
Data, and Linearity. 
producibility of the data was estimated by reducing the 
experimental values of Cy, Ci, (bsb/2V)/6 for the first 
1.665 § 
1.669. 


Internal Consistency or Reproducibility of the 
The internal consistency or re- 


four models to a common Mach Number of J/ = 
by the theoretical value of the slopes at J/ = 
These reduced values are given in Table 6. 


The weighted average obtained from reference & is 


> C/P.E.? 
>1/P.E.? 


= weighted average 





and the P.E. of the weighted average, reference §, is 
' weg [W117 + Weto” + W3d3? + WyV4" 
P.E. = 0.6745 
(n — 1) (Wy Ws W3 + Ws) 
where 
w= 1/P.E.? 


difference of coefficient from average value at 


= 
same Mach Number 


n = number of observations 


These values of P.E., given in Table 6, are a measure 
of the internal consistency or reproducibility of the data, 
and their small values indicate that the technique 1s 
promising. 

In plotting Ci, vs. J, it was assumed that C, was 
linear with 6; this assumption may now be verified by 


plotting C, vs. 6, where 


It is seen in Fig. 13 that the assumption of linearity 
is consistent with the experimental results. In Fig. 
14, the effectiveness parameter is plotted vs. 6, and it 


is seen to be linear in good agreement with theory. 
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5. 








SUPERSONIC 


CONCLUSIONS 


(1) The tests establish the validity of the roll equa- 
tion for accurately representing the true motion of the 
missile, as well as the accuracy of the reduction meth- 
ods employed, since the P.E. of the residuals of the 
“fitted” equation is a maximum of 0.7° and is of the 
same magnitude as the estimated error in the roll angle 


measuremen L. 


(2) The tests show that the rolling moment coeffi- 
cients due to cant and the effectiveness parameter de- 
pend linearly on the angle of cant for the range of angles 
employed in the tests. 

(3) The tests show that the absolute theoretical 
values of the coefficients C;, and C;, as obtained for the 
16 per cent thick wing are, respectively, from 16-22 
and 7-17 per cent less than the experimental values 
over the range of Mach Numbers 1.6 to 2.7. 


(4) The tests establish the good reproducibility of 
the aerodynamic coefficients as shown by the weighted 
P.E.’s of 0.3 per cent for C;,, 0.6 per cent for C;,, and 
0.5 per cent for (pb/2V)/6 when the values of the 
coefficients are reduced to a common Mach Num- 


ber. 


AppENDIX I—ROLLING MOMENTS DUE TO ANGLE OF 
ATTACK a@ AND ANGLE OF YAW B 


These two moments are identical in nature because 
of the symmetry of the missile surfaces. The magni- 


tude of such a moment is given by 
L = (0C_,/0B)B(pV?/2)Sw. .b 


where (OC,/08) is termed the stability derivative and 
Sw, 1s the area of the wing, tail, or both. Then, 


FE (’ —B? 3+ m)| v pV? _ 
= — —— : tee 
B? AB 3A°B? y Z 


where the expression for 0C;,/08 has been obtained from 
reference 9 and 8 now is replaced by the ratio v/V, 
where v is the sideslip velocity. In Eq. (I) A is de- 
fined as the aspect ratio of the surface considered. 
Now, when a, the angle of attack, is replaced by its 
equivalent w, 1’, where w is the component of missile 
velocity in the vertical plane, Eq. (I) becomes 


| | (’ —-B? 3+ ~) | dS... 
L= — me (vw) 
B?*\ AB 3A °B? 2 


But (vw) for the motions resulting from these experi- 
mental tests is the product of two small disturbance 
velocities compared with V’, and, as such, the moment 


(II) 


is of second order and negligible. 
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STEADY-STATE FIN-TIP HELIX ANGLE PER DEGREE CANT ANGLE 


14. The effectiveness parameter (p,b/2V)/6 plotted as a 
function of the flight Mach Number 


_ 
= 


REFERENCES 


1 Charters, A. C., Some Ballistic Contributions to Aerodynamics, 
Journal of the Aeronautical Sciences, Vol. 14, No. 3, p. 155, 
March, 1947. 

2 Charters, A. C., and Thomas, R. N., The Aerodynamic Per- 
formance of Small Spheres from Subsonic to High Supersonic 
Velocities, Journal of the Aeronautical Sciences, Vol. 12, No. 4, 
p. 468, October, 1945. 

3 Bolz, R. E., and Nicolaides, J. D., A Method of Determining 
Some Aerodynamic Coefficients from Supersonic Free-Flight Tests 
of a Rolling Missile, B.R.L. No. 711, April, 1949. 

‘ Liepmann, H. W., and Puckett, A. E., Aerodynamics of a Com- 
pressible Fluid, Chapt. 9; John Wiley & Sons, Inc., New York, 
1947. 

5 Evvard, J. C., Distribution of Wave Drag and Lift in the Vicin- 
ity of Wing Tips at Supersonic Speeds, N.A.C.A. T.N. No. 1382, 
1947. 

® Snow, R. M., and Bonney, E. A., 
istics of Wings at Supersonic Speeds, Bumblebee Report No. 
55, The Johns Hopkins University, Applied Physics Labora 
tory, March, 1947. 

7 Whittaker, E. T., and Robinson, G., The Calculus of Observa- 
tions, Blackie & Sons, 1926. 

8 Scarborough, J. B., Numerical Mathematical Analysis; The 
Johns Hopkins Press, 1930. 

® Harmon, S. M., Stability Derivatives of Thin Rectangular 
Wings at Supersonic Speeds, N.A.C.A. T.N. No. 1706, Novem- 
ber, 1948 


Aerodynamic Character- 











Base Drag and Thick Trailing Edges 


SIGHARD F. HOERNER* 
Air Materiel Command 


ABSTRACT 


The negative pressure originating behind the flat base of pro- 
jectiles presents a drag component that is consequently termed 
“base drag.”’ As illustrated in Fig. 1, the same kind of drag 
exists for every elongated body the flow around which is sepa- 
rated from its rear side, such as airfoil sections, the trailing edges 
of which are cut off or thickened. The same type of drag is also 
found behind certain sheet-metal joints. 
late all the various base drags with each other by considering the 
effect of the boundary layer originating from the forebodies. 
Analyzing competent measurements, simple formulas are found 
with which to calculate the drag of tail turrets and blunt trailing 


It is possible to corre- 


wing edges. In the case of thicker sections, such as are used for 
the blade roots of propellers, the lift is also affected, and the lift- 
drag ratio can be improved remarkably by thickening the trailing 
At supersonic speeds the base drag approaches the nat- 
Depending upon 


edge. 
ural limit as given by zero pressure (vacuum), 
the body shape, the base pressure assumes an approximately 
constant percentage of the theoretical limit. 


SYMBOLS 


coordinate across direction of flow 


y oe 
l = length of bodies in direction of flow 

C = wing chord 

h = height or thickness of trailing edge 

d = diameter of the base or supporting rod, respectively 
D = maximum body diameter 


thickness of the boundary layer 
v local velocity within boundary layer 


V = velocity of flow or flight 

q = dynamic pressure 0.5p V? 

PI = area to which the drag coefficient is referred 

Cy coefficient indicating drag due to skin friction and 
surface protuberances 

Cpr = base drag, based upon base area 

Cy = body drag, based upon maximum cross-section area 

Cpo = profile drag coefficient, based upon wing area 

Ap = static pressure difference 

B = subscript indicating the base 


In several figures, a ‘“‘dot”’ is used as a subscript, indicating the 
drag coefficient based upon the frontal area. Cpp is eventually 
used instead of Copp. 


(1) BASE DRAG OF THREE-DIMENSIONAL BODIES 


(A) Mechanism of Base Drag 
B* EXAMINING THE FLOW PATTERN around and be- 
hind the base of a projectile (Fig. la), it is found 
to be that of a jet pump. The jet—that is, the outer 
flow——placed like a tube around the space behind the 
base, mixes with the dead air and tries to pump it away. 
However, since there is no supply from which additional 
air can be drawn, the pump has no delivery, and its 


Received October 17, 1949. 
* Dr.-Ing. habil. Wright-Patterson Air Force Base. 
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only effect is to reduce the static pressure on the body 
base. . 

The base drag of projectiles is often quoted to ly 
around Cp, = 0.2, equivalent to a base pressur 
Aps/q —0.2. However, by studying available 
pressure or drag measurements, the base pressure jg 


recognized to depend largely upon the length of th 
forebody, the ratio of base to body diameter, and th, 
surface conditions of the projectile. 

The base pressure is proportional to the dynamic 
pressure of the outside flow. Evidently, however, th 
boundary layer existing at the end of the forebod, 
serves as an insulation sheet, reducing the effectiy; 
dynamic pressure of the jet-pump mechanism. The 
boundary-layer thickness is proportional to the drag 
originating from the surface of the forebodyv. Conse 
quently, in Fig. 2, the drag coefficient 


-— D xin qd Sz = CAS Sp) 


é SB 


is chosen as the parameter upon which the base drag is 
considered to depend. In this equation, C;, indicates 
the skin drag of the forebody, which includes not only 
the friction drag but also any additional pressure-drag 
components originating from the surface, such as the 
drag of driving bands or other protuberances. 

In supersonic cases, the wave drag is not to be in- 
This kind of drag generally does not 
noticeably the thickness, 
since the equivalent momentum is generally carried 
Ot course, 


cluded in Cy. 
increase boundary-layer 
sideways by the respective shock waves. 
there remains some reduction in dynamic pressure be- 
Between ./ = | and 2, however, 
that is, in the order of less 


hind a shock wave. 
this loss is relatively small 
than 10 per cent—when considering projectile bodies 
with pointed and slender noses. Therefore, by evalu- 
*l the wave drag was en- 


ating supersonic tests,'® * 
tirely neglected. 
The available experiments on projectiles and fuselages 
(Fig. 1b), as plotted in Fig. 2, show approximately 
Cop = 0.029, V Cre (2 
As far as the tests were extended, the base pressure 
varies between Ap;/g = —0.03 and —0.30. From 
analysis, the thickness of the boundary layer at the end 
of a projectile is found to be, roughly, 
s/D = Crp (5 
valid for small drag coefficients. Consider, for example, 
the case with Cy, = 1, with the boundary layer extend- 
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Several cases of base drag. 
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ric. 2. Drag originating from the base of three-dimensional 


bodies, depending upon their skin drag. 


ing outside and around the base for a thickness equal 
to approximately two-thirds of the base diameter. 
Here, it seems that the effect of the ‘‘jet pump’’—that 
is, the amount of negative base pressure—is reduced to 
one-tenth of that near to Cy, = 0. 


(B) Effect of Boat-Tailing 

In the case of bodies, such as in Fig. 1b, where the 
base diameter is smaller than the maximum diameter of 
the body, the mechanism of the base drag remains essen- 
tially the same as in the case of Fig. la. Referring the 
base drag to the maximum cross section of the body 
and with C,, = Cp (D ‘d)?, from Eq. (2), 


THIicE 


TRAILING EDGES 


AC, = (5) C. = oom (5) 
_— D eS \ Cp D 


with C, indicating the original skin drag of the body 
withd— 0. Fig. 3 presents measurements on a smooth 
and slender streamline body, the trailing end of which 
Fig. 4 shows the drag increment 
Curves, 


(4) 


was cut off in steps.° 
ACp, depending on the diameter ratio d/D. 
according to Eq. (4), seem to interpret the measured 
values sufficiently well. 

All this consideration of boat-tailing assumes, how- 
ever, a flow pattern attached to the sides of the tail. 
In cases where the tail is too short and the flow is 
separated from its flanks, the whole cross-section area 
of the body has to be considered as base and the proper 


drag has to be determined according to Eq. (2). Round- 




































































0.14 hi ' 
G. LCN Sb 2 Qik) tad Coso-oce 
012 AC, = 09 (A) 
6.10 
0.08 }— — —- 
ow 981 AG... “Suh | @ DVL Body R=2-lc", Ref 
0.06 - __—_+—_+— ae oe ——+——_{ 
0.04 :3 7 _ _— a — —_—— 
+—Af-—- Ep 
0.0% pee] pt > 
= FD 
0 | | | a nee l Nl | 
0 02 o4 0.6 0g aa 10 
Fic. 3. Drag of a fuselage body, the end of which is cut off. 





PI re _4, Fuselage Body, <=? D, Cp,= 0.05 
y e— = -f x from Pressure Distrlbu Aon | Ret 2 


+ from Force Meaturemen 


4 Tube with Lege =l0D, G,=0.1 
yA Pressure Measurements Ref 3 








625 j|—_ -—_—_ - — Equation (4): 
Daase .029 3 
Monat | | 
0.20 Pe 





= 


AC,, C.* os | 
ONS 7 / / 



































0.10 x 7 
¥ AZ 
0.05 a SH 
x AY’ 
‘A 
0 aes 





0 a2 0% 06 O3 40 "= 


Fic. 4. Base drag behind boat-tailed bodies of revolution. 








AERONAUTICAL 











624 JOURNAL OF THE 
0.20 @ Pendulum Method at R= 4 lov 
‘ se © Force Measuvements R= 4-10 
C a i. 
De Tentative Calculation with: 











i 
Ol5 . N Gy, =-0.2 [|-"] 
* ~ 
Due to Change in Flow Pattern 
alo o— > 
~ \ 


\ 
S ~ N 


i Stns Friction = 


0.05 











? 





08 d 10 
D 


Fic. 5. Drag of spheres at Reynolds Numbers above the critical, 
depending upon the diameter of the supporting rod. 


oc 

















0 02 oF 06 


ing the blunt end of a fuselage, for instance, affects the 
base drag only as far as the round edge is able to lead 
the flow somewhat farther into the dead space. Refer- 
ence 2 indicates a reduction of base drag coefficient 
(Cpp = 0.04 to 0.06) by ACpz = 0.005 to 0.010 upon 
adding a half sphere to the cut-off end of a fuselage 
This result applies, however, to a smooth body 
only. If one considers a real airplane fuselage with a 
coefficient two or three times as high (Cp = 0.1 to 0.2), 
the effect of rounding the blunt end is expected to be 
hardly noticeable. 

Besides the experiments in reference 2, there are 
also some points included in Fig. 4 measured behind a 
tube. A collar was attached to the rear end of this 
tube, and pressure measurements were made within the 
dead space behind the collar. The tests show that the 
theory, according to Eq. (2), also holds true in this 
special case with d/D > 1. 

The fuselage of a military airplane has a drag coeffi- 
cient of the order of Cp = 0.15 when accounting for 
the usual surface irregularities and protrusions. The 
base drag caused by cutting off or thickening the trail- 
ing end amounts to AC, = 0.09 for d/D = 1 (case like 
In order to accommodate a tail turret, how- 
0.5 may be sufficient. 


body. 


Fig. la). 
ever, a diameter ratio d/D = 
The corresponding amount of drag AC» is less than 
0.01, which is comparatively small—that is, in the 
order of 6 per cent of the fuselage drag. 


(C) Effect of Supporting Rod 


In order to support projectile models in wind tunnels, 
a sting is usually used leading from downstream to the 


SCIENCES—OCTOBER, 1950 
base. 
Fig. 5, it is known that such a rod affects the drag re. 


From measurements on spheres, as plotted jy 


markably as far as Reynolds Numbers above the critica 
are concerned. At first, one assumes that the drag jg 
reduced because a portion of the rear surface of th 
sphere and the corresponding amount of base drag js 
covered or cut out by the rod. A tentative caleyly. 
tion, however, accounting for this effect does po 
check the actual Evidently, the 
flow pattern of the sphere is changed because of the 


drag reduction. 


addition of the rod. The cross section of the “jg 
pump,’’ as described above, is changed from a circular 
to a ring-shaped form, and, as outlined later, the puny 
effect is thus improved. Consequently, the partial 
flow separation still existing behind the sphere at super. 
critical Reynolds Numbers and the corresponding pres 
sure drag are reduced. 

The base drag behind projectiles and similar bodies 
is different from that of the sphere; the separation js 
fixed at the sharp edge of the base. Starting from 
d/D = 1, a ring-shaped base area is formed upon ée. 
creasing the rod diameter. The corresponding flow 
pattern can be considered as a two-dimensional one, 
and the drag can be properly calculated from Eq. (9), 
At d/D = 0, however, Eq. (2) for three-dimensional 
base drag has to be applied. Between this point and 
the two-dimensional curve one has to guess as to the 
course of the base pressure. Tentative calculations, 
the thickness of the boundary layer and the effective 
dynamic pressure of the “‘jet pump’’ being considered, 
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Fic. 6. Variation of base pressure due to rear supporting rod 
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Drag of a sheet-metal joint; experimental points from 


Fic. 7. 

reference 5. 
suggest a slight increase of the negative amount of base 
pressure upon increasing the rod diameter beyond 
d/D = 0.1. Model tests, especially in high-speed and 
supersonic wind tunnels, usually employ rods heavier 
than according to this permissible ratio. 

Considering the base-pressure coefficient AC» in the 
upper part of Fig. 6, the variation is not so great as that 
of the pressure coefficient. Contrary to Fig. 5, how- 
ever, there is no decrease of base drag to be expected as 
far as diameter ratios smaller than 0.5 are concerned. 

(2) DRAG OF SHEET-METAL JOINTS 

The flow pattern behind the sharp-edged sheet-metal 
joint, as shown in Fig. le, is similar to that past the 
thick trailing edge of an airfoil section (Fig. ld). In 
both cases, the corresponding pressure drag is, by prin- 
ciple, the same kind as that originating from the base of 
a projectile. 

The general distribution of the velocity across a 
turbulent boundary layer is 


(d) 


The drag of a small body, placed within the boundary 
layer, is proportional to the effective or average dy- 
namic pressure g,, corresponding to its height h. 


q./q = 0.7 (h/s)” (6) 


Consequently, in Fig. 7, the drag of sheet-metal joints, 
similar to that in Fig. le, follows more or less accu- 
rately the '/; power of / as far as the range of hs < 
| is concerned. 
Cy = 0.18 Vh/s (7) 

This trend is proved by many other measurements on 
various surface protuberances.° 

With the velocity distribution as given in Eq. (5), the 
skin-drag coefficient of one side of a plane plate is 
approximately 


C; = 0.2 ts, l); Crp = 0.2 (s h) (S) 
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when considering the “‘base”’ of a sheet-metal joint. By 
introducing h/s = 0.2/Cy, into Eq. (7), the following re- 
sults: 


Con = ().] fe (9) 


This function is plotted in Fig. 8, together with the ex- 
perimental points as taken from Fig. 7. 


(3) CHARACTERISTICS OF THICK TRAILING EDGES 


(A) Drag at Trailing Edge of Two-Dimensional Bodies 


Fig. 8 presents the available measurements on the 
additional drag, originating from airfoils, the trailing 
edges of which are either thickened (Fig. Ic) or cut 
off (Fig. 1d), thus forming a surface equivalent to the 
base of projectiles. Considering such tests, the value 
Cy» is found from 


Crp = 2 Cf(c/h) (10) 


with C, indicating the skin drag proper to all of the body 
in front of the base—that is, including any superve- 
locity effects and the drag of surface roughness and 
For the experiments at 


protuberances. most of 


smaller Cygz values, 


Con = 0.135/WCyp (11) 


At the same coefficient Cyg, if one compares Figs. 2 
and 8, the pressure difference behind the trailing edge 
of airfoils appears to be at least three times as high as 
the base pressure of three-dimensional bodies. This 
difference can be understood upon studying the veloc- 
ity distribution across the wake behind the body. 
Assume a boundary layer, with a distribution as given 
in Eq. (5), as arriving at the edge of the base, and the 
character of this distribution may be more or less main- 
tained upon entering the space behind the base, as far 


as the two-dimensional case is concerned. In the three- 
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dimensional case, however, the boundary-layer jy, “ig 
a) 2-Dimensional Plate cage ; ri , det <F mate. ! é 
sim - rial is transferred from a ring-shaped to a Circular eros. inc 
Cy,= 2.0; FP =-13 section. In doing this by velocity and volume, jt jg 
nee found by mathematical means that the power in Eq 
(5) is changed from '/¢ to '/3. Consequently, the djs. Th 
tribution of the dynamic pressure is changed from the oti 
= '/3 to the '/;.5 power of the distance or the radius, re. (9) 
Cp.= 1.5 spectively. Confirming this trend, the power in Eq 
(2) is 1/5 as compared to '/; in Eqs. (9) and (11) CL 
C) 3-Dimensconal Plate ov Disk s d | ? . r ys ; At . 
nN the same time, the magnitude of the average dynamic Cp 
= si i pressure within the boundary layer is reduced tg 
21.2; SPB=-0. . ° . “i : ° . 
Gah 4i a roughly one-third of that in the two-dimensional case. 
Fic. 9. Flow pattern past various plates. In (b),Aps/g = 0.8. | Most of the difference in base drag is thus sufficiently 
explained. . 
x: . ° a ‘O 
In Fig. 8, most of the measurements on airfoil sec. y 
; : sla 0. 
tions show somewhat higher drag coefficients [Eq ; 
+ > Ee. re! 
0020 (11)| than the sheet-metal joints [Eq. (9)]. This : 
J 1 di 
¢ 4c, = (0% 0.135 difference is probably attributable to the flow pattern, ; 
3 van oe . ? , roa ab 
Op , tre, The mixing of the outer flow with air from the dead 
. col 
00S | Pa space behind the base does not proceed so smoothly as 
: ae "iis ne 
0 A illustrated in Fig. la. Instead, the pattern is more ' 
° ° ° ° ° ° ° O 
A = (4)% 6.10 irregular, forming bigger vortices, intruding deeper 
a te SC) . f , 2 wi 
a | : (om into the space once in a while, and carrying away air 
0.010 | particles in lumps. Behind two-dimensional bodies, 
. S10 
the flow usually assumes the regular stable pattern of 
x DVL 0018, R=2.6:105 Ret 8 . hie 7 aed semen wane mi 
Munich 2415, R=210% Ree 9 the von arman vortex-street, as illustrated in Fig. 9a. “a 
Because of the vigorous and extended up-and-down 
0.005 movement of this street, the exchange of momentum - 
between the outer flow and the air particles in the dead th 
space is intensified. On the other hand, in the case ol 
of the sheet-metal joints, any larger vortex movement | 
OQ J - § dr. 
S sce 72 , ‘ > sec c cca y > we 
0 / 2 3 4eh SF is slowed down and more or less suppre ssed by the wall | ‘ai 
Cc along which the movement takes place. Experiments - 
Fic. 10. Profile drag, at C, = 0, of airfoil sections, the trailing jn a water tunnel!’ prove that in the absence of the 
edges of which are cut off or thickened. : : ; pe 
vortex street—that is, with a flow pattern as illustrated a 
in Fig. 9b—the drag coefficient of a plate is reduced th 
from approximately 2.0 to 1.5. This means a reduc- TI 
09 tion of the negative ‘‘base’’ pressure from Ap/q * ta 
—1.3 to —0.8. Approximately the same ratio is found oe 
a in Fig. 8S between the respective two curves (constants ; 
8 Split- Flap Tests j 8 I ° ex 
° o Reference 15° 0.135 and 0.10). re 
, © NACA Ref IT as ae —" a 
a7 . ] Ihe base-drag coefficients in Eqs. (9) and (11) are pe 
e e . 
oy h=c,., sin converted into profile-drag coefficients by multiplying m= 
= __ flap ‘ ‘ | 
0.6 ee t them with h /c. , 
e | t 
of Pd ; 0.10 /h\” 0.135 (h\" in 
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Vv Cpoo Vv Choo c 
04 4 , al , 
rs oe a where Cpw = Cy (h/c) indicates the profile drag ol the 
4 
0.3 : al section for h/c = 0. Fig. 10 (like Fig. 8) shows a 
MC marr 2 |~ "a | transition of the experimental points from one of the i 
0.2 f x | calculated curves to the other. Obviously, for small ol 
. 4 DVLOOIS Ref 8 : : i alba ' ne 
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Junkers 0040 Ret! : fie 
Of s oman as roughly one or one-half times the boundary-layer thick- as 
a ness—a regular vortex street does not originate. = 
oO ef 
Oo i = 6 §& (0 12 MH” 6 (B) Effects of Trailing Edge on Lift 
c 


. ; : : ; a, In the case of wings, thickening the trailing edge not 
Fic. 11. Increase in maximum lift due to thickness of trailing 8S, g° § 5 ‘ 
edge. only affects the drag but also the lift. As evident from 
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Fig. 11, the thick edge works similar to a split flap, 
increasing the maximum lift approximately by 
ACz mez, = 4(h/c) (13) 


This effect is in the order of one-third of that of an 
jinary split flap deflected to 90°. Combining Eqs. 


ore 
(9) and (13), the efficiency of an airplane is given by 
Ci maz. ’ + 4 h ° 
= ( L mar. 6 : bs 
( Dmin . 


(14) 


C ‘ 0.1 (“) 
D min. 0 Pls c 


For conventional values, such as Cy, maz. = 1.4, Cp aia. = 
0.025, and Cpoo = 0.01, this equation has its maximum 
between h/c = (0 and 0.4) per cent. This means, 
with h/c = 0.4 per cent, the ratio Cy maz./Cp min. iS 
about the same as at h = (. Four-tenths per cent is 
considered permissible, and the corresponding thick- 
ness of the trailing wing edge of any airplane is found 
to be between !/, and | in., depending upon the size of 
wing and plane. 

The effect of a blunt trailing edge on the lift-curve 
slope is similar in relative magnitude to that upon the 
maximum lift.7 > This effect may be important in the 
case of control surfaces. 

Fig. 12 suggests another application of the thick 


trailing edge. The lift characteristics of a 40 per cent 


thick section are remarkably improved above C;, = 0.4 
when using a thick trailing edge; the maximum lift 
drag ratio is increased by some 100 per cent. Such 


sections 
are employed at the blade roots of propellers. 
pears to be favorable, in this case, to manufacture the 
trailing edge with h/t in the order of 20 per cent. In 
this way, the minimum drag is not necessarily increased. 
The flow past thick sections is usually somewhat de- 
tached from the trailing edge. Consequently, thick- 
ening this edge means filling out the dead space. The 
experiments in reference 14 lead one to assume a slight 
10 per cent or h/c = 4 
Such sections 


that is, profiles with high thickness ratios 
It ap- 


reduction of drag with 4/f = 
per cent, for a 40 per cent thick profile. 
may be especially advantageous in the case of turbo- 
prop engines, where the combustion air usually is taken 
in between and around the blade roots. 


(4) BASE DRAG AT SUPERSONIC SPEEDS 


At first there is no fundamental change to be ex- 
pected in the mechanism of the “jet pump" because 
The 
negative base pressure finds a natural limit, however, 
The so-defined 


ol approaching and exceeding the speed of sound. 


as the condition of vacuum is reached. 
maximum amount of base-pressure and base-drag co- 
efficient is a function of the Mach Number. 


" (= r? 1.42 
Cor nes. ) = oo . 
qd J max yM Mv 


(15) 


THRICE 


TRAILING EDGES 627 




































































42 | | 
te - Ciel ~ ; 
el 1 an 
10 [with 04+ Trailing Edge F if 
. 4 T s d 7 
i | ! 1 | 
* } . : 
of |} | fom died 
H y é ot 
06 - whith Thin Trailing Edge le 
> As, Resi] ~ I Jt ts 
4 wd m , } “ 4 
0: : o4t bY a +—- doy 
: rTP 
4 . . r | 
02 zt ice a 4 —ESS 
o|< | | | 
~ h al G 0.2 “10 0 10 20 a°3o 
i sof 
02} +} natin - : . * 
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I 


different thickness ratios of the trailing edge.'' =4- c= 


5.9 in.; Ress. = 5 XK 10. 


with the constant y = 1.405. As shown in Fig. 15, this 
upper limit is never likely to be encountered by any base 
pressure in the subsonic field. At supersonic Mach 
Numbers, however; the theoretical limit comes down 
remarkably. A series of base-pressure measurements 
carried out on projectiles and similar bodies suggests a 


certain percentage of the theoretical amount 


Con = 0.5 Cop mar. ~ 0.7/M* (16) 


Of course, this drag is true only for a certain boundary- 
that is, for a specific value Cp as used 
In the case of Eq. (16) the cor- 


layer thickness 
in Fig. 2 and Eq. (9). 
rect subsonic value of Cpg is in the order of Cpz = 
0.25. 

Experiments recently published in reference 19 con- 
firm the hypothesis as expressed in Eq. (16). Further- 
more, these tests are remarkable because they are the 
only known ones extending through the speed of sound. 


In the vicinity of 17 = 0.9, the tests show a minimum 
in the order of Cpg = 0.08. Considering the flow 
for this Mach Number, the low 


pattern 
base-pressure difference can be explained when as- 
and .\/ = 1 most of the 
Correspondingly, 


proper 


suming that between J/,,;:. 
drag is caused by flow separation. 
the effective dynamic pressure of the “‘jet pump’’ mech- 
anism is highly reduced. In order to account for 
Cpr = 0.08, it is sufficient to assume in Fig. 2 a fore- 
body drag-coefficient Cy, = 9.15. The flow pattern 
changes rapidly, however, when approaching J = 1. 
Here, the loss of dynamic momentum, as equivalent 
to the drag of the forebody, is evidently shifted from the 
boundary layer into the outer space. With most of the 
momentum being carried away from the body by means 
of pressure waves, the dynamic pressure is not remark- 
ably reduced to any extent around the base, and, as 
done before when evaluating the experiments of refer- 
ence | in Fig. 2, only the ‘“‘skin drag’’ of the forebody is 
to be considered as affecting the base pressure. 

When the subsonic base drag of any body is known, 
it is possible to estimate the course it will take when 


approaching the upper theoretical limit [Eq. (15)]. 
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In order to assist in such estimating, Fig. 14 was pre- 
pared. Including data as found for plates, the rear 
side of which can also be considered as base, three ex- 
perimental points are available. From these points, 
the zero point, and the upper limit as given by Cpg + 
Can ines. 1, a curve is defined which permits one to 
make reasonable estimates of the base drag at super- 
sonic speeds. 

Comparing the base drag to the wave drag as originat- 
ing from the rear of double cones and double-wedge sec- 
tions, limiting Mach Numbers are found beyond which 
single cones or wedges, with half the nose angle and 
flat base, are superior as to their total drag over the 
double shapes. For a 10 per cent thick wedge section, 
for example, this point is reached between \/ = 2 and 3. 


RESULTS 


(1) The base drags of projectiles and of fuselages 
with blunt or cut-off trailing end are a distinct function 
of the forebody drag (Figs. 2 and 8). 

(2) From the drag characteristics of sheet-metal 
joints (Fig. 7), it is possible to calculate the drag orig- 
inating from the blunt or cut-off trailing edge of airfoil 
sections (Figs. 8 and 10). 

(3) <A thick trailing edge is suitable to improve the 
lift and the lift/drag ratio of thicker airfoil sections 
(Fig. 12). Such sections have practical interest for the 
blade roots of propellers. 


(4) At supersonic speeds the base drag is related 
to the maximum possible pressure differential, which 
is that between vacuum and the ambient Pressure 
(Fig. 13). The supersonic drag behind any body can be 
estimated (Fig. 14) when the subsonic coefficient jg 
known. 
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Lateral-Control Devices Suitable for Use 
with Full-Span Flaps 


FRANCIS M. ROGALLO,* JOHN G. LOWRY,? ano JACK 


FISCHEL? 


Langley Aeronautical Laboratory, N.A.C.A. 


SUMMARY 


A brief history of the development of lateral-control devices 
for use with full-span flaps is presented with special 


suitable 
Recent data on the per- 


emphasis on spoiler-type controls. 
formance of spoiler-type controls on unswept wings at high 
subsonic speeds and a limited discussion of spoiler controls on 
swept wings are included. The several advantages shown for 
spoiler-type controls high lift through the use of full-span flaps, 
good aileron effectiveness throughout the speed range, good 
yawing-moment characteristics, low wing torsional loads, and 


the possibility of low control forces—appear to warrant their 


consideration for use on all types of airplanes. 

Application of spoiler ailerons as air brakes while simul- 
taneously functioning as ailerons is discussed, and their per- 
formance is shown to compare favorably with thet of other 


currently used air brakes. 
INTRODUCTION 


I 1941, THE N.A.C.A. published Baker's paper! 
covering the development of the full-span flap 
arrangement for the Vought Kingfisher, and, in 1946, 
the N.A.C.A. published a similar paper by Ashkenas* 
on the development of the full-span flap and aileron 
arrangement for the Northrop Black Widow. Each 
of these papers presented coherent stories of a par- 
ticular development that was carried through to com- 
pletion by a particular individual or group in a rela- 
tively short time. The N.A.C.A. story of its own 
work on lateral-control devices suitable for use with 
full-span flaps, however, covers the work of many in- 
dividuals during the last 20 years. Partly because of 
the different points of view of these individuals but 
mostly because of the ever changing picture of aero- 
nautics in general during these 20 years, the N.A.C.A. 
story of lateral-control devices suitable for use with 
full-span flaps is really several stories, not always 
closely connected. In fact, during part of this time 
the use of full-span flaps was not contemplated or 
was of second-order importance. Some of the devices 
developed, however, were found to be suitable for use 
with full-span flaps and are therefore included in this 
story. 

In 1931, the N.A.C.A. began a systematic wind- 
tunnel investigation of lateral control with special 
reference to the improvement of control at low air 


Presented at the Aircraft Design Session, Eighteenth Annual 
Meeting, I.A.S., New York, January 23-26, 1950. 

* Aeronautical Research Scientist, Stability Research Division. 

t Aeronautical Research Scientist, Stability Research Division. 

t Aeronautical Stability and Control Scientist, Stability Re- 
search Division. 
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speeds and at high angles of attack. Many different 
lateral-control devices were subjected to the same 
systematic investigation in the 7- by 10-ft. wind tunnel, 
and the most promising devices were tested in flight. 
Two of the devices developed in this program-——the 
retractable aileron and slot-lip aileron, shown in Fig. 
l—would have been suitable for use with full-span 
flaps. 

At the outset, no high-lift flaps were used, but, after 
the inception of the systematic research program, 
flaps for the reduction of landing speed or the increase 
of wing loading came into general use, and the research 
program was modified to include lateral-control de- 
vices on wings with partial-span and full-span flaps. 
The effect of flap span on allowable gross weight or 
landing speed is indicated in Fig. 2, where it is seen that 
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Fic. 1. Lateral-control devices on wings without flaps. 
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Fic. 3. Lateral-control devices on wings with full-span split 
flaps. 
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Fic. 4. Lateral-control devices on wings with full-span slotted 
flaps 


reducing the flap span from 100 to 50 per cent of the 
wing span, as might be required for conventional 
ailerons, would either reduce the gross weight by 22 
per cent or increase the landing speed by 12 per cent. 
Before 1935, the plain aileron with flap retracting 
ahead and the retractable aileron or spoiler located at 
80 per cent of the wing chord (see Fig. 3) were developed 
and proved in flight to be suitable for use with full- 
span split flaps.* Slotted flaps soon began to displace 
split flaps, however, and the program was enlarged to 
include such arrangements. The N.A.C.A. investi- 
gated a large number of different configurations in the 7 
by 10-ft. wind tunnel*~ and, in 1941, recommended 
that the arrangements shown in Fig. 4 be flight-tested. 
Subsequent flight experience did not disclose any serious 
difficulties, and it appeared that any of these arrange- 
ments would provide satisfactory lateral control on 
airplanes with full-span slotted flaps. 

In the work reviewed thus far, emphasis was placed 
on the control and performance characteristics of the 
aircraft at the lower end of the speed range, but, as the 
maximum speeds of aircraft increased, control problems 
were also encountered at the upper end of the speed 


range. The first of these was the difficulty of providing 

enough aerodynamic balance to keep the stick or whe¢| 
forces within the physical limits of the pilot. Although 
the provision of power boost was an obvious solution 
the possibility of getting light forces by using a thiy 
retractable aileron was apparent, and this arrangemey 
was found to be satisfactory in the Northrop Black 
Widow night fighter.2, The low wing torsional loa 
of this aileron arrangement on the Black Widow wer 
a welcome bonus, but this characteristic was not of th 
major importance that it assumes when thin wings 
especially those of high sweep or high aspect ratio, ay 
used. The problem of wing twist or low aileron rf 
versal speed has, in fact, become so great that de 
signers are now considering the use of spoiler-tyy 
controls for the improvement of lateral control at high 
speed, regardless of their effect on the low-spee 
control or performance characteristics. That they 
control devices allow the use of full-span flaps may ofter 
be considered incidental, as was true in 1931 when the 
N.A.C.A. began the systematic investigation of such 
devices. 

While the lateral-control research reviewed above 
was ill progress, spoiler-type devices were also being 
investigated for an entirely different purpose—air 
brakes. This parallel line of research would have 1 
connection with our present story if it were not for the 
fact that often a particular device may be used for 
either lateral control or air brakes. An investigation’ 
has been made, moreover, of the application of spoiler- 
type devices as air brakes while simultaneously func 
tioning as ailerons. 

In the present paper, a brief history of the develop- 
ment of lateral-control devices suitable for use with 
full-span flaps is presented with special emphasis on 
spoiler-type controls. Recent information on_ spoiler 
ailerons on unswept wings, including their  simul- 
taneous use as air brakes, and a brief discussion oi 
spoiler ailerons on swept wings are included. More 
complete information on all phases of the work may 
be found in the available references, many of which are 
listed. 


INVESTIGATIONS ON UNSWEPT WINGS AT LOw SPEEDS 


Most of the interesting low-speed test data on lateral 
controls for unswept wings with full-span flaps has been 
published by the N.A.C.A. as a collection,’ and some 
discussion of results is also included in an N.A.C.A 
summary of lateral-control research.!” A large number 
of configurations has been investigated, and several 
appear to be aerodynamically satisfactory. To give 
an overall picture of the work, geometric and aero- 
dynamic characteristics of the devices shown in the 
previous figures and problems encountered during the 
development of those devices will be discussed briefly. 
In the work leading up to the ailerons of Fig. 1, un- 
vented spoilers located well forward on the upper 


surfa 
ward 
have 
vices 
extel 
man 
TI 
facto 
A fe 
requ 
on a 
arise 
affec 
able 
aero 
but, 
hing 
arral 
thos 
W 
way 
(Fig 
shov 
had 
of tl 
port 
stud 
grac 
behi 
cent 
trac 
flap 
this 
the 
behi 
ailet 
tou 
(see 
spor 
effe 
ate 
slot 
defl 
roll: 
ach 
teri 
hau 
and 
of | 
fou 
fore 
aile 
Th 
apy 
stri 
7 
plu 
(see 





Viding 
"Whee 
hougt 
lution 
a thin 
ement 
Black 
loads 
’ Were 
of the 
Witlgs, 
O, are 
MN fe 
it de. 
type 
- high 
speed 
these 
often 
n the 
such 


bove 
Deing 

air 
e 10 
r the 
1 for 
tion’ 
viler- 
unce- 


slop- 
with 
5 on 
oiler 
nul- 
n of 
lore 
may 

are 


EDS 


eral 
eel 


yime 


ber 
oral 
rive 
*TO- 
the 
the 
fly. 
un- 
per 


LATERAL-CONTROL 


ace had prohibitive lag, but the vented or rear- 
1 were found to 


yore located spoilers shown in Fig. 
j It should be noted that the de- 
vices of Fig. 1 had a rather high drag because of the 
rnal hinges of the retractable ailerons and the per- 


have acceptable lag. 


xte 
said open slot of the slot-lip aileron. 

The devices of Fig. 3 were concluded to be satis- 
factory from flight tests* on a Fairchild 22 airplane. 
A fairly complicated flap-retracting mechanism was 
required to allow use of the plain aileron, however, and 
on a larger or faster airplane some difficulties might 
arise from the fact that the hinge moments are greatly 
affected by the movement of the flap. The retract- 
able aileron, on the other hand, had extremely small 
aerodynamic hinge moments regardless of flap position, 
but, as already mentioned, the drag of the external 
hinges was high. In the flap-down condition, both 
arrangements had yawing characteristics superior to 
those of plain ailerons on unflapped wings. 

When the retractable aileron was modified in such a 
way that the hinge axis was within the wing contour 
(Fig. 5), an investigation in the 7- by 10-ft. tunnel‘ 
showed that with the split flap deflected the aileron 
had no effect until it was projected more than 2 per cent 
a significant 
Tuft 
studies showed that, as a result of the favorable pressure 
gradient established by the flap, the flow was reattaching 
behind the spoiler at spoiler projections less than 2 per 


of the wing chord above the wing contour 
portion of its total projection range (see Fig. 5). 


cent wing chord. In the same investigation, a re- 
tractable aileron on a wing with a full-span slotted 
flap gave similar results (Fig. 5). It was found that 
this ineffective region could be greatly reduced by 
the addition of a properly shaped slot through the wing 
behind the spoiler. To achieve this slot action during 
aileron operation without disruption of the wing con- 
tour with the aileron neutral, the N.A.C.A. plug aileron 
(see Fig. 6) was designed.’ The powerful effect of the 
spoiler slot is indicated in Fig. 7, which compares the 
effectiveness of plug and retractable ailerons at a moder- 
ate angle of attack. The beneficial effect of the aileron 
slot increases with lift coefficient. The 
deflection of the aileron (Fig. 6) has little effect on the 
rolling and yawing moments but greatly simplifies the 


downward 


achievement of good stick-force or wheel-force charac- 
teristics. Ailerons of this type have undergone ex- 
haustive wind-tunnel investigations by the N.A.C.A.° 
and have been flight-tested on several different types 
of airplanes. From these investigations it has been 
found that the maximum lift, rolling, yawing, stick- 
force, and wing twisting characteristics of the plug 
aileron and slotted flap combination are satisfactory. 
The principal disadvantage of the N.A.C.A. plug aileron 
appears to be that its installation complicates the wing 
structure more than a trailing-edge aileron. 

To avoid the structural complication of the N.A.C.A. 
plug aileron, the plain and slot-lip aileron combination 
(see Fig. 4) was developed by the N.A.C.A.®!!~" and 
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is satisfactory if one is willing to accept two separate 
aileron systems with a changeover from one to the other 
as the flap is deflected. One of these aileron systems, 
moreover, has restricted down-travel and high up- 
floating tendency, a combination that further compli- 


cates the stick-force picture. 
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for a plug aileron and a flap-type aileron on an aspect ratio 5.76 
wing. 


The plain and slot-lip aileron development by the 
N.A.C.A. and the Vought development of the system 


1 


for the Kingfisher' had much in common, and the 


conclusions drawn are in good agreement. In addition 
to the development of the combination systems, in- 
cidentally, both the N.A.C.A.'§ and the Vought Com- 
pany' investigated the possibility of utilizing the 
trailing-edge aileron for control with the flap at various 
initial settings (drooped ailerons) and found that the 
reduction of aileron rolling-moment coefficient and the 
increase of adverse aileron yawing-moment coefficient 
with flap deflection caused a progressive deterioration 
of control that became prohibitive with the flap de- 
flected more than 20°. 
interconnected to offset the adverse aileron yawing 


Although the rudder might be 


moment in such a system, it seems likely that the avail- 
able lateral control would decrease with reduction of 
airplane speed rather than increase as it may with some 
spoiler-type controls, a more desirable variation. 

The plain aileron and retractable flap combination of 
Fig. 4 is another arrangement that was developed by 
the N.A.C.A. in an effort to retain the trailing-edge 
using full-span slotted flaps. Two 


aileron while 
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versions of the arrangement were investigated, as 
indicated by flap positions A and B. When the flap- 
deflected condition is such that the nose of the flap jg 
near the aileron hinge line, position A, this flap jg 
essentially a split flap, and the characteristics of the 
arrangement are very much like those of the corte. 
The portion of the 
span covered by the aileron has a high drag with the 


sponding arrangement of Fig. 3. 


flap extended, but this drag may sometimes be tolerated 
if a low-drag flap, such as a Fowler’ or an N.A.C.A. 
slotted flap,’ is used inboard of the aileron. Because 
of the difference of flap configuration at aileron and 
inboard locations, this arrangement has generally been 
called full-span duplex flaps. 

During the investigation of the duplex-flap arrange. 
ment, it became apparent that, if the flap were moved 
back to the wing trailing edge, position B of Fig. 4, 
the aileron would be a slot-lip aileron, a device already 
found to have good lateral-control characteristics, 
The flap, moreover, would then be essentially a Fowler 
flap, an arrangement known to have good lift and drag 
characteristics. Such a combination might be used 
over the entire wing span, resulting in some simplifica- 
tion relative to the duplex-flap arrangement. With 
these thoughts in mind, the N.A.C.A. investigated a 
full-span retractable flap in combination with a full- 
span aileron of small chord.'*:'* The arrangement was 
found to have good lift, drag, and control characteristics 
with the flap either retracted or fully extended, although 
a reduction of lateral control was experienced at inter- 
mediate flap positions. Obvious disadvantages of the 
arrangement are the required flap-supporting system 
and the rather high wing pitching moments resulting 
from flap deflection as compared with the pitching 
moments of other flap arrangements. '® 

From the briefly reviewed work on unswept wings at 
low speed, it appears that any of several different 
aileron arrangements will provide satisfactory control 
on airplanes with full-span flaps and that the choice 
for a given application will be determined by the overall 
design of the airplane and the individual preference 
of the designer. For high-speed airplanes, however, 
it will be shown that spoiler-type ailerons may have 
definite advantages over trailing-edge ailerons. 


INVESTIGATIONS ON UNSWEPT WINGS AT HIGH 
SUBSONIC SPEEDS 


Because the results of any investigation of flap-type 
trailing-edge ailerons would be applicable in the design 
of full-span flap arrangements in which such ailerons 
are used in the flap-retracted condition (see Fig. 4), 
most of the high-speed investigations of devices in- 
tended specifically for use with full-span flaps have been 
investigations of spoiler-type ailerons.*» *~** In some 
of these investigations, however, flap-type trailing- 
edge ailerons have been included for direct com- 


parison.** ° 
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Some investigations of spoiler ailerons at low speed 
given the impression that spoiler controls located 


have 
between 0.60 and 0.80 chord are ineffective at small 
projections in the high-speed conditions. The results 


of several high-speed wind-tunnel investigations,” 
however, show that the variation of effectiveness with 
projection at the higher Mach Numbers is almost 
linear (see Fig. 8). These high-speed investigations 
also showed that the spoiler-aileron effectiveness in- 
creased with increase in Mach Number (see Figs. 8 and 
9) and that about the same effectiveness was exhibited 
at all Mach Numbers by both plug and retractable 
ailerons (with flap retracted). Inasmuch as spoiler 
ailerons also produced a linear variation of effectiveness 
with projection at high angles of attack and an airplane 
having even moderate performance characteristics 
would be flying at high speed at low angles of attack 
and at low speed at high angles of attack, it would 
probably have rolling characteristics that vary almost 
linearly with spoiler projection throughout the speed 
range with flap retracted. Flight-test data obtained 
with two spoiler-aileron arrangements on the Black 
Widow P-61 airplane® ** (see Figs. 10 and 11) and with 
two spoiler-aileron arrangements on a modified Vought- 
Sikorsky OS2U-2 airplane,** in general, substantiate 
the wind-tunnel results. Both the flight and wind- 
tunnel results show that the rolling effectiveness of 
spoiler ailerons increased with increase in Mach Number, 
as contrasted to generally opposite effects obtained with 
flap-type ailerons (see Figs. 9 and 11). 

Only recently have the high-speed hinge-moment 
characteristics of spoiler ailerons been reported.” *! 
In these wind-tunnel investigations, generally the same 
effects and trends have been noted as were previously 
noted for spoiler hinge-moment characteristics at low 
speed*—namely, the hinge moments were generally 
somewhat irregular over the aileron projection range, 
varied directly as the area moment of the top edge of 
the aileron about the hinge axis, and could be altered 
somewhat by any of several possible small modifications 
to the aileron. As is the case with flap-type ailerons, 
however, power-boost systems may be used with 
spoiler-aileron arrangements and would obviate any 
necessity for special treatment of spoiler ailerons to 
obtain acceptable control forces. 

Data obtained on the P-61 airplane®* (Figs. 10 and 
11) show that the control forces of the spoiler aileron 
varied smoothly and in the proper direction with 
spoiler projection and air speed. These data and the 
wind-tunnel results* °4 also show the small hinge 
moments of the spoiler controls and the small or 
negligible effects of Mach Number on these spoiler 
characteristics, in contrast to the large hinge moments 
of flap-type ailerons and the adverse effects of increase 
in Mach Number on the hinge moments of flap-type 
ailerons. 

When considering the lateral-control characteristics 
of ailerons at high speeds, it is necessary to investigate 
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the wing pitching moments produced by the ailerons, 
because the resulting wing twist reduces the control 
effectiveness. For example, in an analysis of the 
twisting effects of flap-type ailerons on two fighter-type 
airplanes employed in World War II, it was found that 
wing flexibility was responsible for a 31 per cent loss in 
aileron effectiveness of one production model of the 
P-47 airplane and for a 65 per cent loss on the British 
Spitfire at an air speed of 400 m.p.h."” 

In an investigation of the section pitching moments 
of flap-type and spoiler controls,** it was shown that, 
for normal spoiler-aileron locations (about 0.70 chord), 
the spoiler pitching moments are much less than those 
produced by plain flaps of equal effectiveness. The 
twisting effects of these control surfaces are perhaps 
better illustrated in Fig. 12 by a comparison of the 
twisting moments about an assumed wing elastic axis 
at 0.35 chord produced by spoiler ailerons and flap-type 
ailerons at high speed on an aspect ratio 5.76 wing, for 
which rolling-effectiveness data of the rigid wing are 
presented in Figs. 8 and 9. These data show that, for 
equal rolling effectiveness on a rigid wing, a spoiler 
control located at 0.70 chord produced favorable 
twisting moments that were relatively unaffected by 
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on swept wings. 


changes in Mach Number and, on a nonrigid wing, 
would tend to augment the rolling effectiveness of the 
control, whereas a 0.20-chord plain aileron produced 
adverse twisting moments that became slightly larger 
with increase in Mach Number and would tend to 
reduce the aileron rolling effectiveness. These effects 
are such as to make the comparative rolling-effective- 
ness data of Fig. 9 appear even more favorable for the 
spoiler ailerons. 


INVESTIGATIONS ON SWEPTBACK WINGS 


The problems involved in designing ailerons for swept 
wings are the same as, but differ in magnitude and im- 
portance from, those on unswept wings. Thus, we 
have some or all of the following characteristics to 
consider: aileron effectiveness as measured by rolling 
moment produced by a given control on a rigid wing; 
variation of aileron effectiveness with speed; yawing- 
moment characteristics of the aileron; wing twist and 
the corresponding loss of aileron effectiveness; and 
last, but not least by far, the hinge moments of the 
control. Aileron effectiveness is of primary importance 
The effects of 
adverse yaw and wing distortion, however, are of greater 


on both swept and unswept wings. 
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importance on thin sweptback wings than on relatively 
thick unswept wings, where both have often been taken 
care of in design calculations by a simple 0.8 factor 
Thus, on swept wings either the adverse yaw effects 
at low speed or the wing distortion at high speed may 
be the determining factor in the choice of a lateral. 
control system. All of the factors must, therefore, he 
considered when laying out a new design, and jo 
control can be discarded just because one of the charac. 
teristics is poor, for the other characteristics may out. 
weigh this deficiency. 

The lack of experience or of a reliable design pro. 
cedure necessitated investigations of controls on yn- 
flapped swept wings to provide design data for the 
high-speed condition with little regard, as yet, for the 
As a result, this discussion 
Of the 
many devices already mentioned as satisfactory for use 


flap-deflected condition. 
will be limited to the flap-retracted condition. 


with full-span flaps (see Fig. +), only the spoiler provides 
control in the flap-retracted condition by other than 
flap-type ailerons. Some of the characteristics of 
spoilers will be discussed and compared with those of 
flap-type ailerons. 

Theoretical studies of rolling effectiveness, pb 2), 
of flap-type ailerons on swept wings” ** indicate that 
the rolling effectiveness for a given aileron configuration 
on wings of constant aspect ratio decreases with in- 
creasing sweep (Fig. 13). For comparative purposes, 
data from a few spoiler configurations” * are spotted 

, 


on Fig. 13. 
wing plan forms, the results indicate that spoilers of the 


Although these data are for different 


same span as the flap-type aileron will give lower 
effectiveness on rigid wings at low speed and will show a 
more marked effect with sweep than flap-type ailerons. 
These adverse characteristics of spoilers may be allevi- 
ated somewhat at high speed, since unswept-wing data 
show that spoiler effectiveness increases more rapidly 
than the effectiveness of flap-type ailerons with in- 
creases in speed. Further increases in effectiveness 
can often be obtained by moving the spoiler inboard or 
by changing its alignment.*” For example, rotating the 
spoiler so that it is normal to the free stream instead 
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LATERAL-CONTROL 


of along a constant chord line increases the effectiveness 
throughout the angle-of-attack range® (see Fig. 14). 
This particular arrangement may not be the most 
desirable from a production point of view, but it does 
show what can be done by changes in spoiler location. 

As was pointed out previously, the adverse yawing 
moments of ailerons can be detrimental to the rolling 
effectiveness. Studies of airplane motions*' * have 
shown these effects to be exaggerated by large amounts 
of dihedral. Since sweptback wings, in general, have 
relatively large values of effective dihedral,*’ adverse 
yaw would be expected to be even more serious on 
sweptback wings than on unswept wings of low effec- 
tive dihedral. A comparison of characteristics of a 
flap-type and a spoiler aileron™ (see Fig. 15) shows that 
the yawing moment of the spoiler is favorable over most 
of the lift range and is never so adverse as that of the 
flap-type aileron. Since adverse yaw does, in addition 
to reducing the rolling effectiveness, change the airplane 
heading in the wrong direction, the favorable yaw over 
most of the lift range for the spoiler is an advantage 
to be considered when evaluating different aileron 
systems on swept wings. 

The effect of elastic deformation of a sweptback wing 
on the rolling effectiveness is more serious than the 
corresponding effect on an unswept wing, since the 
wing bending changes the angle-of-attack distribution 
such as to add to the torsional deformation in reducing 
the rolling tendency. To further aggravate the prob- 
lem of elastic deformation, some experimental results 
obtained on a series of untapered wings of the same 
span” (see Fig. 16) indicate increases in twisting mo- 
ment with increased sweep. These results show that the 
twisting moments of the flap-type ailerons increase in 
the order of six times the value at zero sweep when the 
wings are swept 45°. Spoilers of the same span show 
the same trend, but the twisting moments are about 
one-third the magnitude of those shown for the flap- 
type ailerons. Since the data of Fig. 16 are for low 
speeds and the effect of twist on aileron performance 
becomes more serious with increasing speed, it would be 
desirable to see how serious these adverse twist effects 





a 


YAWING MOMENT 


ROLLING MOMENT © Er a 
ee 
“he 
i. 
7, 
“ti 
™ 


a 








O 4 ~ 12 16 20 24 
ANGLE OF ATTACK NACA 


Yawing-moment characteristics of spoiler and flap-type 
ailerons on a 42° swept wing of aspect ratio 4. 
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Twisting moments about the 0.35 chord line for spoiler 
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Fic. 17. 
on an elastic 35° swept wing of aspect ratio 9.4 


would be at high subsonic speeds. To illustrate the 
effects, use has been made of calculations made by the 
U.S. Air Force.** In the Air Force paper, a method of 
estimating rolling effectiveness for elastic wings at high 
speeds is presented, and some sample calculations for 
flap-type ailerons on practical swept wings are given. 
Utilizing the method presented and section character- 
istics of spoilers,** *° some calculations were made for a 
spoiler of the same span as the flap-type aileron and 
having the same effectiveness at Mach Number of 
0.3. The results of these calculations, Fig. 17, show 
the spoilers to be considerably better than the ailerons. 
In fact, the spoilers lose only about 15 per cent of their 
effectiveness as the speed is increased to a Mach 
Number of 0.8, the reversal speed of the flap-type 
aileron. These differences in effectiveness may be 
more pronounced than usual because the wing aspect 
ratio is 9.4, but similar results, although maybe not of 
this magnitude, would be shown for other wings. One 
of the factors contributing to the performance of the 
spoilers is the increase in effectiveness with speed (see 
Fig. 9), while the effectiveness of flap-type controls, in 


general, shows a decrease. 
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Any discussion at this time of hinge moments of 
spoilers on swept wings and their variation with speed 
would be pure conjecture, since no data are available. 
It appears certain, however, that, if a thin circular-are 
spoiler is used, the hinge moments will probably be 
small, since only the edges of the spoiler can produce 
moments about the hinge. If flap-type ailerons are 
used, however, the hinge-moment picture changes, and 
large forces might be expected. 


SPOILER AILERONS USED AS AIR BRAKES 


On many of the high-performance aircraft currently 
in operation or in the design stage, the need of suitable 
devices for use as aerodynamic speed brakes, glide- 
path controls, or both has become well recognized, and 
several applications of such devices to current aircraft 
Speed brakes and _ glide- 
path controls are beneficial for aircraft under various 


have already been made. 


normal or emergency operating conditions, such as a 
rapid descent from high altitude while airplane speed 
is being limited, reducing speed rapidly in level flight 
to increase firing efficiency of fighter aircraft, landing 
on short runways over obstacles, etc. A brief history 
and some of the background in the development of air 
brakes has been traced in a recent article,** and it is not 
the purpose of this discussion to repeat this information. 
However, as shown in reference 34 and the air-brake 
summary data in a recent N.A.C.A. publication,® 
most of the research on air brakes has been on auxiliary 
devices mounted on the wings of airplanes, although 
recently some research and employment of air brakes 
Almost a 
decade ago, during the course of investigating high-lift 
lateral-control configurations at the Langley Aero- 
nautical Laboratory of the N.A.C.A., 
obviate the necessity of including additional devices 


has been on fuselage-mounted devices. 


in order to 
on the airplane, the use as speed brakes of spoiler ai- 
lerons, either alone or in conjunction with slotted flaps, 
was investigated at low speed and was shown to be 
satisfactory.*® More recently, an investigation up to a 
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Mach Number of 0.7 of plug and retractable ailerons 
for use as speed brakes or glide-path controls produced 
similar results,* and other high-speed spoiler-aileroy 
investigations also showed the effectiveness of spoiler 
controls in producing drag.*? In order to illustrate 
the comparative drag-producing effectiveness of several] 
brake devices through a range of Mach Numbers, some 
drag data previously presented in N.A.C.A. publica. 
tions*” * are reproduced in Fig. 18 and show the 
superior drag characteristics, based on brake area, of 
spoiler-type devices as compared with currently used 
shown® 


has been 


that the drag of spoiler brakes increased linearly as the 


fuselage devices. Moreover, it 

a much-desired characteristic, 
In order to illustrate the utility and one of the ad- 

vantages to be gained from the increments of drag 


brakes were extended 


produced by spoiler ailerons used as air brakes on an 
airplane, the computed descent characteristics from an 
altitude of 40,000 ft. of a high-performance propeller 
driven airplane with and without glide-path controls, 
originally presented in a report of the N.A.C.A$ 
The 
airplane had a wing loading of 63 Ibs. per sq.ft., the 


on spoiler air brakes, are reproduced in Fig. 19. 


four engines were assumed to deliver an effective thrust 
of zero with flap retracted, and the glide-path controls 
assumed for the airplane were half-span plug ailerons 
The 
airplane descent was assumed to start at an altitude of 
40,000 ft. and a Mach Number of 0.7, and this Mach 
Number was maintained until an indicated air speed of 
450 m.p.h. was reached. 


projected 8 per cent chord above both wings. 


The latter air speed was 
then maintained for the remainder of the descent to sea 
level. Fig. 19 shows that projection of the plug ailerons 
for use as glide-path controls appreciably decreases the 
time and horizontal distance required to reach sea 
level. This saving of time in descending from high 


altitudes would be particularly important for an 
emergency condition, such as fire or failure of cabin 


pressurization, and also for normal operating conditions, 
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such as for rapid descent from operating altitude at the 
termination of a long-distance flight.*” Supplementing 
their action as glide-path controls, spoiler-aileron brakes 
provide the added advantage of decreasing wing 
bending moments and the possibility of wing structural 
failure by moving the spanwise center of loading in- 
board on the wings.* 

Projection of the ailerons with the flaps deflected in a 
landing approach would also substantially aid the air- 
plane in landing over high obstacles and would ap- 
preciably decrease the length of landing run.’ The 
use of spoiler ailerons as speed brakes to limit or reduce 
airplane speed in a dive or to reduce airplane speed 
rapidly in order to increase firing efficiency of fighter 
aircraft is also feasible, as shown by spoiler-aileron drag 
data and performance calculations for these maneu- 
vers. * 

Some question may arise regarding the adequacy of 
rolling effectiveness of ‘the spoiler ailerons when used 
as glide-path controls or speed brakes. To answer 
that question, the rolling characteristics of spoiler 
ailerons from an initially projected position (as speed 
brakes) have been investigated and were found to 
compare favorably with the effectiveness of the con- 
trols from the fully retracted position.* 

Although a detailed comparison of the characteristics 
of spoiler-aileron brakes and of other brake devices 
beyond the comparison shown in Fig. 1S—is not pre- 
sented herein, several advantages of the aileron brakes 
are readily apparent. These advantages include: 
the variable braking control permitted by the aileron 
brakes as compared with the inflexibility of control 
of some of the other devices; the use of spoiler-aileron 
brakes would eliminate the necessity of including sepa- 
rate braking devices on an airplane; the aileron brakes 
may be used, retracted into the wing, and immediately 
used again, but a parachute brake, for example, can be 
used only once before disposal or repacking (on the 
ground) and is inflexible in control; and, also, the 
spoiler-aileron brake would not adversely affect ad- 
joining wing controls, whereas other devices may. 
In addition, spoilers, when used as air brakes on un- 
swept wings, probably would not cause tail-buffeting, 
inasmuch as the ailerons are placed outboard, near the 
tip, and the wake formed by them would be outboard 


of the tail surfaces. 


DESIGN CONSIDERATIONS 


In the design of a lateral-control system, in order to 
compute the rolling moment produced by a control, 
the effective change in the angle of attack Aa (section 
effectiveness) produced by a given control deflection 
(or projection for spoiler ailerons), as well as the span- 
wise effectiveness, C,/ Aa, of the control, must be known 
because the rolling-moment coefficient is the product 
of C,/Aa and Aa. Inasmuch as spoiler ailerons on 
unswept wings were found to have the same spanwise- 
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effectiveness characteristics as flap-type ailerons,** 
the spanwise effectiveness, C,/ Aa, currently used in the 
design of conventional ailerons*®*~*' may also be used in 
the design of spoiler ailerons. The section effectiveness 
of spoiler ailerons is, however, a complex function of 
aileron projection, wing angle of attack, aileron geom- 
etry, chordwise location of the aileron, and the wing 
profile.’ Since it is impossible to get simple design 
charts of spoiler effectiveness, Aa, similar to those of 
flap-type ailerons,® it is necessary to use data from 
configurations similar to the one under study. The 
lack of data on spoiler controls makes it impossible to 
give even a preliminary design procedure for spoilers 
on swept wings. It is necessary, therefore, to use the 
available data and interpolate or extrapolate the re- 
sults to encompass the proposed configuration. 

In addition to the aerodynamic design, it is immedi- 
ately evident that it will be difficult to get sufficient 
internally hinged circular-are 
Two of the numerous solutions 
The double spoiler retains all 


projection from an 
spoiler on a thin wing. 
are shown in Fig. 20. 
the advantages of spoilers by using two segments that 
are geared together. As shown schematically, the two 
segments stack within the wing and extend in the 
manner shown. This arrangement allows projections 
of about twice the wing thickness, which should, in 
general, be enough. Another solution requires a com- 
promise in the yawing characteristics of spoilers but 
retains the advantages of low wing twist. In this 
arrangement, the spoiler would be able to deflect 
either up or down (Fig. 20). In normal operation, the 
spoiler on one wing would deflect through the upper 
surface and the spoiler on the other wing would deflect 
through the lower surface of the wing. In addition, 
this arrangement has the advantage of allowing the use 
of a simple linkage, since equal up and down deflec- 


tions may be used. 


CONCLUDING REMARKS 


From the briefly reviewed work on unswept wings at 
low speed, it appears that any of several different 
aileron arrangements will provide satisfactory control 
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on airplanes with full-span flaps and that the choice 
for a given application will be determined by the overall 
design of the airplane and the individual preference of 
the designer. 

For high-speed airplanes, however, it appears that 
spoiler ailerons have certain advantages over flap-type 
ailerons in that they give increased effectiveness with 
speed, have low hinge moments throughout the speed 
range, and decrease the aeroelastic problems. 

From the data presented it would appear that there 
are no insurmountable problems associated with 
spoilers on thin sweptback wings. 

It has been shown that spoilers can be used as air 
brakes while simultaneously functioning as ailerons. 
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An Initial Approach to the Overall 
Structural Problems of Swept Wings 
Under Static Loads 


E. E. SECHLER,* M. L. WILLIAMS, ano Y. C. FUNGt 
California Institute of Technology 


ABSTRACT 


For the past 2'/2 years, under the sponsorship of the Air 
Materiel Command of the U.S. Air Force, GALCIT has been in 
vestigating the subject of the stresses in thin swept wings. The 
investigation has utilized theoretical, experimental, and electrical 
agencies on this problem. The problem has been reduced to a 
form amenable to theoretical analysis by assuming the wing re 


placed by a thin cantilever plate of uniform thickness 

A theoretical solution of a thin cantilever plate of rectangular 
cross section with varying angle of sweep has been obtained. Ex 
pansion of the deflection in terms of the normal modes of vibrating 


bars is employed, the coefficients being determined by the 


Rayleigh-Ritz method. Experimental data for both stress and de- 


flection are also presented. An approximate solution yielding the 


qualitative behavior of the stresses in such plates is included 

The equations and circuits for the use of the electrical analog 
type of solution of plate problems are indicated, and a comparison 
is made between theoretical, experimental, and analog data, the 
latter being obtained by the use of a twelve-point electrical net 


work 


INTRODUCTION 


—— IN 1947, THE STRUCTURES GRouP at GALCIT 
became interested in swept-wing design, and, 
through the cooperation and financial assistance of the 
Aircraft Laboratory of the Air Materiel Command, the 
initial effort and thought on this problem has been 
magnified in both theoretical and experimental direc- 
tions. While it was deemed advisable to restrict the in- 
vestigation to static loadings and postpone aeroelastic 
and flutter considerations until a later date, the scope of 
the program was, nevertheless, quite broad and included 
theoretical solutions for idealized wing structures of 
high solidity, with special emphasis upon stress and de- 
flection variation with sweep angle, and experimental 
data for a scale model of a typical wing. It was desired 
to compare these latter data with theoretical results if 
they were obtainable and, if not, to have sufficient data 
to represent adequately the phenomena taking place. 
As an added objective, the effects of taper in plan form 
and thickness were to be investigated. Finally, a study 
of existing design methods being employed in the air- 
craft industry was to be made. 

Presented at the Structures Session, Eighteenth Annual Meet- 
ing, I.A.S., New York, January 23-26, 1950 

* Professor, Guggenheim Aeronautical Laboratory 

t Lecturer, Guggenheim Aeronautical Laboratory. 

t Research Fellow, Guggenheim Aeronautical Laboratory 


It should be stated immediately, however, that the 
answers to all the questions posed in the contemplated 
program have not been obtained, although it is believed 
that satisfactory progress has been made. To date, four 
reports pertaining mostly to the introductory work have 
been submitted, and it is anticipated that several ad- 
ditional ones will be completed in the near future. 


An exploratory study of possible simplifying assump- 
tions necessary to reduce the swept aircraft wing to an 
idealized structure amenable to rigorous mathematical 
treatment, while still retaining the essential characteris- 
tics of the problem, was made by C. B. Ling, many of 
whose observations are included herein. Several modes 
of attack were proposed, but each one was essentially 
reduced to the study of a thin cantilever plate of con- 
stant thickness for a given plan form. For example, 
some of the idealized forms originally suggested were a 
rectangle, parallelogram, triangle, skew parabola or 


ellipse, and sector. 


THE IDEALIZED WING 


While it is not the purpose of this paper to discuss at 
length the reasons for choosing the idealized wing to be a 
thin plate of constant thickness (this matter is covered 
elsewhere’), it is apparent that the assumption of a thin 
wing is justifiable on the basis of the aspect and thick- 
ness ratios used on current high-speed aircraft. On the 
other hand, the chord cross section is not necessarily 
constant or solid, although in the high solidity wings 
being investigated the structural section is approxi- 
mately rectangular. Aside from this, the mathematical 
complications introduced by a variable thickness plate 
were too involved to consider at this stage. Finally, 
thin plate theory is fairly well established, and extensive 
background material exists in the literature. 


The wing has thus been replaced by a thin (cantilever) 
plate of uniform thickness under a normal loading. 
According to standard plate theory, a mathematical 
boundary value problem can be formulated as one in 
which the deflection function w(x, y) is to be found 
from the solution of a fourth-order partial differential 
equation 
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for a given plan form where part of the boundary is 
clamped and part is free. 

Except for certain isolated cases, exact solutions to the 
problem are difficult to obtain. The problems of the 
circular and rectangular plate with simply supported 
boundary conditions are relatively simple and were 
solved early in the nineteenth century. The solution 
for the rectangular plate with clamped edges is consider- 
ably more difficult and was solved during the latter part 
of the same century. When “mixed” boundary condi- 
tions corresponding to some edges free and some edges 
fixed are imposed upon the problem of the rectangular 
plate, the solution becomes much more elusive, al- 
though two authors have arrived at essentially the same 
solution independently' * using energy methods. 


APPROXIMATE SOLUTIONS 


The exact solution of the problem in question is 
generally difficult in the sense that such a solution is 
usually intractable mathematically. Presumably, the 
problem would be made easier if one or more of the 
conditions to be satisfied were relaxed. Suppose, for 
instance, that the condition of satisfying the differential 
equation exactly is relaxed and, instead, an alternate 
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condition-—namely, the condition of minimum energy 
is added. Under such less stringent requirements, 
system of functions can often be found which satisfies 
the so-called “rigid” boundary conditions only.* The 
resulting solution is then expressed as a linear combing. 
tion of the system of these functions, and the arbitrary 
constants associated with the functions are adjusted to 
satisfy the condition of minimum energy. The minimijz. 
ing process presents no special difficulty mathematically, 
This method is the well-known Rayleigh-Ritz energy 
method. 

It may be noted that in this method the boundary 
value problem outlined previously is, in fact, substi. 
tuted by a variational problem. Such a change in the 
nature of a problem usually serves to provide a lead to 
solution originally intractable. 


* The “rigid’’ boundary conditions are those not derivable from 
the corresponding variational integral, such as the imposed specific 
edge deflections or slopes. On the other hand, those derivable 
from the variation of the energy integral are called the ‘‘natural’ 
boundary conditions, to which belong the imposed specific edge 
In the Rayleigh-Ritz method, only the 
rigid boundary conditions must be satisfied by the constituent 
The natural boundary conditions will be satisfied in the 
limit as a result of the variational procedure. 


moments and reactions. 


solutions. 


The energy in a thin plate is given by the following integral: 
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where the double integration extends over the entire region R of the plate and the last two line integrals are ex- 


tended over the free edges of the plate. 


The first integral represents the total strain energy stored in the plate, 


while the others represent the amount of work done by the external loading acting on the plate. 


If the elementary solutions are infinite in number and 
form a complete system in the mathematical sense, the 
solution found by the Rayleigh-Ritz method is exact in 
the limit. 
mate. Such approximation is usually deemed necessary 


Otherwise, the solution obtained is approxi- 


in order to simplify the actual calculation difficulties. 
By properly choosing the set of functions, the results are 
usually substantially correct, even though only the first 
few functions of the system are used. 

If each of such functions belongs to an orthogonal 
system, great simplification will result when minimizing 
the energy integral. In the procedures used in reference 
1, the constituent solutions are chosen as the products of 
two systems of vibration modes of a bar of uniform cross 
section with suitable end conditions. In this manner, 
certain orthogonality relations are utilized to facilitate 
the computations. To be explicit, the deflection surface 
is assumed to be of the form 


w(x, VY) = > Amn Wmn(X, ¥) 


m,n 


where 
Wmn(X, V) = f(x) gn(y) 


In this expression, f,,(%) is the mth free vibration mode 
of a thin elastic bar of length L, clamped at the end 
x = 0, and free at the end x = L; g,(y) is the nth mode 
of a bar of length a with both the ends y = + a/2 free. 
For rectangular plates, x is the spanwise dimension 
measured from the root, and y is the chordwise dimen- 
sion measured from the mid-chord. For swept plates 
(x, y) are replaced by the oblique coordinates (£, 7). 


Upon substituting w(x, y) in this form into the energy 
expression, the first variation of the energy may be 
taken with respect to the coefficients @,. On equating 
this first variation to zero, a doubly infinite number of 
equations in the same number of unknowns are ob- 
tained; these, in general, cannot be solved. An 
approximation is made by arbitrarily taking only the 
first few terms. Assuming rapid convergence of the 


series that was satisfactory for the cases investigated (ol 
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SWEPT WINGS 
the order of ten terms), the a,,,can be obtained and the 
approximate deflection function determined. This 
work has been carried out in detail, including tabulated 
values for necessary calculations and complete proce- 
dural steps for both the swept and unswept rectangular 
plate.’ 

It may happen that the solution is still intractable 
evel after relaxing the differential equation. Alter- 
nately, a solution may be constructed which satisfies the 
differential equation but leaves some of the boundary 
conditions imperfectly satisfied. The latter procedure 
is usually described as the Trefftz method. In this 
procedure, the arbitrary constants involved in the 
solution are adjusted such that the line integral taken 
around the whole boundary—or the square of the error 
arising due to imperfectly satisfying the boundary con- 
ditions becomes a minimum. This implies, in fact, 
that the solution thus obtained is the best solution ob- 
tainable under the circumstances. 

A procedure applying this method to the deflection of 
arectangular cantilever plate under static loading, using 
the elementary solutions of the plate equation in rec- 
tangular coordinates, in the form of products of trigono- 
metric and hyperbolic functions, has been worked out.! 
The actual calculation is laborious and the convergency 
rather uncertain. Another example, which is un- 
published, using the elementary solutions in polar 
coordinates (for plates with a sector plan form), in 
the form of products of a polynominal or logarithmic 
function and a trigonometric function, also indicates 

Thus, the effective- 
though well demon- 


slow unsatisfactory convergence. 
ness of the Trefftz procedure, 
strated for harmonic equations (such as the torsion 
problem), is not readily extended to problems involving 
biharmonic equations of the cantilever plate. In this 
connection, one may reflect that, in the former case, 
many exact solutions are obtained directly as elemen- 
tary solutions, while, in the cantilever plate case, no 
exact solution is known expressible as such simple ele- 
mentary solutions in finite terms. Hence, before some 
special functions particularly suitable to the cantilever 
plate are demonstrated, the relaxation of boundary con- 
ditions may be regarded as impractical except in iso- 
lated cases. 

It is to be observed that in all the preceding solutions, 
exact as well as approximate, the general procedure 
followed is to find the deflection*function first and then 
to derive the stress from the deflection function. As 
the stress is at least a second differential coefficient of 
the deflection function, any error introduced in the de- 
flection function, either inherently or accidently, would 
be magnified in the stress expression. Since it is the 
stress that is of primary concern in practical design, one 
naturally prefers to estimate the stress directly, instead 
of indirectly, from the deflection. 

For this purpose, the compatibility conditions of the 
plate may be relaxed in order to find a system of solu- 
tions which satisfies the equations of equilibrium and the 
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Fic. 1. 


Coordinate system. 


boundary conditions of imposed external tractions. To 
choose among the infinite number of solutions the one 
that solves the problem, the complementary energy 
theorem (or in particular cases with suitable boundary 
conditions, the Castigliano theorem) is applied so as to 
minimize the total strain energy in the plate. The re- 
sult thus obtained will give a direct approximation to 
the stresses. 

Such a direct method, systematically developed by 
introducing some stress functions,' however, is still too 
lengthy to be elegant. A different formulation of the 
plate problem using stress functions suitable for the 
cantilever plate problem, to be solved by the relaxation 
method, is more attractive. This will be mentioned in 
the last section of this paper concerning future work. 
In these procedures, stresses are obtained by a single 
differentiation of the stress functions, instead of the 
second and third derivatives of the deflection function. 


A SIMPLIFIED ANALYSIS OF THE STRESSES IN SWEPT 
PLATES 


A simplified analysis,’ aiming at an approximate pre- 
diction of the stresses in a cantilever swept plate, is 
based on the complementary energy method. A swept 
cantilever plate of fairly large aspect ratio and of the 
practical thickness ratio is usually regarded as a simple 
cantilever beam. In such plates, the stress pattern 
differs appreciably from that given by the simple beam 
analysis only in the neighborhood of the clamped edge. 
To take the effect of sweep into account, the simple 
beam analysis is modified by introducing two parame- 
ters into a system of solutions which satisfies the equa- 
tions of equilibrium but modifies the stress distribution 
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near the clamped edge. The parameters are deter- a set of rectangular coordinates (s, m) be introduced, 


minable by minimizing the strain energy integral. Let with the s-axis coinciding with the center line of the 
plate and the n-axis passing through the inner corner o/ 
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(c) @y CIRCUIT 


FIG, 


cent chord line. It shows the damping of & to the limit- 


ing value of unity as s/a increases. In the root sector 
(ABC in Fig. 1), the stresses decrease quickly as the dis- 
tance from the re-entrant corner is increased. 

These conclusions cannot be regarded as quantita- 
tively accurate. However, they offer a simple picture of 
the stress distribution in a swept plate. 


THE USE oF ELECTRICAL ANALOGIES 


Aside from the application of direct mathematical 
methods for solving plate problems, the use of electric 
analogies has also been investigated by means of the 
Electric Analog Computer in the Analysis Laboratory. 
It has been found possible to devise an electric circuit 
consisting only of passive elements, including trans- 
formers, which is suitable for dynamic and transient 
studies as well as for the calculation of static load de- 
flections. Furthermore, plates of varying thickness and 
irregular plan form present no basic difficulty. 

In this approach, deflections and slopes are analogous 
to voltages, while the shears and bending moments are 


represented by currents. It should be noted that, in- 





(b) ©, CIRCUIT 


DYNAMIC ANALOGY FOR A 
VARIABLE THICKNESS PLATE 
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asmuch as all the physical quantities desired may be 
measured directly, an inherent advantage is apparent 
over the straight mathematical approach in which the 
stresses are obtained by two differentiations with a re- 
sultant loss in accuracy. The network used, applicable 
for variable thickness plates, is given in Fig. +, while the 
representation of clamped, supported, free, and elas- 
tically constrained edges are worked out in reference 5. 
As a matter of record, it may be mentioned that a con 
venient list of finite difference equivalents, useful for 
plate work in general and electrical analogies in par 
ticular, has been prepared by the Defense Research 
Laboratory.‘ 

Because of the limited number of high-quality trans- 
formers available, it was possible to represent only 
With this 
limitation, solutions were obtained for the following 


twelve points of a finite difference network. 


specimens: (1) simply supported rectangular plate, (2) 
rectangular plate clamped on all edges, (3) cantilever 
rectangular plate, and (4) swept cantilever rectangular 
plate. 

The results obtained appear satisfactory considering 
the few number of points available in the network, al- 
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uniformly distributed along the tip. Note: Deflection in inches. 





CHORDWISE DEFLECTION AT TWO SPANWISE POSITIONS FOR VARYING 
SWEEP ANGLES 
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thickness swept back 40°, with a normal shear force of 600 Ibs. 
evenly distributed along the tip. 
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perpendicular to the tip. 
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1-in. 
thickness swept back 40°, with a normal shear force of 600 Ibs. 
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though the work is not being continued because of lack 
of equipment. Additional information in much greater 


detail than included herein is furnished in reference 5, 


EXPERIMENTAL INVESTIGATION 


While it was not part of the original purpose to enter 
into an extensive experimental investigation of the 
stresses and deflections in swept wings or the idealized 
structure, it was thought advisable to obtain enough 
data to check any proposed theoretical solution. Also, 
should particular phases of the analysis prove intracta. 
ble, it was desired to have at hand sufficient experi- 
mental data to admit of a qualitative understanding of 
the plate behavior and perhaps even furnish semj- 
empirical design formulas based, for instance, upon 
beam theory—or, alternately, show the magnitude of 
the beam theory approximation. 

For these reasons, therefore, a fairly careful experi. 
mental program was planned, and subsequently carried 
out, primarily upon a plate specimen of 24 ST, whose 
dimensions outside of the support in the unswept posi- 
tion were 50 by 10 by 1 in. In addition to its unswept 
position, the plate was tested at various aspect ratios 
and several angles of sweep up to 60°. Strain gage 
rosettes were liberally placed over the surface of the 
plate upon which a rectangular grid had been inscribed. 
It was thus possible to record both stress and deflection 
readings for all sweep angles tested. In general, three 
different types of loadings were used: uniformly dis- 
tributed, concentrated tip, and tip torque load. 

Without going into detail regarding the technique or 
apparatus used in the experimental setup, the accuracy 
of the measurements may be summarized with respect 
to the stress and deflection data. For the former, it is 
believed that the overall error is within +5 per cent, 
with approximately half of this figure due to possible 
zero drift of the readings. Accuracy of the deflection 
data cannot be quoted too easily in percentage because 
of the small deflections near the root. In this region, 
care was taken to hold the measurements +0.003 in., 
but further outboard a deviation of +0.005 in. was 
allowed since here the percentage error was small. 

In the report covering this phase of the investigation,* 
the deflection data were presented in the form of con- 
tours. For the most part, the data confirmed earlier 
qualitative ideas regarding the deflection pattern over 
the plate and furnished the qualitative information. A 
typical deflection contour for the case of a 40° swept 
plate under a uniformly distributed shear at the tip is 
given in Fig. 5, while the variation of deflection with 
sweep angle at two spanwise positions is given in Fig. 6 
for the same loading. Similar curves are presented in 
Figs. 7 and 8 for the case of a pure torque applied at the 
tip of the plate. 

In keeping with the desire to present primarily a 
factual report,® it was decided to present the strain gage 
data in terms of principal stress direction and magnitude 
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24 ST aluminum plate specimen of uniform 1-in. thickness swept back 


Fic. 9. Directions and magnitudes of principal stress. 
Solid lines at gage points indicate directions 


40°, with a normal shear force of 600 Ibs. uniformly distributed along the tip. NOTE: 
of principal stress; the numbers indicate magnitude in kips. 
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Fic. 10. Contour of maximum principal stress. 
normal shear force of 600 Ibs. uniformly distributed along the tip. Nove: 
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24 ST aluminum plate specimen of uniform 1-in. thickness swept back 


Fic. 11. Directions and magnitudes of principal stress. 
Solid lines at gage points indicate directions of principal 


40°, with a 45,000-in.lb. torque applied perpendicular to the tip. Nore: 
stress; the numbers indicate magnitude in kips. 
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Fic. 12. Contour of maximum principal stress. 
Note: The numbers indicate magnitude in kips. 


45,000-in.lb. torque applied perpendicular to the tip. 
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at the points where the gages were located. As one ex- 
ample, Fig. 9 shows the data as recorded for the case of a 
40° swept plate under the tip shear load. The stress 
data may also be cross-plotted and presented in the form 
of contours of maximum principal stress magnitude, as 
has been done in Fig. 10. Here, it should be added that 
additional root stress data have been used in order to 
extend the contours to the fixed edge. It is believed 
that contours of this type, coupled with the known 
theoretical knowledge of the direction of the principal 
stress, will prove to be of the most value, since one’s 
intuitive feeling of stress concentrations is borne out by 
the pattern of the contours in the vicinity of the re- 
entrant corner. Similar stress data for the torque load- 
ing are given in Figs. 11 and 12. 

In some respects, the stress data are more interesting 
than the deflection data, although one’s point of view is, 
of course, influenced by the use to which the data will be 
put. From the strength point of view, consideration of 
the stresses will be the more important, and, since there 
are SO many ways one may wish to cross-plot stress 
variations with sweep angle, no attempt has been made 
to do so, leaving the matter to those interested. In 
concluding the remarks upon this experimental work it 
should be briefly noted that the root stress variations 
with sweep, being possibly of the greatest interest, have 
been made the subject of a separate theoretical and 
experimental investigation, results of which are not 
ready for presentation at this time. 


FUTURE WORK 


As was indicated earlier, the work upon this program 
is still under way, and it is appropriate to conclude this 
paper with some remarks relative to work in progress. 
First of all, relaxation methods have been under con- 
sideration for some time, with applications to triangular 
and sector plan forms and plates of variable thickness, 
because of the considerable amount of experimental and 
theoretical work already completed upon the cantilever 
sector plate. Also, some numerical examples for low 
aspect ratio, swept rectangular plates and hollow, thick- 
walled, boxlike sections are contemplated. Then, finally, 
some experimental data are being accumulated for a 
swept box beam with thick cover plates. As work on 
various phases of this project is completed, additional 
reports will be prepared under the same Air Force 
Technical Report number—namely, AF-TR-5761. 

One of the more important parts of this swept-wing 
program pertains to a design method analysis wherein it 
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is hoped that a correlation between practical design ay 
mathematical theory can be achieved. While the styqy 
of plates may lead far afield from the complex Wing 
structures in present-day use, it is nevertheless an aq. 
vantage to treat one example—even a limiting case—jy 
a rigorous fashion and thus provide at least one chee, 
upon the analysis or approximation. In order to carr 
out this design study, GALCIT, with the cooperation oj 
the Air Materiel Command and other Government agey 
cies, has requested aircraft companies and research 
organizations to provide reports and references upon 
the subject. Their assistance in this undertaking jg 
greatly appreciated. 


CONCLUSIONS 


The purpose of this paper has been to outline, jy 
general terms, the problems not yet solved and some oj 
the methods of solution that have been tried. Inasmuch 
as four formal reports have already been submitted 
upon the preliminary phases and are available for the de. 
tailed study of those interested in this problem, it was 
thought more relevant to outline briefly their content 
and indicate how they fit into the broad integrated pro- 
gram than to try and cover the detailed mathematics of 


specific analyses. 
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SUMMARY 
The results of earlier papers are used to calculate the lift and 
moment coefficients for a rectangular airfoil that experiences a 
sudden change in angle of attack (without rotation) in an other- 
wise uniform supersonic stream or which enters a sharp-edged 
The results are valid when the Mach lines from the lead- 


gust 


ing-edge corners do not intersect the opposite side edges. Nu- 
merical results are given in the form of curves. 


(1) INTRODUCTION 


pe INDICIAL ADMITTANCE of a system is defined as 
its response to a step function excitation. In the 
present situation, it is desired to calculate the pressure 
distribution along a spanwise strip x of a rectangular 
airfoil due to an angle of attack along a second spanwise 
strip £, which has the value zero before ¢ = f(£) and 
unity thereafter. This result, which is essentially the 
Green's function for the problem, may then be used to 
find the response of the airfoil to an excitation that has 
an arbitrary distribution in space and time. 


NOTATION 


A.R. = aspect ratio (ratio of span to chord) 
Cr = lift coefficient, Eq. (2.10) 
Cy = moment coefficient, Eq. (2.11) 
J = Bessel function, first kind, zero order 
i! = Mach Number (ratio of free-stream velocity to sonic 
velocity 
l = free-stream velocity 
T = time required to reach steady state, Eq. (4.3) 
a = indicial admittance, Eqs. (2.7), (3.4), and (3.6) 
chordwise velocity distribution, Eq. (2.3) 
k = Green's function for harmonic excitation, Eq. (2.2 
m = (M — 1)/M 
n = (W+1 V 
p = gage pressure 
t = time 
t = time at which disturbance arises at point &, Eq. (2.5 
= (x — £)/U(t — to) 
x = chordwise coordinate (Fig. 1) 
a = downwash angle (Fig. 1) 
8 = M? am 1)} 2 
Y = dimensionless pressure, Eqs. (2.8) and (2.9 
6 = Dirac delta function, Eq. (2.4) 
€ = small, positive, real constant, Eq. (2.5b) 
6 = Mach angle (sin~!1/M, Fig. 1) 
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g = running value of x 

p = air mass density 

T = dimensionless time, Eq. (4.4) 

w = angular frequency 

The indicial admittance for the two-dimensional 


supersonic airfoil has been obtained by several writers! ~ 
through a variety of approaches. The method to be 
used in the present analysis is that of reference 4, 
wherein the indicial admittance was determined from 
the transfer function—which is to say that the response 
to a step function excitation was obtained from the 
response to a harmonic excitation. This method is 
chosen principally because of the availability of the 
latter result® but also because of the pedagogic emphasis 
that it places on the equivalence of the two results, 
either of which may be used to analyze arbitrary dis- 
tributions of the excitation in time, by virtue of the 
Duhamel and Fourier methods of superposition, re- 
spectively, for linearized systems. 


(2) STATEMENT OF PROBLEM 


The chordwise (gage) pressure distribution—i.e., the 
force per unit chord acting on a spanwise strip at x—on 


a rectangular airfoil of aspect ratio A.R., leading edge 
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x = 0, and trailing edge x = 1 in a supersonic flow The function a[x — & ¢t — &(&)] is defined as the 
directed along the positive x axis with velocity U and  indicial admittance of the airfoil and gives the (dime. 
Mach angle 6 (Mach Number, M = cse @),asshownin _ sionless) pressure at station x due to a downwash incre 
Fig. 1, due to a downwash angle distribution ment at station £, which vanishes for time earlier than 
; i,(&) and has the magnitude f(£) for time later than 4(¢ 
a(é, t) = a(&) exp (tot) os ’ ‘ : ; s 
rhe first term, ado, gives the two-dimensional response 
is given by® and the second term, a, gives the defect in the three. 
; dimensional response, relative to its two-dimensiona| 
12 x P 
pl* : . siliiiee diiaias thee, Siteliaetat ailit ‘Wain ie te 
ie a = ( yf k(x — t)a(é, Adt (2.1) value, due to induction effects. Since the two-dimep. 
B 0 sional problem has been solved previously,!— ; 
suffices herein to consider the three-dimensional effects 


R(x) = Ro(x) — (BA.R.)—1Ry(x) (2.2a) as described by a. 
a ry In the remainder of this paper it will be convenient to 
Ro(x) = 6(x — &) + (5 a *) x make the pressure dimensionless by writing 
x 


|exp Cae) Tp (a) (2.2b) p(x, t) = —(pU?/B)y(x, 2b) (2.8 
Bl B7L 1 | 


y(x, t) = yo(x, t) — (BA.R.)—! yi(x, 2) (2.9a 


ra) tw ~% | Fi 
aati 4 - r) 7 vo, a(x, t) = | do, [x — & ¢ — bo(E))f(OdE (2.9 
| (5) (“37 *) ; (SF) 
exp = sin =e (2.20) a 
wM Bl Bl The problem to be solved may now be posed as: pl 
8 = (Mt 1) (2.3) Given f(&), find y(x, t); this pe has been reduced be 
to two successive quadratures, as indicated by Eqs, 
lim f "be — Qdt = 1 ww ero ieee ay 
comet Jy —e While some applications may require a direct knowl- 
where Jo is Bessel’s function of the first kind and order edge of the ee pee the airfoil, the lift and saint , 
zero and &(x — &) is the Dirac delta function. The ment coefficients are generally the items of principal 
unit of length is taken as the chord of the airfoil, the interest. The increments in these coefficients are given 
notation having been slightly modified from that of ad 
references 4 and 6. The result is linearized and holds 7] 
only when the Mach lines from the leading-edge cor- —9C:(t) = (pl*/2)~! / [—2pi(x, t)|dx (2.10a ' 
ners do not intersect on the wing. The latter restric- a ti 
tion is met by demanding the inequality BA.R. > 2 = (4/6°A.R.) J v1 (x, t) dx (2.10b) | 
although it can be shown® that the final results for the , 
lift and moment coefficients are valid provided only that aa ” ‘ 
BA.R. > 1, which implies that the aforementioned —6Cu(t) = (pU*/2) | [—2pi(x, #)] (1/2 — x)dx 
Mach lines do not intersect the opposite wing edges. (2.11a 5 
The downwash distribution to be considered is given Ps. | ” 
by = 59 = (ax) f xyil(x, t) dx (2.11b 
a(& #) = CE) Tle — Wf) (26a) The moment coefficient is taken about the mid-chord, 
s f(é) lim Z. luke +- di)~Y exp lewlt — &)I and the reference moment arm is the chord (= 1). 
_ “eset ame tw)~* exp [tw 0 
(2.5b) (3) INDICIAL ADMITTANCE 
1 (¢ — t) is the Heaviside step function, and Eq. Consider first the integration of Eq. (2.7b), assuming 
(2.5b) introduces its Fourier representation in the w that reversing the order of the integrations with respect 
spectrum. This distribution is introduced in Eq. to ~ and w is justifiable, which may be verified a poster- 
(2.1), and the pressure may be expressed in the form ori. Carrying out the indicated differentiation 1 | 7 


Eq. (2.2c) and substituting in Eq. (2.7b) yields 


rs ' t 
p(x, t) = — (’ yf a[lx-—ét— ty(é) | f(e)dé (2.6) : Q 
B a,(x, ft) = (1/2r) lim [ (e + tw)! X 


‘ 
a(x, t) = do(x, t) — (BA.R.)—'a,(x, ft) (2.7a) : Mx t 
| _ . _— | exp E (: — Fer | x t 
do 1(x, t) = on lim dw (€ + iw) ~' ky 1(x) exp (iwf) me ( ; ) es} . 
). . cos as sin a dw (3.1) | 
(2.7b) Brl M Bl 
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Fic. 2. Two-dimensional: lift coefficient for sudden 


change in angle of attack. 


If the trigonometric functions in Eq. (3.1) are re- 
placed by their exponential equivalents, the result may 
lim 


T e— 0+ —« 
" exp i (: _ = )| + n exp ie (: — = )| dw 


(3.2) 


be written 


a(x, t) = 


(3.3a) 
(3.3b) 


m = (VM — 1)/M 
n= (M+ 1)/M 


The integrals may be directly interpreted as step func- 


tions [ef., Eq. (2.5)|, whence 


m , x n ; x 
a(x, t) = l (uv _ ) + 1 (ui — ) (3.4) 
2 n 2 m 


It is evidently convenient to introduce the single 
variable (x/ Ut) and restate Eq. (2.9) in the form 


yi(x, t) = ft ay Se L reyas (3.5) 
aii 0 LU[t — w(e) f°” 

a(w) =~ 0, w<O (3.6a) 
=i. 8< 8s mw (3.6b) 
=m/2, miawein (3.6¢c) 
= wes (3.6d) 


(4) SuppEN CHANGE IN ANGLE OF ATTACK 


The simplest transient problem of practical interest 
is the sudden change of angle of attack (without rota- 
tion) over the entire wing. [Physically this might 
arise from a suddenly imposed, vertical (downward) 
velocity, but the primary value of the solution is as a 
basis for Duhamel integration.] Consider, therefore, 
the downwash distribution [ef. Eq. (2.5)], 


f(é) =1 
to(&) = 0 


(4.la) 
(4.1b) 


SUPERSONIC 


















































RECTANGULAR AIRFOILS 649 
.06 
> 
O 
O5 |ralst r™ M= 1.2 
o4 - 
osl_{ | 
.02 ie 
ol J 
O o & tS 8 1.0 
Fic. 3. Two-dimensional moment coefficient for sudden (unit) 


change in angle of attack 











LO 
> M=4 
8 4 h- 
ay U2 


| |G 
































4 

P 

oO a 4 T 6 8 1.0 
Fic. 4. Increment of lift coefficient due to induction effects for 


sudden (unit) change of angle of attack 
for which Eq. (3.4) reduces to 
y(x, tl) = 


It is evident that the steady-state pressure will exist 
over the entire airfoil (0 < x < 1) after a time 7 given 
by 


T = M X (chord)/(M — 1)U (4.3) 


and it is convenient to introduce the dimensionless time 
parameter 


(t/T) =m Ut = ((M—1)/M]Ut_ (4.4) 


T= 





By substituting the indicial admittance from Eq. (3.5), Eq. (4.2) yields | 
m(x,t) =x, *S 7 (45a 
m n 
= >*¥ +o, Tix (n/m)r (4.5b 
= nt, x > (n/m)r (4.5¢ 


The lift and moment coefficients due to the pressure distribution [Eq. (4.5)] are given by Eqs. (2.10) and (2.1; 
Adding to these results the two-dimensional coefficients of reference 4 (allowing for the differences in notation), ; 


is found that 


Px n m 
Cy, = (mn)” — (2BA.R.)—! | 2n7 — ri. rs (4.6a 
4 m n - 
1 (mn)? m l a / 
= -—cos-! [M(1 — mr)] + cos~!| M{1— + — (2nr — n*r? — mn)” — (2BA.R.)—! x 
T T 7 Tv 
m nm. m ; 
= or — 2h {75.1 (he 
= ] — (28A.R.)—', De (4.6¢ 


p .. m~'n'*(1 — m)? n\ , nN es m 
Cy = tr? + (12BA.R.)—' | 3 rT? 2 -) (2 — 2m + m*) * Il, 7 < (4.74 
| m m* n 


m”~“n’*(1 —m)? . m l T ie ; 
= r*>cos~'| M1 — Sa (1 — (2nr — n*r? — mn)” + 
T do m 


tor 
; m on F m £ ace 
(12 BA.R.)— (" + ge lina nr), <7rS1 (4% 


(4.7¢) 


a 


= (128A.R.)—', +> 1 


The two-dimensional results C,, and Cy,,, as taker are 





from reference 4, are plotted in Figs. 2 and 3. The 
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sharp-edged gust. Eq. (3.5). 
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Cwo-dimensional lift coefficient for uniform (unit) gust 
entry 
accordance with the relations 


(4.Sa) 
(4.Sb) 


Cr = Cr, + 6C; 
Cy = Cu + 6Cy 


are plotted in Figs. 4 and 5. 


(5) Gust LOADING 


A second special case of practical importance is the 
entry of the airfoil into a sharp-edged gust. In this 
case, the downwash distribution, which is depicted in 


Fig. 6, is given by Eq. (2.5), where 


(5. la) 


(5. 1b) 


f(g) = 1 


Ut(é) =E 


and Eq. (3.5) reduces to 


* — 
v(x, t) = / a ( : ) ae ($2) 
Jo "\Ut —f; 


In carrying out the integrations in Eq. (5.2), it is ex- 
pedient to examine the behavior of w in the (Ut, &) 


Adding the results of substituting this result [Eq. (5.4)] in Eqs. 
and moment coefficients of reference 4* yields 
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Two-dimensional moment coefficient for uniform (unit) 
gust entry 
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plane—viz., 


w= (x — £)/(Ut — &) (3 


Ut = [x + (w — 1)é]/w (5.3b) 


oa) 


Eq. (5.3b) one-parametric family of 
straight lines, the lines of dominant importance being 


shown in Fig. 7, together with the domains of a)(w) 


represents a 


which they delineate. 


Carrying out the integration indicated by Eq. (9.2) 
with the aid of Fig. 6 yields 
v(x, #4) =x, xSf (5.4a) 
l m l n ; n 
= s+ ;, rs 2s T 
2\m — 1 Z2\n — 1 m 
(5.4b) 
=), xs & (s/n)r (5.4¢) 


(2.10) and (2.11) to the two-dimensional lift 


(5.5a) 


Wk m m n . 
= -cos~! [M(1 — nr)] + ( r cos-'| M{1 — — (2BA.R.) uM (- + nr — - rT’ 
1 aw \m T 2 2 


= 1 — (26A.R.)-', +r >] 


m/n<71r< 1 (5.5b) 


(5.5c) 


* > (< ) \ ° o,* e . a* . . . 
Eq. (28b) of reference 4 should have contained an additional factor of 1/x in the last term, for which reason Fig. 6 therein is 


inc rect 
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= (128A.R.)—', r>1 
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The results of Eqs. (5.5) and (5.6) are plotted in Figs. 8-11, in accordance with Eqs. (4.8). 
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The Gyroscopic Effect of a Rigid Rotating 
Propeller on Engine and Wing Vibr ation 


M odes 
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TRUMANt 
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SUMMARY 


In many wing vibration analyses it is found necessary to take 
into account the effect of Hence, 
it is reasonable to ask what vibratory gyroscopic effect this flexi 
An engine 


flexibly mounted engines 


bility may give rise to when propellers are whirling 


mount may be thought of as a horizontal beam cantilevered 


wing, 
were 


vertical flexibility 
then, 


from the having both horizontal and 


If this beam 
vibrated, the gyroscopic moments induced in the propeller due 


infinitely rigid horizontally, when it 
pitching motion of its axis would not produce 
lateral 
actually be 


to the resultant 


propeller axis yaw However, engine-mount stiffness 


which allows propeller axis yaw) may less than ver 


tical stiffness, so that gyroscopic effects will play a role as the 
propeller axis undergoes pitching vibrations at the tip of the 
The purpose of this paper is to in 


assumption that the propeller itself 


cantilever engine mount 
vestigate this role under the 
isa rigid disc 

rhe 
with classical gyroscope theor y 


Part (1) deals briefly 


Part (II) presents engine vibra 


is divided into four parts 


paper 


tion mode studies —experimental photographic techniques on a 
model gyroscope mounted at the ends of two different cantilever 
beams. Part (II1) presents the theory of the coupled motion of 
an elastic wing upon which a gyroscope is mounted to simulate 
Part (1V) consists 


in which, 


an engine-propeller system on an airplane 


of an example of the theory of Part (II1), by taking 


what thought to be reasonable parameters, results are ob 


tained showing how the whirling of a rigid propeller may mate 


are 
rially affect wing normal mode shapes and frequencies 


(I) Baste GyRoscoPpE THEORY 


Pp R A GYROSCOPE (Fig. 1) with one point 
at the origin 0 of a rectangular coordinate system 
0-X-Y-Z fixed in space. Let 0 also be the origin of a 
rectangular coordinate system (0-1-2-3 fixed in the gyro- 
scope. Assume that the gyroscope has mass moments 
of inertia J), Js, 
and 3 and that it is symmetric about axis 3 
h=1, It that the of 

then the 3-axis and that the gyroscope is spinning about 


respectively, about its axes 1, 2, 
that is, 
figure” 


and 7/s, 


is assumed “axis is 


this axis. Under these conditions, the position of the 
gyroscope at any time ¢ is uniquely defined by the three 
Euler angles — 6, ¢, and y. 

It is of interest to the present discussion to conceive 
at 0 with a built-in spring sys 
restraint 


of a “universal joint”’ 
tem having the sole effect of providing a 


Presented at the Aeroelasticity Session, Eighteenth Annual 
Meeting, I.A.S.. New York, 

* Associate Professor of Aeronautical Engineering 
Now, 


General 


January 23-26, 1950 


Development Engi 
Com 


T Instructor, Physics De partment 


ueer, Thermal Power Systems Division, Electric 


pany 


653 


against nutation. Let the spring constant there be A 


so that the restoring torque becomes 
M, = —Ké (1) 


If no other torques act on the gyroscope, the equations 
of motion become (in the absence of gravity) 


16+ Ke = lhe¢ sin 8 cos @ — Ing sin 6 (¥ + ¢ > cos 6) 
pe, = const. = J; sin? 0g + TI; cos 6 ( (¥ + ¢ cos 8) 7 (2) 
py, = const. = 13(~ + ¢ cos 8) 


At time ¢ = 0, let the following representative initial 


conditions hold: 
6=6+ 0 
4 = 0 » 
er fo = 
y=y=2+0 


Then, initially, also 


Info = _ [;Q = po 


Po cos A 


ay 0 
Pm - 


The solutions of Eqs. (2) for @ and g, under the as 
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sumption of small nutation angle 6, are 


| 
0 = 0 : OS wl 
whe iar 


si — poo ( l NC sin “) 4 a 
Ty sin 1+ P Ww ~— 


where P bo?/ KI, is nondimensional and w? 
(K/1,)(1 + P). Thus, the small motions of nutation are 
sinusoidal about some @-value, and the value of w? indi 
cates that nutation frequency increases with increasing 


(+4) 
(5) 


gyroscope spin. 

Two limiting cases are of interest: first, when P > | 
and, second, when P < 1—4.e., where py > V ‘KI, or 
po< VKI,, respectively. In the first case, nutation 
is almost nonexistent, since 6 maintains a near-constant 
value %. In the second case, the gyroscope nutates 
about a 6-value near zero with appreciable amplitude. 

Therefore, with nutation small, the gyroscope pre- 
cesses in a quasi-steady fashion, the variation about a 
steady state being sinusoidal. In the limiting cases 
above noted, when P is large, the rate of precession is 
diminished, while the variation from a steady value 
becomes faster and of lesser amplitude. At extremely 
large values of P, there is effectively no precession. 
On the other hand, when P is small, precession is more 
rapid and has oscillatory motion of large amplitude 
superimposed. 

The equations of motion given above may be replaced 
by equations in another useful set of coordinates. This 
proves to be fruitful for two purposes: The first is the 
discussion of the motion of, say, the c.g. of the gyro- 
scope as seen by an observer looking down along the 
i.e., the projection of the motion of the c.g. 
the second is the discussion of the 


Z-axis 
on the x-y plane; 
gyroscope motion when the previous “universal” spring 
constant in the nutation direction is replaced by two 
independent spring constants acting in the mutually 
perpendicular fixed planes z-y and 2-x, respectively. 
Such spring constants will be alluded to in the practical 
examples to follow. 

In terms of the approximate “‘projection coordinates, ”’ 


x= KOsmeo 


(6) 
— R86 cos of 


y = 
the equations of motion of the gyroscope can be shown 


to be 


NX + py + Ks = O 


; ‘ 7) 
hv — pyx + Ky = Of M4 


where now the universal nutation spring of stiffness 
K has been replaced by the two stiffnesses A, and AK, 
in perpendicular planes (A, and A, are measured, for 
example, in inch pounds per radian). 

The general solutions of Eqs. (7) are each sums of 


t tavl 
lies “and may be 


four terms having factors like e€~", « 
expressed in the form 
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x = A cos (ul — «) + B cos (vt — &) 


with a similar equation for y, where A, B, €, © are ap 
propriate constants. 

The entire discussion to this point follows classieaj 
lines'~* and is included to serve as introduction to sub- 
sequent parts. 


(11) APPLICATIONS OF CLASSICAL GYROSCOPE THRopy 


The questions that originally motivated the present 
study were to determine the dynamic effect of spin 
on an elastically mounted rotating body and, in par 
ticular, the effect of a spinning propeller on the vibra. 
tion modes of an aircraft engine mounted on a cant; 
lever engine mount. 

Bennett* has examined this problem partially for g 
typical small, high-speed, shock-mounted motor, such 
as an aircraft accessory. However, it is felt that the 
present discussion improves and extends his treatment 
He does not discuss the vibratory modes, which are of 
greatest interest to the problem; nor does he evaluat. 
properly the role of the spring supports, which are not 
generally oriented simply in the nutation direction as 
he assumes. Finally, intermediate steps in his numeri- 
cal evaluation of the increase in nutational frequency 
due to spin appear to be in error.* 

For the discussion of the case of the aircraft engine, 
the z-axis, Fig. 1, will be assumed horizontal pointing 
forward in the of The 
x-axis will be taken positive vertically downward, while 


direction propeller thrust. 
the y-axis is horizontal positive to the left of an ob- 
server looking in the direction of positive z. 

The cantilever engine mount may be assumed effec- 
tively as a simple cantilever beam. The aircraft or 
aircraft wing supporting the system will be considered 
completely rigid and stationary in space, as for example 
a rigid test stand. It is readily shown that a single 
force or moment acting at the free end of a cantilever 
beam causes the beam to slope so that the tangent line 
to the beam at its free end passes through an effective 
“rotation point’ about which the end of the bent beam 
appears to pivot for varying load. Associated with 
this rotation point is, then, an effective radius R from 
the pivot to the beam tip. When both forces and 
moments act simultaneously on the beam and vary 
with time as well, they must bear a constant propor- 
tionality to each other in order to maintain, at all times, 
the same effective rotation point for the beam tip. 
In harmonic oscillations, this constant proportionality 
is satisfied, and it will be assumed approximated in the 
practical cases to be discussed here. The validity of 
this approximation is akin to that of the usual assump- 


* Bennett appears to give an incorrect value C/A = 13/li = 
(0.019 for the rotor in his example. The present authors find this 
value to be 0.0085, which, taking K = 0 in w? = (K/I,) + (bo/I 


yields f = w/2r = 0.85 cycles per sec., a value comparable to 


Bennett’s 0.95 cycles per sec 
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tions of effective elastic axes of wings or empennages in 
flutter analyses, etc. 

Thus, the propeller and engine parts that rotate will 
assumed to be an effective rigid rotor; this rotor, 


be 
ie for nutational motion about the effective rota 
he) 


yvotil 
a point described above, then forms a classical sym- 
metric gyroscope, with point 0 of Fig. | now at the 
pivot point. The effect of gravity on the horizontal 
gyroscope call be ruled out by considering the equilib- 
rium position to be the actual steady position of the 
system with rotor spinning and gravity acting. Under 
these assumptions, the effect of the spinning propeller 
on the vibrational characteristics of the system may be 
calculated, employing the theory of Part (1). 

A model was built to study qualitatively what the 
gyroscopic effect on engine vibration mode shapes is 
likely to be (Figs. 2 and 3). The gyroscope spin axis 
was hung vertically dgwnward, rather than horizontal 
as for an engine-propeller system, thus providing the 
opportunity to employ relatively flexible supports in the 
role of gyroscope mount without unwanted distortion 
due to gravity. This allowed correspondingly large 
gyroscopic motions so that the modal shapes could be 
clearly ascertained. 

Under the model and attached to it, as shown in 
Figs. 2 and 3, was mounted a strontium are lamp 
that provided an effective point source of light. A 
portrait camera mounted vertically beneath the lamp 
photographed traces of the motion over 
specified time intervals, the camera shutter being left 
Figs. 


gyroscope 


open in darkness during the entire interval. 
!-10 show typical traces obtained in this manner. As 


these traces actually lie in the x,y or “projection” 
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coordinates, they illustrate motion of the type repre- 
sented by Eqs. (7). 

Fig. 4 trace of the motion of the gyroscope 
mounted on a beam of square cross section and given an 
initial displacement (nutation) with no initial preces 


Figs. 5 and 6 show the corresponding motion of a 


is a 


sion. 
beam of rectangular section, not square, initially dis 
placed in the direction of lesser and of greater stiffness, 
respectively. Fig. 7 indicates motion under more gen 
eral initial conditions, the gyroscope having an initial 


displacement together with an initial precessional ve 
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locity. Figs. 4, 6, and 7 covered a period of 1.5 sec., 
while Fig. 5 was for 3.0 sec. 

Figs. 8, 9, and 10 show damping of the oscillations, 
resulting in a decrease of amplitude as the characteristic 
figure precesses. This damping has not been con 
sidered in the theory presented herein but may be taken 
into account by the introduction, for example, of a dis- 
sipation function into the kinetic potential. 


trace of the same motion as that of Fig. 4 for a period of 


Fig. Sis a 


15 sec., while Figs. 9 and 10 correspond to Figs. 5 and 
6, respectively, for a period of 10 sec. The small os 
cillations observed on all experimental traces result 
from vibrations in the gyroscope at its spin frequency. 
It should be realized that the traces presented herein 
are by no means general for all systems but depend for 
their shape on the various parameters of the particular 


OCTOBER, 


L950 


system—namely, J, J;, 2, A,, and A,—as well as the 
It was found possjbj, 
to predict qualitatively the motion indicated in Fig. 4 


No the 


initial conditions of the motion. 


for example, using the theory of Part (1). 


retical obstacle would prevent the point-for-poin 


verification of Figs. 4-10 if this were desired, althoug 


the process is laborious. 


(III) THrory oF THE GyROSCOPE WITH ELAsticayy 
MOUNTED PIvorT 


The gyroscopic theory and experiments of Parts (| 


and (II) provide a firm basis for undertaking a furthe; 


problem—that of the gyroscopic beam system mounte 


on an elastic, rather than a rigid, base. The system 0 
Parts (I) and (II) simulates a propeller-engine-moun 
system attached to a rigid jig or foundation. Her 
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ENGINE AND WING 
this foundation becomes elastic, playing the role of the 
aircraft wing. In particular, the gyroscope will be 
taken in the orientation of a rigid propeller with axis 
horizontal as described in Part (II). However, the 
pivot point (origin of x, y, and z coordinates) will now be 
assumed to move vertically, with deflection /, as part 
of a vibrating aircraft wing (Fig. 11). 

If, as in many modern vibration analyses,’ the wing 
is approximated by a series of N individual masses 
m,, then the kinetic energy of the wing with elastically 
mounted engine, but with propeller not rotating, is 

N 
Tw = (1/2) > mAh? (8) 
k=1 

Although a greater number of normal modes may 
enter, for simplification it will be assumed here that 
only two normal modes of the wing are brought into 
action by the presence of the spinning propeller. Thus, 


hy = Iyadi + Ayabe (9) 


where h,, and hy. are the nondimensional values of the 
displacement of mass m, in the first and second normal 
modes, respectively, while £, £ are the associated gen 
eralized coordinates (with dimensions of length). 


Therefore, 


N 
Tw = (1/2) > my(ya2és? + Aye*&?) = 
k=1 


(1/2)(AiE,? + Ask?) (10) 


since orthogonality between normal modes implies that 


N 
+» mAh = 0 (11) 
I 


Then the potential energy of this motion is 


(1/2)4?A 181? + (1/2)we*A oko? (12) 


Uw = 


where w; and w, are the circular frequencies of the first 
two normal modes. 

Let the linear horizontal engine mount spring con- 
stant be k, = K,/R?, so that energy stored by y deflec- 
tion is (1/2)k,y*. Let h, be the vertical displacement 
in space of the engine-propeller mass and h, be the verti- 
cal displacement in space of the pivot point of the gyro- 
scope system. Under these assumptions, the kinetic 
potential of the entire system is 


$ RSs ng v2) 4 E — yx 
= x* y Ds 
aR? *~ TY) T Pool ve 4 ye 
i ~. ale ee 
( ~ 2R? ) eg Te oe 
where 


x= h, = ha = ak; + Ako (14) 


Thus, @ is a function of &, £, and y. 
Lagrange’s equations in §, &, and y yields three equa- 
tions of motion, the coefficients of which are constants 


Application of 
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of the system. The equations of motion are thus ex- 


pressed as 
aye, + dye + ay + dk = r, 
Gok, + Dobe + co + dot. = (15) 
eV + fe = i gto + ky a a 


where 
= + A, “on ae | 
— hdyaz — 1,a,* ce. 
R? " ke 

apy ay YW, 

ae aR - (16) 

d, = "Aj, dy = we"A» 

e = 1,/R? 

f = -—pya@/R 

g = —p,, a/R 

k =k, 


It may be noted that the case of zero gyroscope effect 
present can be obtained by theoretically moving the 
propeller off to infinity in the z-direction—i.e., R— o. 
The equations then reduce to those of two noncoupled 
normal modes, and the value of horizontal engine pro- 
peller displacement y is found to be zero, as expected. 
The terms with denominator R? are thus the gyroscopic 
coupling terms in this instance. It is seen that their 
effect is generally increased if any of the parameters in 
them cause the terms themselves to increase. In par- 
ticular, large engine-propeller pitching inertia (J;) 
about the rotation point, or large propeller rotational 
inertia (J), and small distance to effective pivot point 
(R) will tend to increase the gyroscopic coupling effect 


of the rotating propeller. The solutions 
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TABLE lb 


Gyroscopically Modi- 
fied Modes of XR60-| 


TABLE la 
Normal Modes of XR60-1 Wing 





Wing 
fi: = 122.8 Cycles per Min fi = 122.0 Cycles per 
Min ; 
fo = 144.1 Cycles per Min fy = 134.3 Cycles pe 
Min: ; 
Mass hy hey hy h 
1 4. 6664 —0.7938 3.2270 —1.93 
2 4.6757 —0.4114 3.3346 —() | 
3 2.4974 —().6292 1.6731 —(). 87 
4 2.5740 —(0. 1365 1.8594 —() 352: 
5 1.3136 —0.3741 0. 8687 —().5 
6 1.3826 —(.06373 1.0013 —() 
ij 0.4604 —().3758 0.2400 —() 
8 0.6701 +0.4345 0). 6080 +() 
y 1.0000 +1, 0000 1. 0000 + | 
10 —(), 1423 —(0. 2562 —(.1723 —(). 2654 
11 —0.1167 +0.01035 —(0. 08321 +0) 
12 —(0. 4662 —0 2304 —(),. 4041 —().2] 
13 —(0.4872 +(). 09644 —(). 3334 +() 
14 —0.6962 +(0).4338 —0. 3983 +0). 5245 
15 — 0.6379 —(0. 09191 —(). 4940 —() 
TABLE 2 
Inertia Properties of Idealized Wing 
(Masses in Slugs) 
m = 13.20 mg = 45.45 
Mm. = 13.20 mg = 210.8 
m; = 20.05 my = 228.1 
ms, = 20.05 my, = 228.1 
ms; = 214.6 Me = 88.9 
me = 214.6 m3; = 88.9 
m, = 45.45 my = 210.8 


ms = 1215.3 


may be assumed for these equations, from which are 
obtained the characteristic values \—1.e., solutions of 


the sextic in A 


a," + dy bd? OA 
Ay? bod? + do Cod =@Q (Is 
fr gn ed +k 


This sextic proves to be a cubic in \*, with all roots 


negative. Hence, the solutions are of the form 

£) = &1e a (19 
the motion 
The three 


etc., where @ is a characteristic frequency; 
is thus seen to be stable simple harmonic. 
frequencies @, a, and a; represent the frequencies of 
the gyroscope-coupled system. 


VIBRATION OF A WING WITH A SINGLE 
ROTATING PROPELLER 


(IV) EXAMPLE 


Reference 5 presents a convenient example for study 
in the XR60-1 wing, normal vibration modes of which 
have been calculated assuming flexibly mounted en- 
gines. 

The plan form of the masses as given in reference 5 
is shown in Fig. 12. Table la gives the first two nor- 
mal modes as calculated in reference 5, both modes 
being normalized on mass my. Table 2 gives the values 
of all masses in slugs. 
it will be assumed that only 4 
although there are tw 


For simplification, 


single propeller is rotating, 
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ENGINE AND WING 
engines and propellers on the wing. The outboard 
propeller will be assumed as rotating; this is located 
at my. It will be further assumed that m7 is the location 
of the pivot point for the gyroscope (rigid propeller). 
Thus, R = 3.28 ft. The value of x is given by x = 
hy — hs or 


xX = Abi + Anke (20) 


where 


hoy hs, A = ly — Igo 


a = 
The example to be worked will contain the following 
parameters ; 


0.3299 


a . 

ds = 0.5655 

L, = 200 slugs-ft.? 
I; = 235 slugs-ft.? 


The determinant [Eq. (18)] becomes 


2,502.102A7 + 413,433) 


—1,207.411A 


The solutions for \ are 


Ay = 12.78162 = 1a) 
As = 75.60102 = 103 


and, therefore, the characteristic frequencies are 


f; = 122.05 cycles per min. 
fs = 134.33 cycles per min. 
f; = 721.93 cycles per min. 


The first two of these may be compared to given normal 

mode frequencies of 122.8 and 144.1 cycles per min. 
Upon entering \, into the equations of motion [Eqs. 

(15)], the first modified mode is obtained, wherein it is 


found that 


£9 = +0.358£10 


whence the new mode is represented by 
hy = (hia + 0.358hy2)é; 


Hence, the new first mode contains about 36 per cent 
of the original second. Therefore, in the modified 
second mode, which is obtained by solving the equations 
of motion with \» known, there is obtained 


hy. = [hya(—0.074) + hy] & 


Thus, the new second mode contains 7 per cent of the 
original first mode. The new modes as modified are 
given in Table 1b. 

The frequency shift is seen to be negligible in the first 
mode and approximately 7'/2 per cent in the second. 


3.468172 
3.46817\? (352.179A? + 78,841.3) 
— 2069.69 
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R = 3.28 ft. 
h = 91,875 lbs. per ft. (taken as 75 
per cent of vertical stiffness) 
A, = 2,500 .079 
A» = 346.254 
Po = 39,375 (propeller speed is 1,600 
r.p.m.) 
fi = (60/27) o 122.Scycles per min. 
fy = (60/2m)e. = 144.1 cycles per min. 


In assuming propeller inertias, the values of /;, and 
especially of /; (which are based on simple geometric 
estimates), are probably much less than those actually 
to be The thus 
thought to be on the conservative side. 
R cannot be known with accuracy except when more 


encountered. calculated effect is 


The value of 


detailed data are given. 


1207.411X 
2069.69A = 0 
(18.59012A? + 91,875) | 


This is of the order of magnitude encountered in many 
cases wherein engine vertical flexibility alone is deemed 
worthy of inclusion in the vibration calculation. 

The magnitudes of the effects shown in the present 
example must not be taken too seriously, since the as 
sumptions involved are strong. The point of the situ- 
ation appears to be, however, that it is not unreason- 
able to expect propeller gyroscopic effects, say in flight, 
to modify the results of either calculations or ground 
vibration tests. Whether the modifications will be 
large or small depends on the particular system param- 
eters involved. The present analysis suggests that the 
modifications may, on occasion, be large if some of the 
(a) large values of 
propeller-engine pitching inertia; (b) small effective 
radius to ‘“‘pivot point’’ of gyroscopic system; (c) large 
propeller rotational inertia; (d) more than a single en 


following factors are encountered: 


gine running; and (e) low engine-mount lateral stiff- 


ness. The fact that the propeller is not rigid may be 
an alleviating factor. 

The analysis can be extended directly to cases of two 
engines running and more wing modes if the higher de- 
gree equations that result are no barrier; what is needed 
in this problem is merely the additive collection of two 
gyroscope kinetic potentials plus the other usual ener- 
gies in analogous fashion to Eq. (13), from which the 
equations of motion are subsequently obtained as be- 
fore. 

The casting of the gyroscope problem in the x-y- 
coordinate system rather than the usual Euler coordi- 
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Approximate Solution of Compressible 


Fluid-Flow Problems 


by Galerkin’s Method 


CHI-TEH WANG* anp R. F. BRODSKY? 
New York University 


SUMMARY 


The Galerkin’s method has been used in the literature for the 
approximate solution of elasticity and vibration problems. In 
this paper, its application to compressible fluid-flow problems 
is carried out. It is found that, for compressible flow problems, 
Galerkin’s method has two advantages over the Rayleigh-Ritz 
method, which was previously worked out by one of the authors. 
First, it is not necessary to formulate variational integrals for 
various new problems; therefore, these problems can be treated 
without any mathematical difficulty. Second, there is no re- 
striction on the value of the ratio of specific heats, y, to be used, 
a fact that is necessary in the discussion of some irrotational flow 
problems.‘ The nonlinear problem of compressible flow passing 
a circular cylinder is carried out as a numerical example. The 
results are found to be in good agreement with other methods. 
An extension of Galerkin’s method to the general case of flow 
passing arbitrary bodies is briefly indicated. For these irro- 
tational flow problems, however, Galerkin’s method seems to 
suffer the disadvantage of requiring more work in carrying it out 
than the Rayleigh-Ritz method. The application of the Galer- 
kin’s method to laminar boundary-layer problems is discussed, 
and the boundary-layer problem of a compressible flow past a 
flat plate is formulated. A numerical example is carried out to 
illustrate the method. 


SYMBOLS 


a = velocity of sound 

by = specific heat at constant pressure 
é = specific heat at constant volume 
Cy = skin friction coefficient 


H,h = enthalpy 


I = the variational integral 
k = coefficient of heat conduction 
M == Mach Number 
p = pressure 
R, = 1p:X,/u1, Reynolds Number at x = x, 
R = specific gas constant 
r,@ = polar coordinates 
\ i = absolute temperature 
u = velocity component in x-direction 
U = free-stream velocity 
v = velocity component in y-direction 
x,y = rectangular coordinates 
= ratio of specific heats = c,/c, 
M = the coefficient of viscosity 
p = density 
a = Prandtl Number = wcp/k 
T = viscous stress 
¢ = the velocity potential 
w = defined by ux 7° 
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of Engineering. 


Subscripts: 
0 = 
1 = 


cs = 


free-stream conditions 
quantities measured at outer edge of boundary layer 
quantities measured at wall 


INTRODUCTION 


- PREVIOUS PAPERS,'”? one of the present authors 
has applied the variational method to the com- 
pressible fluid-flow problems following the Rayleigh- 
Ritz procedure. The results show surprisingly good 
agreement between the variational method and the 
other approximate methods. A study of the results 
also reveals that the variational method gives good ap- 
proximate solutions at both high and low Mach Nun- 
bers and for flows past both thick and thin bodies. 
The variational method has been applied to steady, 
nonviscous, irrotational flow past arbitrary bodies with- 
out any difficulty. A variational procedure for tran- 
sonic flows with shock waves has been proposed,’ but 
numerical examples are yet to be worked out. For 
viscous and other flows, it is necessary to formulate 
new variational principles that are not always easy to 
accomplish. 

On the other hand, in carrying out the variational 
method, the ratio of the specific heats, y, has to be 
such that y/(y — 1) is an integer. For simplicity in 
computation, y = 2 was used in such calculations. 
This assumption leads to a slight increase in the final 
velocities’ but has practically no effect on the final 
pressure. Nevertheless, this is a point open to objec- 
tion, especially when it is necessary to use some con- 
stants in which y plays an important role.‘ 

In order to overcome these difficulties, Galerkin’s’ 
method has been applied by the authors to compressible 
flow problems.* Galerkin’s procedure does not require 
the formulation of any new variational principles, and 
there is no restriction in the value of y to be used. 
The only difficulties appear to be the choosing of a 
function that satisfies the boundary conditions for 
flow past arbitrary bodies and the method of solution 
of the resulting nonlinear simultaneous algebraic equa- 
tions. Methods to surmount these difficulties have 
already been discussed in detail in references 1 and 2 
When the authors first worked on the method,® it was 
thought that the application of Galerkin’s method to 
compressible flow problems was a new venture. It was 
discovered later that the method was used previously 
by Akita.’ However, it appears that the method was 
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incorrectly carried out in Akita’s paper because of an 
oversight in the integration process. 


GALERKIN’S METHOD 


The Galerkin method has been used in the literature 
jor the approximate solution of elasticity and vibration 
problems. However, in these problems, other meth- 
ods for their solutions are available, and Galerkin’s 
method has not received wide attention. 

Consider the case of a steady, irrotational compres- 
sible flow; the proof for Galerkin’s method may be 
carried out as follows: 

In reference 1, the problem of taking the first vari- 


ation of the variational integral 


[ = | / F(o,, ¢,)dx dy (1) 
R 


and setting 6J = 0 wasseen to lead to 
IS Li¢)n(x, v)dx dy = 0 (2) 
R 


(OF, /oy) = 0 (3) 


where 
L(o) = —(OF, /ox) 


is the Euler equation associated with the variational 
problem and n(x, y) is an arbitrary function or the first 
variation of @. After differentiation and simplification, 


Eq. (3) becomes 


dy" )byy — 
2o:,Or, = 0 (4) 


L(¢) = (a? — $,7)¢,; + (a? — 


where 

a? = ao? + [(y — 1)/2] (LU? — ¢,? — ¢,”) 
In the above equations, ¢ is the velocity potential, a 
is the velocity of sound, L’ is the free-stream velocity, 
the subscript 0 refers to quantities at infinity, and the 
subscripts x, y denote partial differentiation. 

Eq. (2) can be interpreted as stating that the Euler 
expression L(@) is orthogonal to every variation en(x, y), 
where ¢€ is an arbitrary small parameter. Assume that 
¢ can be approximated by the series 

\ 


ov(x, ¥) = Do a filx, y) (5) 
1 


' 
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where the functions f;(x, y) satisfy the boundary con- 
ditions imposed on ¢(x, y) and ¢; are undetermined 
coefficients. Eq. (2) will, in general, no longer be valid 
for an arbitrary variation en(x, y), since this would 
imply that ¢x(x, y) satisfies the Euler Eq. (4) and, hence, 
is the exact solution of the boundary-value problem. 
Although the variation in / does not vanish for an 
arbitrary variation en(x, y), the constants ¢; in Eq. 
(5) can be so chosen that the variation in / is zero for the 
N variations ¢f;(x, y). If, now, N becomes infinite 
and the functions f;(x, y) form a complete system of 
functions, then the set of all relations 


SS Liox)fidx dy =0 (i = 1,2,..., 0) (6) 


becomes equivalent to the arbitrary relation (2) and 
ox (x, y) becomes an exact solution. 

The above proof is essentially the one given by 
Sokolnikoff.* Duncan,® on the other hand, has given 
an interesting proof that shows that the Galerkin’s 
method can be derived from the method of least mean 
square error. Duncan's proof, however, is valid only 
in the case of linear problems. 

The Galerkin process may now be summarized as 
follows: First, assume a series of functions that satisfy 
the boundary conditions exactly but with undetermined 
parameters. For example, in subsonic irrotational 
compressible flow, the velocity potential @ may be ex- 
pressed as follows: 


1) = do + | > A mnPmn (7) 


m n 


where ¢o and ¢,,, satisfy the boundary conditions and 
Am, are the undetermined parameters. Substitute 
the assumed expression into the differential Eq. (4). 
Eq. (4) vanishes if Eq. (7) is the exact solution of the 
problem. Otherwise L(¢y) is a function of x and y and 
therefore can be taken as a measure of the error asso- 
ciated with the approximation ¢y. Multiply L(@y) by 
mn, integrate the resulting expression over the domain, 
and equate to zero. If there are N coefficients, N alge- 
braic equations are obtained which may be solved 
simultaneously, and the coefficients A,,, are thus 


determined. 


IRROTATIONAL COMPRESSIBLE FLOW PAst A CIRCULAR CYLINDER 


As an example, the problem of compressible flow past a circular cylinder has been carried out. 


The fundamental 


differential equation of a compressible fluid in polar coordinates is as follows: 


os i le tet Pe » oe oo | 
L(¢) = & +2 (u: — $7 - *) - 62] rr + a ‘= (. —~— ) 7 4 -* 
3 a . : 


ems l y s de” by” d, dre , 
| a +- _ (1 a -) + | a 26,6 “ie O (8) 
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where the subscripts 7, @ indicate partial differentiation. The boundary conditions are 


[¢,2 + r—°o,?]'"* = Uatr = o and ¢, = Oatr = 1 (9) 
The radius of the cylinder is taken to be unity. To satisfy the boundary conditions, ¢ may be assumed as 
? l K : : l 1 : 
@= U({rt+ ) cos @ + den >» y = - | (Amn cos 00 + By, sin n6) (10 
r 2r m=1 n=1L™r (m -+- 2)r"™t 2 


where x is the strength of circulation and A,,,, B,,, are undetermined parameters. Then, the Galerkin equations 
for determining the parameters A ,,,, and B,,, are 


"r= a "O=2e \ 
T | 1 ) 
/ | L(y) | _ -|cos nOrdrd6=0 | 
r=1 6=0 mr” (m a 2" +2 
a x "9=27 1 i 
L(y) — | sin nO rdrd@=0 | 
-— o=0 mr™ (m -+- 2)r™*? } 


where L(y) is the resulting Eq. (8) when the assumed ¢ with NV parameters [Eq. (10)] has been substituted. If 
the assumed ¢ is the exact solution of the differential Eq. (8), then L(@y) is equal to zero. In general, of course, 
¢ cannot be the exact solution; then, L(@y) will not be zero. 

In the case of the flow past a circular cylinder without circulation, the flow must be symmetrical. The un- 
symmetrical terms in ¢ can therefore be discarded. Taking four parameters, we have 


] l l ] l 
@=U (- + - ) cos 6+ Ay ( —- :) cos 86 + Ais ( _ : cos 30 + 
r r 3r® r 3r3 


l l l 
Ax ( o> a, oe ) cos 6 + A33 ( oe :) cos 36 (12) 
or* 5r° ors ee \ 


Eqs. (11) then become j 


1 1 ei I l ) 
IS L(¢4) ( —-- ;) cos 6r dr dé = 0, / / L (4) ( —_ ) cos 30 r dr dé = 0 
r 3r8 r 3r3 | 


(11) 


: . p (13) 
; 1 l l l 
L(¢s) | == — =~) cos Or dr dé = 0, L (4) — ~~) cos 30rdrd0=0 ! 
ors 5r5 ; 3r? a J 
Carrying out the integration and taking y = 1.405, the following equations are obtained: 
(0.51851844a9? — 0.6819993U2)An + 0.4205556U2Ai3 + ) c 
" 
(0.0814814a9? — 0.0902122U?)A3, + 0.0271777U?A33 = Ky, | t] 
+ 0.7841668U2A1n + (3.3333328a9? — 3.5691657 U2) Ais + n 
| tl 
0.1182662U*Aa + (0.3999998a9* — 0.4778126U")Ass = Kz | ¥ 
; (14) | 
(0.081481 4a? — 0.0928529U2)Ay, + 0.0772793U?Aj3 + t 
(0.0225185a9? — 0.0152367 U*)A 3 + 0.0081639U"A;; = Ky 
{ a 
+ (.0972529U? An + (0.3999998a92 — 0.4950988 U2)A 13 + 0.0174196Aq + | | 
(0.0693329ao? — 0.0789673U2)Ass = Kg } 
1 
— or k 
where Kj, Kz, K3, and AK, are sums of nonlinear termsin Use the values of A’s obtained to compute K’s, and solve d 
Au, Ais, Asi, and A;3; their complete expressions are for these parameters again. This gives the second ap- t 
given in Appendix I. The order of Eqs. (14) follows proximation. Repeat the cycle until the values of ( 
from that of Eqs. (13), respectively. These equations A’s computed in two consecutive cycles are the same. | 
can be solved as follows: First, assume the values of An example of the calculation may be found in refer- | é 
A’s to be zero in the K’s. Eq. (14) then becomes a encel. At My = 0.4, the solution of Eq. (14) is 
system of linear simultaneous equations which can be 
solved by the Crout method.” The values of A’s so An/ao = 0.09664, A13/ao = —0.02194 | (15) 7 
° ° , P sa a 0 ; 
computed may be taken as the first approximation. A3:/do9 = —0.04563, A33/ao9 = —0.01975 | 


PROBLEMS 663 





SOLUTION OF COMPRESSIBLE FLUID-FLOW 
TABLE 1 
(9 Velocity Distribution on the Surface of a Circular Cylinder at MJ) = 04 
9) a _—_——$3. $$ $$ 
Variational Rayleigh-Janzen 
Galerkin Method (y = 1.405) ; Method Method 
q/l Linearized y¥=2 y = 1.405 
ry 1 Term 2 Terms 3 Terms 4 Terms solution (4 Terms) 
(10 0 0 0 0 0 0 0 0 
10 0.3714 0.3125 0.3119 0.3078 0.3106 0.3022 0.3104 
20 1.0696 0.9534 0.9517 0.9435 0.9483 0.9347 0.9588 
10Ns 45 1.5127 1.5108 1.5088 1.4956 1.4962 1.4937 1.4154 
60 1.8526 1. 8608 1.8584 1. 8583 1.8541 1.8663 1.8355 
80 2.1068 2.2208 2.2183 2.2254 2.2142 2.2460 2.2269 
90 2.1393 2.2697 2.2672 2.2753 2.2632 2.2979 2.284 
11) 
TABLE 2 
Velocity Distribution at 6 = 2/2 from Various Methods 
If Variational Taylor Rayleigh- 
i Galerkin Method (y = 1.405) Method Experimental Janzen 
se, q/l Linearized y¥=2 Method Method 
r/b 1 Term 2Terms 3Terms 4Terms solution (4 Terms) y = 1.407 y = 1.405 
1.0 2.1393 2.2697 2.2672 2.2753 2.2632 2.2979 2.284 
un- 1.04 2.0582 2.1833 2.1810 2.1887 2.1770 2.210 2.13 2.197 
1.2 1.8059 1.9104 1.9099 1.9144 1.9052 1.931 1.89 1.920 
1.5 1.5235 1.5976 1.5990 1.6000 1.5940 1.610 1.60 1.603 
2.0 1.2979 1.3427 1.3412 1.3437 1.3406 1.349 1.35 1.345 
3.0 1.1335 1.1544 1.1560 1.1547 1.1534 1.157 1.18 1.155 
2) Tan S As mentioned previously, the velocity distributions 
Values of Amn When One, Two, Three, and Four Parameters Are computed by the Galerkin procedure show good 
Taken : ; . ‘ 
i . $ agreement with those obtained by the Rayleigh-Ritz 
Ay Ais Au A procedure. This confirms the fact that the use of 
| 0.0836000 i . . 
0.0892318 —0.0242009 y = 2 in the computation has only a slight effect on 
0.0969708 —0.0242472 —0.0469989 igi : il it eo _ eer 
3) 0.0966415 —0.0219423 —0.0456325 —0.0197478 the magnitude of the velocities. Now let us compare 
—_——_ the necessary amount of computational labor in these 
: 0.09352 —0.02158 —0.04462 —0.01425 two methods. It is found that Galerkin’s procedure 
requires much more work, especially when flow past an 
The velocity distribution on the surface of the @tbitrary body is to be studied. Therefore, in the case 
cylinder and that at @ = x/2 computed by Galerkin’s of irrotational steady flows past arbitrary bodies, the 
method are tabulated in Tables 1 and 2, together with variational method following Rayleigh-Ritz procedure is 
c c c c » « “~y) ¢ 
the Taylor! electrical analog results, the third approxi- recommended. If only the velocity and pressure dis- 
mation of the Rayleigh-Janzen method by Imai,'? and tributions are to be computed, the linearized variational 
‘ é anz j 4 
the variational method following Rayleigh-Ritz pro- integral may be used to give good approximate results 
, cedure.'' 2 The results show close agreement with with only a small amount of work. This fact has al- 
ready been discussed in reference 4. As an additional 


those previously obtained. In order to get an idea of 
the convergence of the series, the velocity distributions example to show that the method of linearization sug- 
at these points are also computed by taking only one 8¢sted in reference 4 also gives good results in the pres- 
and three ¢"t approach, the velocity distributions computed 


parameter A,,; two parameters Ay, Ai3} 


from the linearized solution are also included in Tables 
1 and 2. It can be seen that the agreement between 


these results is good.* 


parameters Ay, Ai3, As. It may be mentioned that 
if only A; is taken, the first equation of Eqs. (14) must 
be used for solving for Ay. If two parameters, A, and 
A};, are taken, the first two equations of Eqs. (14) are 
to be used, etc. The values of A,,,, are different when 
different combinations of parameters are taken. At 


My = 0.4, the parameters A »,,/do are given in Table 3, 


LAMINAR BOUNDARY-LAYER PROBLEMS IN 
COMPRESSIBLE FLow 


The motion of a compressible fluid in laminar 


——— 


and the velocity distributions computed by taking one, 
two, and three parameters are tabulated in Tables 
1 and 2 for comparison. It appears that the rapidity 
of convergence of the values for velocities is satisfac- 
tory. 


boundary layer is to be determined from the following 
governing equations: 


* The authors are indebted to Messrs. P. C. Chou, R. N. House, 
Jr., and J. Van Name for their help in carrying out and checking 


the integration. 
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(1) Equations of motion, 


puu, + ply = (uty), — pz (16) 
pb, = 0 (17) 
(2) Equation of continuity, 
(pu), + (pv), = O (18) 
(3) Equation of energy, 
pull, + pv, = [(uH,),/o] + wu,? + up, (19) 
(4) Equation of state, 
b = RoeT (20) 


where y is the coefficient of viscosity, H = fc, dT is 
the specific enthalpy, c, is the specific heat at constant 
pressure, JT is the absolute temperature, o = uc,/k is 
the Prandtl Number, & is the coefficient of heat con- 
duction, and R is the specific gas constant. For a per- 
fect gas, c, is a constant and H = ¢,T. Eqs. (16) 
to (20) are five equations in terms of the five variables 
u,v, p, p, and 7. 


The coefficient of viscosity u is a function of T only. 
An approximate relation between 7° and uw may be 
written as 


pa T° (21) 


where w varies according to the range of JT in which we 
are interested. By means of Crocco’s transformation, '* 
Eqs. (16) to (20) may be reduced to the following two 
equations in the variables 7 and r: 


u(pu/T)> + Tu 7 pPr(u lt) 4 = 0 (22) 
‘i — Hun - 
( *) Hytr, + = Tt wt = upplly + 
o o 


u(HT, + u) pz = 0 (23) 


where tr = pu, A summary of Crocco’s work may be 
found in a report by Young." 


In the case of the flat plate at zero incidence, Eqs. 
(22) and (23) may further be reduced to ordinary equa- 
tions in terms of a single independent variable. This 
variable is of the form § = cy/+/2x, where c is a con- 
stant. In this case, we may write H = H(u), p = 
p(u), uw = u(u), and 


r= pu, = put t, = f(u)/V/2x 
Eqs. (22) and (23) 


(24) 


where f(u) indicates a function u. 
therefore become 


ff” + upp = 0 (25) 
(H" + of + (1 —o)H’f’ =0 (26) 


where primes denote differentiation with respect to 
u. Note that p, = Ointhiscase. Then, introduce the 
following nondimensional parameters: 


u/u = %, H/Hi=h, p/n =3 
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p/n =f FR) = f(uyV2 ‘piu? 


where the subscript 1 refers to quantities in the maip 


stream. From the equation of state [Eq. (20)], we 
have 

B= p/m = 1\/T = M/H = 1/h 
and 


B= w/m = (T/T;)* = h* 


With these relationships, Eqs. (25) and (26) may be 
put in the following form: 


FF" + 2uh “"® =0 (27 
(h” + ou,?/M)F + (1 — o)h’F’ = 0 (28) 


where the primes now denote differentiation with re- 
spect to %, noting that 


u;" ‘Hy, = (y = 1)M,? 
The boundary conditions are 


F’ = 0, h = H,/H, (when ai = 0) \ 
F =0, 1 (when a = 1) J 


29) 
h= 
where the subscript w denotes quantities measured at 
the wall. 


By a laborious trial-and-error method, Crocco inte- 
grated Eqs. (27) and (28) step by step for different 
values of ¢, w at various M;. Approximate solutions, 
however, can be more easily obtained by the Galerkin 
method. In this case, Galerkin’s method may be car- 
ried out as follows: First, assume F in terms of a series 
of functions that satisfy the boundary conditions but 
with undetermined parameters. Then, Eq. (28) be- 
comes a linear ordinary differential equation in / and 
can be solved analytically. Now, the solution of h is 
obtained which contains the undetermined parameters 
assumed in F. These parameters can be determined 
from Eq. (27) by means of the Galerkin method. 


Define the skin friction coefficient cy to be 274/ piu’; 
the local skin friction coefficient at a point distant x 
from the leading edge of the plate may be computed 
from the following equation: 


cy = F(0)/VR, (30) 


where R, = up:x/, and F(0) denotes the value of F 
at “7 = 0. Eq. (30) may be easily derived from Eq. 
(24). Once F(a) is determined, the relation between 


y and u is given by Eq. (24), and, hence, 
y= JS udu/r(x, u) 


which reduces to 


Pt 
mpiy/m = 2/R; f h*di/ F(a) (31) 

As a numerical example, assume 
F = A(1 — a?)'” (32) 


Eq 
tu 
ob 


wi 


an 


ao 


Main 
], we 


29) 





SOLUTION OF 


where A is an undetermined parameter. Note that 
Eq. (32) satisfies the boundary conditions. Substi- 


tuting Eq. (32) into Eq. (28) and solving, one 
obtains 
— uy? p 
h=1-—h'(0) Cl4é,o) +o D(a, o) = (33) 


iM, 


where 


1 
Clu, oc) _ I. (1 _— u?)‘*~ 1) /2 di 
and 


1 
D(u, a) = i (1 a ui)" 1)/2 x 
i (1 — 42)°- 9 aa] di 


The parameter A may be determined from the Galer- 
kin method by substituting into Eq. (27) the assumed 
expression for F and the expression A given by Eq. 
(33), by multiplying the resulting expression by (1 — 
a?)'”, integrating over the region of a, then setting the 


resulting expression to zero. Therefore, we obtain 


1 
A? = | f. oak “~~ (1 — a)” aa | : 
0 
ne 
| / (1 — wy da] (34) 


For various values of o and w, Eq. (34) can, in general, 
be integrated by numerical methods. 

Ifo = landw = 1, Eq. (34) can be easily integrated 
to give 


A = 0.652, c = 0.652/+/R, 


The value of A computed by the present method is seen 
to compare favorably with the value 0.664 arrived at 
by Crocco for one side of the plate. 


COMPRESSIBLE 
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Ky = 0.5277779 UL? + 0.1578077 UA? + 0.0114401A1*? — 0.0793961A 133A 1? + 0.0031080A9An? — 
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0.4331232 UA 13A 33 + 0.0506732A 1A 13A33 — 0.0144869UA 31A33 + 0.0298569A 1A 134A 33 + OU ao? 


yy 


2 = —0.8333333U% — 0.1672530UA 1? — 0.0085908A 13 + 0.2978344A 9A 1? — 0.0002069A 1A 11? + 


0.0389466A 33A 1? + OUA13? + OAA13? + 1.2948582A13* + 0A 91A13” + 0.6152349A 33413? + 0.0002484UA 3:7? — 
0.0165082A 1A a1? + 0.0106790A 134 12 + 0.0001667A 3:3 + 0.0013032A 3A 1? + OUA33? + OA1A3;” 

0.2591467A 13A 332 + 0A 3A 33” + 0.0054076A 33* + 1.9752402UA11A13 — 0.0380053UAnAn + 

0.2912086 UA 1A33 + 0.2862628UA 3A 3: + 0.1029188A 1A13As + 04134 9A33 + OUAA33 + OAnAAss + 


0.0399478UA3,A 33 + 0.0501962A 114 13A 33 + Ol lay? 
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K; = 0.0788884U? + 0.0229610UA1;2 + 0.00258544 1? — 0.0108586A 134 u2 + 0.0101780A 34 12 — 
0.0108677A 3A 1? + 0.1835371 UA 13? + 0.0387298A 1A 13? + 0Aj53 + 0.1617810A 3A 15? + 0A y 4,3? — 
0.0001661 UA s:2 + 0.0002916A 1A s12 —0.0007270A 13A 1? + 0.0000385A 3:3 — 0.0001291A Ay? 4 
0.0058676 UA 33? + 0.0009858A 1A 33? + 0A13A 33? — 0.0001042A 914332 + 0A35* — 0.0539584UAy Ay + 
0.0052325UAnAs — 0.0097463UA1A33 — 0.0114135UA13A 3 + 0.0479104A 11.4154 3) — 0.0005767A 3A Ay 4 
0.0654494 UA 13A33 + 0.0113346A 1A 13433 — 0.0028457UA51A33 — 0.0010434A 114 1A 33 + OUae 


K, = —0.1000001U* — 0.0276360UA1;? — 0.0017326Ai,° + 0.0500488A 134A 1? — 0.0090943A 3:A 11? + 
0.0260358A 33A 1? + OUA13? + 0A1A13? + 0.2164743A 13 + 0A 31A13? + 0.1892152A :A},? ~ 
0.0003166U'A 3:7 + 0.0002075A 1A 317 + 0.0020334A 313A 17 + 0.0000382A 31° + 0.0003830A 33A 9)? + 
OUA3;3? + 0A11A33? + 0.0186546A 13A 337 + 0A 31A 337 + 0.0010804A 33° + 0.3104708UA A); - 
0.0081622UA As: + 0.0497816UA :A33 + 0.0484180UA 134A 3) + 0.01842744 1A 3A 31 + 0A 123A 31Ay + 
0U'A\3A33 + 041A 13A 33 + 0.0073674UA 3A 33 + 0.0031616A 1A 1A 33 + OU ay? 
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The Gyroscopic Effects of a Rigid Rotating Propeller 
on Engine and Wing Vibration Modes 


(Continued from page 659) 


nates is advantageous for other problems. For example, 
the effect on aircraft stability in flight of the high-speed 
rotor of a jet engine can readily be treated. Whether 
these effects will be large or small has not been investi- 
gated by the present authors. Could gyroscopic 
moments be a partial cause of some “‘snaking”’ effect in 
jet aircraft observed while trying to hold them on 
course? 

The application of the considerations of this paper 
to other spring-mounted, rotating machinery is also 
direct. It should be noted, in passing, that the rotor 
problem treated here is not the same as the well-known 
one of the gyroscopic effect on critical speeds of over- 


hanging shafts, etc. In such problems, rotation and 


precession are usually locked together, while nutation 
is nonexistent. 
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First U. S. National Congress of Applied Mechanics 


Plans have now been formulated to hold the First U.S. National Congress of Applied Mechanics on 
The Illinois Institute of Technology will be host to the Congress. 
Members of the Institute, one of the sponsoring societies, are cordially invited to submit papers, 
which should constitute original research in applied mechanics. 
words or the equivalent in equations, tables, and diagrams.) Abstracts must be submitted to the Chair- 
man of the Editorial Committee before April 14, 1951; complete papzrs must be submitted prior to June 


Additional information about the Congress may be obtained by writing to N. M. Newmark, 
Secretary of the Congress, University of Illinois, Urbana, III. 


(Their length must not exceed 5,000 


All members are invited to attend. 
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On Boundary-Layer Growth in a Compressible 
Flow Between Radial Walls 


J, Conrad Crown 
U.S. Naval Ordnance Laboratory, White Oak, Silver Spring, Md. 


June 28, 1950 


i ig SOLUTION FOR the laminar flow field about a flat plate in 
a uniform stream is well known. However, it is also of 
€.g., ones in which streamwise 


We shall consider here the com- 


interest to study other flows 
pressure gradients are present. 
pressible laminar radial flow between radial, two-dimensional 
walls! and through an axially symmetric conical nozzle. In par- 
ticular, we shall concern ourselves with the displacement thick- 
ness of the boundary layer. 

For two-dimensional steady flow, the equation of continuity in 
polar coordinates is 


(Opur/dr) + (Opv/d00) = 0 


where u and v are the components of velocity in the directions 
rand 6, respectively, and p is the density. 
If the flow is radial,v=0. Then 


pu = f(0)/r 


where f(@) is so far undetermined except that it is a function only 
of 6, Reynolds Number, and wall divergence angle. Let @ = 0 
be the centerline of the flow and such quantities be designated 
with a subscript “1.” Let @ = +a be the boundary along which 
u=U. 

The displacement thickness is defined by 


a 1 
5* =f (: ia mY yap 
0 pill, 
and for radial flow becomes 
sc =r ff. {1 — [f(6)/f(0)]} de 


Thus it can be seen that the boundary-layer displacement thick- 
hess grows linearly with the radius. The proportionality con- 
stant depends on the Reynolds Number and wall divergence 
angle. 

Similarly, if we consider radial, axially symmetric flow through 
a pipe or nozzle, we find that the displacement thickness also 
grows linearly with the radius. 

It is interesting to note that, in a supersonic nozzle (with radial 
walls), in spite of the favorable pressure gradient, the boundary 
layer grows faster than over a flat plate with no gradient. In 
this comparison it should be borne in mind that in the nozzle the 
displacement thickness is based on centerline flow characteristics 
that vary along the nozzle while that of the flat plate is referred 
to free-stream conditions that are constant. 


No restrictions have been imposed here on Mach Number. 
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On Virtual Mass and Wave Drag at Subsonic 
Speeds 


John W. Miles 
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University of 


po PROBLEM OF IMPULSIVE MOTION of a body in an ideal, in- 
compressible fluid has been discussed by Lamb, Larmor, 
Kelvin, and others.!. As might be expected from the concept of 
virtual mass,? an impulsive force is required to generate impulsive 
motion, insofar as the fluid is assumed to be incompressible. 
However, if the compressibility of the fluid is taken into account, 
the force required to generate impulsive motion may be expected 
to remain bounded and to vanish only asymptotically with time. 
Moreover, the concept of virtual mass, while still valid for the 
description of the energy and momentum of the fluid flow in the 
steady state, is of no direct aid in computing the transient re- 
sponse. In particular, we find that the work required to generate 
impulsive motion is not entirely recoverable, so that the process 
is irreversible, and some part of the starting force must be re- 
garded as ‘‘wave drag.’’ (Indeed, the present remarks were 
motivated by a brief discussion of just this question with Prof. 
Th. von Karman, who raised the question of the validity of the 
virtual mass concept in a compressible fluid. Our thanks also 
are due to Prof. H. K. Forster, with whom we discussed the same 
question. Professor Forster pointed out the similarity between 
the present situation and the acceleration of an electron, with 
which is associated radiation.) 

As a simple example, we consider certain aspects of the sudden 
acquisition of a velocity v by a body in an ideal, nearly incom- 
pressible fluid. First, after equilibrium is attained, the flow will 
be potential, and there will be no resultant force on the body 
Second, if m is the virtual mass of the body in this potential flow, 
the final momentum of the associated flow will be mv, whence the 


total impulse required to impart the motion is given by 


I= fi Poa = mv (1) 


(Of course, if the fluid is truly incompressible, the time required 
is zero and the force infinite, but this indeterminacy does not affect 
the integrated results.) Third, since this force acts at a constant 


velocity, the total work done is given by 


W = f, F (t)v dt = mv? (2) 


Fourth, the kinetic energy of the associated flow, after equilibrium 
has been established, is given by 

K.E.: = (1/2) mv? (3) 
It follows that an equal amount of energy [(1/2) mv*] has be- 
come unavailable during the period of transient motion. More- 
over, if the body were to be stopped suddenly, the remaining 
(1/2)mv? also would be rendered unavailable. The agency of 
dissipation, in both cases, is radiation, and the associated force is 
appropriately designated as ‘‘wave drag,’’ although this term 
generally is applied only to supersonic configurations. 

We remark that a result somewhat analogous to the foregoing 
arises in the rather more transparent problem of a linear system 
of one degree of freedom. Thus, if an ideal capacitor is charged 
by a suddenly applied constant voltage source, the final energy 
stored in the capacitor is only half that delivered by the source, 
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even though there is no (ohmic) resistance in the circuit. In 
both cases, we have irreversibilities, albeit without direct degra- 
dation of energy. 
In seeking an explicit solution to the foregoing problem, di- 
mensional consideration leads us to write 
F(t) = mv(c/))f(ct/l) (4) 
where / is a characteristic length and (l/c) is a characteristic time, 
the characteristic velocity being chosen as the sonic velocity c. 
Similarly, if a constant force F is suddenly applied to a body, we 
write for the resulting acceleration 
(d/dt)v(t) = (F/m)g(ct/l) (5) 
In the case of a circular cylinder of radius a, moving normal to 
its axis, we find on the basis of the acoustic approximation 


(v/c< 1) 


1 5 
f(x) =1— o* — _ +7" + 0(x*) (6a) 
= —2x-3 — [62 — 24 In (2x)]x7-5 + O(x77 In? x) (6b) 
g(x) = (x) + : + - + O(x?) (7a) 
= 1 — x>? + [4.232 — 3 In (4x)]x~4 + O(x-* In? x) = (7b) 
where 6(x) is the Dirac delta function, and the characteristic 


length is a. It is of particular interest that f(x) changes sign. 
In connection with g(x), we remark that the actual apparent mass 
is initially zero, jumps immediately to 2m, decreases in accordance 
with 2m[1 — (37r/4)] and is asymptotic to m as (1 + 7r~?), where 
t = ct/a. More detailed calculations are in progress but await 
the availability of certain tabulated integrals, to which we have 
been referred by Prof. A. Erdelyi. 

If a constant force acts for time 4 on a circular cylinder, we 
find that the work done for small and large values of 7; is given 


by 


. my 1 l 
W = [27,-! + — +-—7 + O(71?) | (8a) 
2 2 8 
mv? bs 
= 2 {1 + 27; 2 In (271) + Tl ‘— 2774 In (4 T1) + 
O (7176 In? 7) (8b) 


where v is the final velocity. Hence, if a body is accelerated 
rapidly, the starting efficiency is vanishingly small; while, if it is 


accelerated slowly, the efficiency approaches unity. 
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Effect of Wing Parameters on Load Intensity 
in Wing Skin 


B. E. Gatewood 
Department of Mechanics, U.S.A.F. Institute of Technology, 
Wright-Patterson Air Force Base, Dayton, Ohio 


June 19, 1950 


TT” LOAD INTENSITY /, (Ibs. per in.) is extremely useful in 
preliminary wing design as a guide to the type of wing struc- 
At any wing section 


P; = M/hd (1) 


ture that should be used. 


where M is maximum beam moment on the section, / is the aver- 
age height of the structural section, and d is the length of the 
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structural section. The equivalent skin thickness for the top or 
bottom surface of the wing is P; divided by the allowable streg. 
Thus, in preliminary design, a large P; indicates the Possibility 
of using heavy skin without stringers, while a small P; indicates 
thin skin with stringers. Naturally, the shear flow q% (Ibs, per 
in.) from the torque 7, where 


gq = T/2hd (approx.) (9 


must be considered in conjunction with P;. 

Since in many cases in preliminary design the bending moment 
M in Eq. (1) is not known or requires considerable calculation to 
obtain, it is desirable to express M (and, hence, P;) in terms of the 
wing parameters. This can be done without difficulty if the 
section lift coefficient is assumed constant and the wing dead 
weight section coefficient is assumed constant. The P; thus ob. 
tained is sufficiently accurate for preliminary design. The for. 
mula for P; is 


P, = 


W 
(2) can (1+ A) [1 — (2y/b)]? [1+ 2A- (1 - ) (2y/b) 


$04 Ce 
576 ( \ ) [1 — (1 — A) (2y/b)]? 
c/X\e 


where 
y = distance along semispan from fuselage centerline 
b = structural span (ft.) 
ny = taperratio = cr/cr 
cr = tip chord 
CR = root chord 





P; (KAIF Fy E 2a-ti-ar® | 
nb(IAR 576 (48) $8) [ (I-A) |? 
Pi 


nb()AR 
0.00 9.04 0.08 O12 O16 O20 


L's) 
je) 
ee. My 
a A 
u c RS 











150 














Ns ™ 
ee 


A) 


17 | 
ae 
<0 
a 
ig 
i 
SO 





% 
a 





Ry 
MR 








Rar’ 
Poy 





ra 
o= 
x 








= 


| 
0.0 O02 0.4 » os 08 1.0 


FIG.1 EFFECT OF WING PARAMETERS 
ON LOAD INTENSITY 
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FIG.2 EFFECT OF WING PARAMETERS 
ON WING WEIGHT 


W/S = wing loading (lbs. per sq.ft.) 

A.R. = structural aspect ratio = b?/S 

n = effective wing load factor (correct airplane load factor 
for tail load and wing dead weight) 


( = chord 

t,/¢ = effective structural thickness ratio 

¢,/¢ = effective structural chord ratio 

P; = load intensity (lbs. per in.) on upper or lower surface 
of wing 


Fig. | shows a nondimensional chart of P;/nb(W/S)A.R. 
against 2y/b with \ and ¢scs/c? as parameters. 

Since the equivalent thickness of the structural upper or lower 
wing skins is P;/Fay., the weight S,,, of the structural part of the 
wing to take the loads P; (i.e., the beam bending moments) can 


be expressed in terms of the wing parameters 


*b/2 P; 
Ws, = 4D c, dy (4) 
0 Fait. 


where 
P; = P;in Eq. (3) 
Ce = structural chord 


material density (Ibs. per cu.in.) 


1 ( 1 . 1 ) 
Fou. = =2\ Fou, Fall. 


Fait, allowable stress (Ibs. per sq.in.) in upper surface 


Fou, = allowable stress (lbs. per sq.in.) in lower surface 


W 


Ml 


s» = Structural weight of wing (Ibs.) 
The total wing weight will include another component for the 
nonstructural wing items, for the internal support structure (as 
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FORUM 


ribs), and for any webs to take the shear and torque loads. 
Eq. (4) gives upon integration 


W 
nDb* ) 
’ Ay 1 — 9d? + 8A? — BAF In A 
Ws, = (5) 


an te (1 —a)3 
3 Fall. 
" 


Fig. 2 shows a nondimensional chart of Fy. Ws,,/nDb*(W/S) 


against Fay /F with A and ¢,/c as parameters. Here, Fur. 


is the ultimate tensile stress of the material. 


The Basic Theorems in the Plastic Theory of 
Structures 


P. S. Symonds 

ay ety Laboratory, Cambridge University, England (on leave 
rom Brown University) 

July 26, 1950 


pres TWO THEOREMS DISCUSSED RECENTLY in the JOURNAL OF 
THE AERONAUTICAL SCIENCES by Neal! have been applied 
mainly to framed structures of mild steel whose members are 
assumed to exhibit the ‘‘limited strain hardening’’ characteristic, 
as it was termed by the author. Thus, the bending moment is 
assumed to be in the range between a ‘‘fully plastic moment”’ 
M, for positive bending and a corresponding value M,’ for nega- 
tive bending. At either of these values the entire cross section 
would be flowing plastically, with the material at the yield stress 
either in tension or compression. Based on this hypothesis the 
so-called ‘‘plastic methods” of analysis and design have been 
developed,?~* which have attractive features in their simplicity 
and deal directly with failure conditions of a structure. Thus, 
they eliminate some of the uncertainty in the choice of safety 
factor and make for greater economy in use of material. 

Although most of the applications have been to steel struc- 
tures, it should be possible to apply the plastic methods to 
frames composed of certain of the light alloys. The basic hy- 
pothesis that the bending moment cannot exceed a certain magni- 
tude (in one sense or the other), at which large changes in curva- 
ture occur at constant moment, is not strictly valid but would 
seem to be accurate enough for some alloys to allow useful re- 
sults to be obtained from the theory based on it. The differences 
in behavior between frames of steel and those of certain light 
alloys are probably matters of degree rather than of kind. The 
fact that the light alloys do not exhibit the flat yield character- 
istic of mild steel is not of major significance; the strain harden- 
ing at large strains which is neglected in the theory actually plays 
an important role in the behavior of steel frames* loaded to fail- 
ure in the sense of the plastic theory. However, a number of 
problems must be investigated, such as those of stability and of 
possible failure by rupture, before the plastic methods can be ac- 
cepted as fully valid for framed structures of light alloys. 

The writer has been especially interested in problems in whic h 
a number of loads on a structure may be repeatedly applied in a 
more or less independent manner. The criterion of safety here is 
taken to be that the structure is capable of ‘shaking down’”’ to 
some state (when loads are removed) with residual stresses and 
permanent strains such that the further application of any com- 
bination of loads in their prescribed ranges produces no additional 
plastic flow. 

The basic theorem governing this type of problem proved by 
the author can be put in a different form, which greatly facilitates 
the numerical work in solving a particular problem and also clari- 
fies its physical significance. The author showed that necessary 
and sufficient conditions for a structure to shake down are that 
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residual moments m; and parameters /;* exist such that 


mi? + m < M;* + (My); (la) 
Mim + m; > Mi* + (My’); (1b) 
M,* < (Mp): — (My); (1c) 
M;* > (Mp’)i — (My’); (1d) 


In these inequalities the subscript 7 specifies the location of a 
typical cross section of the structure; 9,”¢7- and I,” are 
the maximum and minimum bending moments at the 7th section 
which the given set of loads could produce if the frame behaved 
wholly elastically; (My); and (My’); are initial yield moments 
for positive or negative bending, respectively; and M;* is a 
parameter that depends on the previous straining history and is 
defined so that yield begins in positive bending when 
M; = M,* + (My); 
or in negative bending when 
M; = M;* oa (My’); 


A simpler form of the inequalities (1) can be obtained by elim- 
inating the parameter M;* from them. Consider, for example, 
inequalities (la) and (lb). These can be rearranged as follows: 

Mime + m; — (My)s < Mi* < Mm + m; — (My’); (2) 
When written in this form it is obvious that these two inequali- 
ties can be satisfied if, and only if, the left-hand member of in- 
equalities (2) does not exceed the right-hand member and, hence, 


if 


Mmez — "'™ < (My); — (My’); (3) 
Similar consideration of inequalities (la) and (1c) leads to 
Mm + mi < (My); (4a) 
while consideration of inequalities (1b) and (1d) yields 
Ms + my > (M,’): (4b) 


The only other pair of the original inequalities which can be 
combined in this manner are (lc) and (1d). These require the 
inequality 


(My); — (My’)s S (M,)is — (My 


to be satisfied. This is automatic in view of the definitions, and, 
hence, this result is of no further interest. 

Thus, the necessary and sufficient conditions for shake-down 
are seen to be the existence of residual moments satisfying the 
inequalities (4a) and (4b) and the satisfaction of inequalities 
(3), which depend only on the extreme values of the external loads 
and on the elastic moment range for each member. These sys- 
tems of inequalities are simpler than the author’s, and they can 
be used without difficulty in the numerical solution of specific 
problems.® 

The physical interpretation of the new inequalities is of par- 
ticular interest. 

There are two different and independent phenomena of con- 
tinuing plastic flow when shake-down does not occur.» ® The 
first and most obvious is that in which plastic flow in some re- 
gions is made to occur alternately in one sense and in the other, 
as appropriate load cycles are applied. (This would be a severe 
type of fatigue loading condition.) The second type of con- 
tinuing plastic flow is that in which definite amounts of plastic 
deformation recur always in the same sense, as certain load cycles 
are repeated. Thus, if a sufficient number of cycles were applied, 
deflections of arbitrarily large magnitude could be produced 
(under the hypotheses of the theory) without increasing the 
extreme values of the external loads. In practice, the unlimited 
growth of deflections is halted mainly by strain-hardening, but 
recent tests show that the theoretical critical load limits have 
very real physical significance. 

These two distinct types of continuing plastic flow are gov- 
erned, respectively, by the sets of inequalities (3), (4a), and 
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(4b). Thus, if loads are such that inequalities (3) are not Satis. 
fied, then plastic flow alternating in sense can be made to con. 
tinue indefinitely. On the other hand, if residual moments satis. 

fying inequalities (4a) and (4b) (and the equilibrium equations | 
cannot be found, then deflections of large magnitude ean be 
produced by sufficiently many repetitions of certain loading 


cycles. 
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Note on Plastic Column Behavior 


George Gerard 
Guggenheim School of Aeronautics, College 
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July 3, 1950 


of Engineering 


Srver AL DISCUSSIONS OF THE STABILITY of a plastic column 
which have been concerned with a model of the stress dis- 
tribution in the cross section at buckling have appeared in the 
recent literature. The most recent paper by Duberg and Wilder 
(JOURNAL OF THE AERONAUTICAL SCIENCES, Vol. 17, No. 6, 
pp. 323-327, June, 1950) attempts to go further by redefining 
the critical load for a plastic column “.. . based on the behavior 
of a slightly bent column as the initial lack of straightness 
vanishes.’’ This new redefinition of the critical load appears 
to be a necessary consequence of accepting the tangent modulus 
stress distribution model as valid for the plastic perfect column. 
The tangent modulus model requires that during buckling the 
load be permitted to increase simultaneously with bending. If 
this model were applied to the elastic perfect column, however, 
its stability limit could not be determined, since the commonly 
accepted definition of the critical load is the load at which an 
exchange of equilibrium configurations occurs with no increase 
in axial load. Thus, it appears necessary to examine critically 
column behavior in both the elastic and plastic range to determine 
whether plastic column fits within the framework 
accepted for the elastic column or if a new definition is required. 
It is the purpose of this note to indicate that both elastic and 
plastic column behavior fit within the same framework, permitting 
a consistent treatment of the column problem throughout the 


behavior 


entire stress range. 


Evastic COLUMN 


The stability problem of the perfect elastic column is charac- 
terized by the following features: It is an eigenvalue problem, the 
solution of which determines the point at which an exchange of 
stable equilibrium configurations occurs. The important detail 
to note is that the equilibrium change occurs at a fixed value of 
axial load. Since the solution is for one point only, the value of 
central deflection remains indeterminate. 
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The problem of the initially imperfect column is not one of 
stability at all, since only the initially bent form is stable. The 
yalue of central deflection is known for all values of load, and it is 
found that, in a region below the buckling load of the perfect 
elastic column, the central deflection of the imperfect column 
grows rapidly. To define the buckling load for the latter if it 
were a perfect column, the Southwell method may be used. 

For an elastic column the postcritical behavior of the perfect 
and imperfect column is essentially the same in the region where 
bending predominates. In this region, the failing load of the 
imperfect column is little influenced by small initial imperfections 
approximately coincides with the Euler load for the perfect 


and 
Thus, the Euler Formula can be 


column within close limits. 
used to predict the failing load of the imperfect column for en- 


gineering purposes. 


PLASTIC COLUMN 


The stability theory for the plastic column is considerably 


postulate various models of stress distribution at the critical 
load to determine the bending stiffness in the plastic region. 
In particular, since the behavior in the region of the critical load 
vitally affects the failing load, which was not the case for the 
elastic column, the buckling model assumes great importance in 
the case of plastic buckling. 

To solve the stability problem for the plastic column, the von 
Karman strain reversal model must be considered valid, since it 
is the only one at which an exchange of equilibrium configurations 
occurs at the same value of load. ‘his model results in a critical 
load that is consistent with the stability definition accepted for 
the elastic column. 

For the imperfect column, which is not a stability problem, it 
appears reasonable to assume a model such as the tangent modu- 
lus model. It is to be noted that the tangent modulus model 
cannot be used for the stability problem, since this model re- 
quires an increase in load as bending becomes large. Thus, the 
conclusion is drawn that the strain reversal model is valid for the 
perfect column, whereas the tangent modulus model may be 
valid for the imperfect column. The critical load for an im- 
perfect plastic column may be determined by use of the modified 
Scuthwell method as proposed by Wang. 

To determine the failing load for a plastic column, it is neces- 
sary to construct a model of the failure mechanism. It appears 
reasonable to hypothesize that failure occurs when, or shortly 
after, the condition of strain reversal is reached. Thus, for the 
perfect plastic column, the buckling load and failing loads should 
coincide. The failing load of the imperfect column, however, 
is governed by the magnitude of the initial imperfections, since 
strain reversal occurs in a region where axial loading and bending 
proceed simultaneously. Thus, we can expect the failing load 
of the plastic imperfect column to vary within rather wide limits, 
depending upon the initial imperfection with the reduced modulus 
load as the upper limit for the case of vanishing imperfections. 


Some Relations Between Harmonic and 
Transient Loading of Airfoils 


John W. Miles 
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um WELL-KNOWN DUALITY between the responses of a linear 
system to harmonic and step function excitations and the 
associated duality between the Fourier and Duhamel methods of 
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superposition, respectively, can be used to estaklish some rather 
useful theorems relative to the transient loading of a thin air- 
foil. Let fo, 1, 2 denote (operationally) the responses to the nor- 
mal (ai the airfoil) velocity distributions wo, 1, » in accordance with 


(1) 


Wo (x, y) = a(x, y) > fo La(x, y)} 


= 9 
_ «) 


a(x, y) exp (ikt) > filk, a(x, y)} exp (ikt)  ( 
(3) 


wi(x, y) 


w(x. ¥) a(x, y)1[t — to(x, y)] — felt, to(x, y), a(x, y) }1(d) 


where ¢ is a dimensionless time obtained by multiplying the true 
time by the quotient of a characteristic velocity and a charac- 
teristic length, & is a dimensionless frequency obtained by dividing 
the true (angular) frequency by this same quotient, (x, y) are 
coordinates in the plane of the airfoil, and f is the time at which 
the normal velocity arises at (x, y). We assume fp to be positive, 
definite and bounded. Then, it follows directly from the afore- 


mentioned duality that 


lim felt, to, a} = fil0, a} = fola} (4) 
t— « 
S210, to, a} = lim f,lk, exp (—ikty)a} (5) 


k—> 


In connection with Eq. (5), we remark that linearized theory 
is not reliable for large k, since terms of order (k.W6)? are neg- 
lected, 6 being the airfoil thickness parameter.! 


felt, to, a} dt 


wre re) 
I, [feit, bh, aj — fol a} ] dt = ES fi {k, a} 


= (6a) 


lim (ik)~! f; {k, exp (—ikto)a} 


k—o 
| — fo itoa} 
k=0 


(6b) 


The integral in Eq. (6a) exists only if fo'a} = 0; otherwise, 
Eq. (6b) is to be used. In connection with Eq. (6), we remark 
that the time required to attain the steady state after a suddenly 
imposed disturbance is M(M — 1)~! for a supersonic airfoil if 
the maximum chord is taken as the characteristic length. 

Eq. (5) may be used to establish the initial lift coefficient (re 
ferred to the product of dynamic pressure and wing area) of an 
airfoil under transient loading as 

(C1), 00 = (4/M) (a), <9 (7) 
where is the free-stream Mach Number, and a is the average 
value of a over that portion of the airfoil for which tj = 0. This 
is a generalization of a result obtained by Heaslet and Lomax? and 
is valid at both subsonic and supersonic speeds. Also, the initial 
lift acts through the centroid of the initial distribution of a. 

A result that is useful in connection with Eq. (6) is® 
gol a} +ikM?2(M? — 1)-! X 

[pol xa} — xdola}] + O(R?) 


1 { k, a} 
(8) 
where ¢ is the velocity potential for an airfoil whose trailing edge 
is everywhere supersonic. 

If the results of Eqs. (6) and (8) are applied to the penetration 
of a uniform, sharp-edged gust by a supersonic, rectangular air- 
foil of effective aspect A, we find 


1 M 3A —2 
C, = [: ~ 3 | + \e - "Ye (9) 
Cr, = 4(M? — 1)~'7? [1 — (24)~!] (10) 
a -( M ) 6A? — 6A +1 | sal 
: M + 1/| 9(M? — 1)'*A(2A — 1) 
1 = (M?—1)'7A.R.>1 (12) 


where C,, and Cy, are the average (incremental) lift and mo- 
ment (with respect to the final center of pressure) coefficients over 
the transient period, C,, is the final lift coefficient (increment), 
and the chord is chosen as the characteristic length. Similarly, 
for a delta airfoil of apex angle 2y entering a sharp-edged gust. 


we find 
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(2M? + 1) — 2(fi/fo) . 
| Ci, 
3M(M + 1) 


Cr. © 4 — 1) Se 


a= § 


9 


: (2M? + 1)fo — 2f, 
Cu, 7” " ite 
7 9M(M + 1) (M? — 1) 


0 = wa/2E'(c), o < (16a) 


=1, ¢62>1 (16b) 
= (3/4) o(1 — o?)[(1 — 207)E'(c) + o?K'(a)|~}, 


a0 <1 (17a) 


Py 23 (17b) 


o = (M* — 1)'“ tan y (18) 
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where K’(c) and E’(a) are complementary, complete, elliptig 
integrals of the first and second kinds, respectively. 

Similar results may be established concerning the TeSponse of 
an airfoil to an impulsive (delta function) gust, since this re. 
sponse may be represented as the time derivative of the Tesponse 
toa uniform, sharp-edged gust. 

All of the foregoing results are, of course, subject to the usual 
restrictions imposed in the linearization of the hydrodynamical 
equations. 
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The Quasi-Steady Air Forces for Use in 


Low-Frequency Stability Calculations 


(Continued from page 608) 


The Laplace transform for F(7) is readily shown to be, 
using the notation of Carslaw and Jaeger," 


F(p) = a(p) + jr “\P) _ Put?) 
r=“) 3#=2) 

The steady-state air forces, which proceed according 
to e”’, are evidently given by the portion of the inverse 
transform of F(p) which corresponds to the pole p = 


jr—iee., 


(A4) 


F,(r, X) = jr[w(jr)e" (A5) 


Eq. (A5) obviously holds equally well for both real and 
complex values of the parameter \ and thus leads to 
the desired conclusion: The circulatory steady-state 
air forces for nonharmonic motions of the kind contem- 
plated here can be derived simply by replacing the pa- 
rameter k by in the expressions for the harmonic air 
forces. 

Regarding the noncirculatory air forces, it can be 
shown without difficulty that the previous conclusion 
continues to hold in these cases as well. 


It is to be noted that a stability study based on the 
steady-state air forces, which is the conventional ap 
proach to the problem, does not necessarily describe 
the stability of the system immediately after a dis 
turbance during the period when aerodynamic starting 
transients are still effective. The starting transients, 
which are defined by the remaining poles in Eq. (Ad), 
are all heavily damped aperiodic forces; hence, their 
presence will only affect the stability of aircraft with 
heavily damped stability modes, in cases where the 
stability motions die away in times commensurate with 
the decay times of the transients. The heavy mode) 
damping required for this situation to exist causes the 
question to be of academic interest only, and the use 
of the steady-state air forces alone for stability analyses 
is thus justified. 


It is of further interest to note that Eq. (A5) implies 
a transform relationship between x(r) and F,(7, d). 
This matter has been fully explored by R. T. Jones and 
W. R. Sears, as well as other investigators. 








